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The  theory  of  time  scales  was  introduced  by  Stefan  Hilger  in  his  1988  PhD  dis¬ 
sertation,  [18].  The  study  of  dynamic  equations  on  time  scales  unifies  and  extends 
the  fields  of  differential  and  difference  equations,  highlighting  the  similarities  and 
providing  insight  into  some  of  the  differences. 

In  his  dissertation,  [18],  Hilger  introduced  the  notion  of  the  “delta-derivative”  on  a 
time  scale.  An  analogous  concept,  the  “nabla-derivative”  was  developed  and  explored 
by  Ferhan  Atici  and  Gusein  Guseinov  in  [4].  It  is  interesting  to  look  at  what  happens 
when  these  two  kinds  of  derivatives  are  present  in  the  same  equation.  The  interaction 
between  them  3delds  some  fascinating  behavior,  which  in  some  cases  is  “cleaner”  than 
the  behavior  found  with  only  one  type  of  derivative. 

In  Chapter  2,  we  examine  the  second-order,  self-adjoint  dynamic  equation  which 
contains  both  delta-  and  nabla-derivatives.  We  develop  a  reduction  of  order  theorem, 
explore  oscillation  and  disconjugacy,  and  look  at  Riccati  techniques.  The  material  in 
Chapter  2  has  been  previously  published  in  [23]  and  [22]. 

In  Chapter  3,  we  look  at  solution  techniques  for  linear  dynamic  equations  which 
can  be  written  in  a  factored  form.  We  develop  complete  results  in  the  case  where  the 
equation  contains  only  one  kind  of  derivative.  We  briefly  discuss  the  mixed  derivative 
case  deferring  further  consideration  to  later  work.  The  material  in  Chapter  3  has 
been  previously  published  in  [21]. 

In  the  final  chapter.  Chapter  4,  we  return  to  the  self-adjoint  equation.  Here, 


we  consider  the  matrix  form  of  the  equation.  As  in  the  scalar  case,  we  develop 
a  reduction  of  order  theorem  and  explore  Riccati  techniques,  culminating  with  the 
proof  of  Jacobi’s  Condition. 

Throughout  much  of  this  dissertation,  the  interaction  between  the  delta-  and 
nabla-derivatives  plays  a  key  role.  In  many  cases,  it  is  rather  startling  to  see  how  all 
of  the  pieces  fit  together.  It  is  our  hope  that  this  work  will  inspire  further  exploration 
in  this  area. 
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Chapter  1 
Introduction 


The  theory  of  time  scales  was  introduced  by  Stefan  Hilger  in  his  1988  PhD  disser¬ 
tation,  [18].  Scholars  in  the  fields  of  differential  equations  and  difference  equations 
have  long  been  aware  of  the  startling  similarities  and  intriguing  differences  between 
the  two  fields.  The  study  of  dynamic  equations  on  time  scales  unifies  and  extends 
the  fields  of  differential  and  difference  equations,  highlighting  the  similarities  and 
providing  insight  into  some  of  the  diflPerences. 

In  his  dissertation,  [18],  Hilger  introduced  the  notion  of  the  “delta-derivative”  (or 
A-derivative)  on  a  time  scale.  This  derivative  is  a  generalization  of  both  the  usual 
derivative  firom  differential  equations  and  the  forward  difference  operator  from  dif¬ 
ference  equations.  An  analogous  concept,  the  “nabla-derivative”  (or  V-derivative) 
was  developed  and  explored  by  Ferhan  Atici  and  Gusein  Guseinov  in  [4].  The  nabla- 
derivative  is  a  generalization  of  the  usual  derivative  and  the  backward  difference  oper¬ 
ator.  Study  of  the  nabla-derivative  alone  is  not  particularly  revealing,  since  the  results 
are  usually  directly  analogous  to  results  developed  for  the  delta-derivative.  What  is 
interesting,  however,  is  to  look  at  what  happens  when  these  two  kinds  of  derivatives 
are  present  in  the  same  equation.  The  interaction  between  these  two  derivatives  yields 
some  fascinating  behavior,  which  in  some  cases  is  “cleaner”  than  the  behavior  found 
when  only  one  type  of  derivative  is  considered. 

In  Chapter  1  of  this  work,  we  provide  a  brief  introduction  to  the  calculus  on  a 
time  scale,  including  both  the  delta  and  nabla  derivatives,  for  the  reader  who  may 
be  unfainiliar  with  the  field.  We  also  summarize  some  of  the  key  results  concerning 
“generalized  exponential  functions” ,  which  play  a  key  role  in  the  study  of  dynamic 
equations  on  time  scales.  The  material  in  this  chapter  is  taken  from  existing  sources, 
primarily  [6,  7,  4,  3],  which,  in  turn,  rely  on  previously  published  works  by  Hilger 
and  others,  and  we  refer  the  interested  reader  to  those  sources  for  a  more  complete 
introduction. 

In  Chapter  2,  we  examine  the  second-order,  self-adjoint  dynamic  equation 

[p(t)x^(t)]^ -4- q(t)x(t)  =  0 
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on  a  time  scale.  Although  the  similar  equation, 

\p{t)x^]^  +  q{t)x'^  =  0, 

has  been  studied  extensively,  (see,  for  example  [9,  10,  11,  5,  12,  15,  13,  17]),  little 
work  has  been  done  on  this  equation,  which  combines  both  the  delta  and  nabla 
derivatives.  We  begin  by  establishing  several  results  concerning  the  interaction  of 
these  two  derivatives.  Also  included  in  the  first  section  is  a  theorem  which  shows  that 
under  certain  conditions,  the  generalized  exponential  functions  ep{t,to)  and  er{t,to) 
can  be  related  to  one  another.  We  next  look  at  three  second-order  linear  equations 
and  demonstrate  that  they  can  be  written  in  self-adjoint  form.  The  first  results  which 
are  directly  related  to  the  self-adjoint  equation  are  contained  in  Section  2.3,  which 
culminates  with  a  reduction  of  order  theorem.  We  turn  our  attention  to  oscillation  and 
disconjugacy  in  Section  2.4,  where  we  establish  an  analogue  of  the  Sturm  Separation 
Theorem,  and,  via  the  Polya  and  Trench  factorizations  demonstrate  the  existence  of 
recessive  and  dominant  solutions  of  the  self-adjoint  equation.  The  final  section  of  the 
chapter.  Section  2.5,  discusses  Riccati  techniques  as  they  relate  to  the  self-adjoint 
equation.  The  material  in  Chapter  2  has  been  previously  published  in  [23]  and  [22]. 

We  alter  our  focus  slightly  in  Chapter  3.  There,  we  look  at  solution  techniques 
for  linear  dynamic  equations  on  time  scales  which  can  be  written  in  a  factored  form. 
As  special  cases,  we  obtain  solution  techniques  for  constant  coefficient  dynamic  equa¬ 
tions  and  for  Euler-Cauchy  dynamic  equations.  These  special  cases  are  consistent 
with  results  developed  by  Akin-Bohner  and  Bohner  in  [7,  2].  We  consider  equations 
containing  either  delta-derivatives  or  nabla^derivatives,  and  develop  complete  solu¬ 
tion  techniques  in  these  cases.  We  briefly  discuss  second-order  equations  containing 
both  delta  and  nabla  derivatives,  and  defer  consideration  of  higher  order  equations 
with  a  mixture  of  the  two  kinds  of  derivatives  to  later  work.  We  conclude  Chapter 
3  by  examining  another  equation  which  cannot  itself  be  written  in  factored  form, 
but  is  equivalent  to  one  which  can.  The  material  in  Chapter  3  has  been  previously 
published  in  [21]. 

In  the  final  chapter.  Chapter  4  we  return  to  the  self-adjoint  equation  with  mixed 
derivatives.  In  this  chapter,  we  consider  the  matrix  form  of  the  equation.  As  in  the 
scalar  case,  the  related  scalar  equation  which  contains  only  deltarderivatives  has  been 
studied  extensively  (see,  for  example  [1, 14, 15, 13]).Many  of  the  results  in  this  chapter 
are  analogous  to  the  results  obtained  in  Chapter  2  for  the  scalar  case.  We  begin  by 
establishing  the  Lagrange  Identity  for  this  equation,  which  leads  to  the  development 
of  Abel’s  Formula  and  a  reduction  of  order  theorem.  From  there,  we  move  into  a 
section  on  Riccati  techniques,  where  we  develop  the  relationship  between  the  self- 
adjoint  matrix  equation,  and  it’s  related  Riccati  equation.  Chapter  4  culminates 
with  the  proof  of  Jacobi’s  Condition,  which  is  much  more  difficult  to  establish  in  the 
matrix  case  than  the  scalar  case. 

Throughout  much  of  this  dissertation,  the  interaction  between  the  delta  and  nabla 
derivatives  plays  a  key  role.  Although  there  are  some  theorems  where  the  interaction 
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does  not  have  much  impact,  in  most  cases,  it  is  rather  startling  to  see  how  all  of  the 
pieces  fit  together.  It  is  our  hope  that  this  work  will  inspire  further  exploration  in 
this  area. 


1.1  The  Calculus  on  a  Time  Scale 

A  time  scale  is  simply  a  closed  subset  of  R.  Throughout  this  work,  we  will  use 
the  symbol  T  to  represent  a  time  scale.  A  given  time  scale  is  assumed  to  have  the 
topology  which  it  inherits  as  a  subspace  of  R  with  the  usual  topology.  The  notation 
[a,  b]  is  understood  to  mean  the  real  interval  [a,  b]  intersected  with  T.  Open  and  half- 
open  intervals  are  understood  similarly.  We  begin  by  introducing  the  forward  and 
backward  jump  operators  on  T. 

Definition  1.  Let  T  be  a  time  scale.  For  t  e  T,  we  define  the  forward  jump  operator, 
£r  :  T  — »■  T  by 

a{t)  :=  inf{s  €  T  :  s  >  t}. 

If  {s  G  T  :  s  >  t}  =  0,  (i.e.  if  t  =  maxT),  we  take  a{t)  =  t. 

Similarly,  we  define  the  backward  jump  operator,  p(t)  ;  T  — »  T  by 

p{t)  :=  sup{s  eT  :  s  <t}, 

and,  if  t  =  min  T,  we  take  p{t)  =  t. 

If  /  :  T  ^  R,  then  the  notation  is  understood  to  mean  f{a{t)),  and  /^(t)  is 
understood  to  mean  f{p{t)). 

Points  in  T  are  classified  as  follows:  If  a{t)  >  t,  we  say  t  is  right-scattered.  If 
t  <  supT  and  a{t)  =  t,  we  say  t  is  right-dense.  If  p(t)  <  t,  we  say  t  is  left-scattered, 
and  if  t  >  inf  T  and  p{t)  =  t,  we  say  t  is  left-dense.  Points  that  axe  both  right  and  left 
scattered  are  called  isolated,  and  points  that  are  both  right  and  left  dense  axe  called 
dense. 

Definition  2.  A  function  /  :  T  ^  R  is  said  to  be  right-dense-continuous  ox  rd- 
continuous  provided  /  is  continuous  at  all  right-dense  points  in  T,  and  provided 
lims_t  /(s)  exists  and  is  finite  at  all  left-dense  t  eT. 

We  say  /  is  left-dense-continuous  or  Id-continuous  provided  /  is  continuous  at  all 
left-dense  points  in  T,  and  provided  lims^t  /(s)  exists  and  is  finite  at  all  right-dense 
t  eT. 

Definition  3.  The  graininess  function,  p,{t) ,  is  defined  by 


p{t)  :=  <7{t)  —  t. 

The  backward  graininess  function,  u(t)  is  defined  by 

u{t):=t- p{t). 
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Remark  4.  Recently,  there  has  been  some  disagreement  regarding  the  most  appro¬ 
priate  definition  of  u{t).  In  this  work,  we  retain  the  original  definition,  which  is 
consistent  with  previously  published  literature  on  V-derivatives.  It  is  inconsistent, 
however  with  the  current  work  on  a-derivatives.  When  working  with  a-derivatives, 
the  a-graininess,  fia  is  defined  to  be  //q,  :=  a(t)  —  t.  When  a(t)  =  p{t),  then,  we 
would  have  Hp  :=  p{t)—t  =  —u{t).  This  inconsistency  is  unfortunate,  but  we  feel  it  is 
more  important  that  we  remain  consistent  with  the  way  v>{t)  was  defined  in  previously 
published  work.  To  minimize  confusion,  we  recommend  the  notation  Pp{t)  =  p{t)  —  t 
be  used  in  work  that  is  to  be  interpreted  in  the  more  general  a-derivative  setting. 

Definition  5.  The  set  is  defined  as  follows.  If  T  has  a  left-scattered  maximum, 
M,  then  T'‘  :=  T  \  {M}.  Otherwise,  T''  =  T.  Similarly,  if  T  has  a  right-scattered 
minimum,  m,  we  define  the  set  :=  T  \  {m}.  Otherwise,  =  T. 

Definition  6.  Assume  /  :  T  — >  M,  and  let  t  Then  we  define  the  delta- derivative 
of  f  at  t,  denoted  f^{t)  to  be  the  number  (provided  that  it  exists)  with  the  property 
that  given  any  e  >  0,  there  is  a  neighborhood,  U  of  t  such  that 

l[/(<^(<))  -  /(«)]  -  fHtMt)  -  s)|  <  -  s\ 


for  all  s  G  {/. 

If  f^(t)  exists  for  all  t  G  T*',  then  we  say  /  is  delta-differentiable,  and  we  call 
/^  :  T**  — »•  M  the  delta-derivative  of  /  on  T'‘.  It  is  straightforward  to  show  that  for 
t  G  is  well  defined.  Note,  however  that  if  T  has  a  left-scattered  maximum, 

M,  then  f^{M)  is  not  uniquely  determined.  This  is  precisely  why  we  eliminate  this 
type  of  point  in  our  set  T'*. 

The  nabla-derivative  is  defined  in  similar  fashion  to  the  delta-derivative. 

Definition  7.  Assume  /  :  T  — »  M,  and  let  t  G  T^.  Then  we  define  the  nabla-derivative 
of  f  at  t,  denoted  /^(t)  to  be  the  number  (provided  that  it  exists)  with  the  property 
that  given  any  £  >  0,  there  is  a  neighborhood,  U  of  t  such  that 

,  |[/w<)) - /(s)l - /’(*) W») -4< - s| 

for  all  s  E  U. 

If  /^(t)  exists  for  all  t  G  T,,,,  then  we  say  /  is  nablardifferentiable,  and  we  call 
/^  :  Tk  — K  the  nabla-derivative  of  /  on  T^. 

Note  that  in  the  case  T  =  R,  both  the  delta  and  nabla  derivatives  are  simply  the 
usual  derivative.  If  T  =  Z,  the  delta  derivative  is  the  forward  difference  operator, 
while  the  nabla  derivative  is  the  backward  difference  operator.  Both  the  delta  and 
nabla  derivatives  possess  many  useful  properties.  We  summarize  some  of  them  here. 

Theorem  8.  Assume  /  ;  T  — >  M  is  o  funetion  and  let  t  G  T'^.  Then 

(V  U  f  ^  delta-differentiable  at  t,  then  f  is  continuous  at  t. 
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(ii)  If  f  is  continuous  at  t  and  t  is  right-scattered,  then  f  is  delta-differentiable  at 
t,  and 

(Hi)  Iftis  right-dense,  then  f  is  differentiable  at  t  if  and  only  if  the  limit 

S^t  t  —  $ 

exists  as  a  finite  number.  In  this  case 


s — >t  i)  —  s 


(iv)  If  f  is  delta- differentiable  attj  then 

rw  = /w +M^)/^(0- 

Theorem  9.  Assume  5  :  T  — »  R  is  o  function,  and  let  i  €  T^.  Then 

(i)  If  g  is  nabla- differentiable  at  t,  then  g  is  continuous  att. 

(ii)  If  g  is  continuous  at  t  and  t  is  left-scattered,  then  g  is  nabla-differentiable  at  t, 
and 

ait)  - 

(Hi)  Iftis  left-dense,  then  g  is  differentiable  at  t  if  and  only  if  the  limit 

s-*t  t  —  S 

exists  as  a  finite  number.  In  this  case 

^  ^  ^  s->t  t-  S 


(iv)  If  g  is  nabla-differentiable  att,  then 

gP{t)  =  g{t)  -  u{t)g^(t). 

Looking  at  properties  (ii)  and  (Hi)  in  each  of  the  theorems  above  gives  us  a  more 
intuitive  understanding  of  these  derivatives  than  can  be  gained  via  the  definition 
alone.  In  the  case  of  the  delta-derivative,  if  f  €  T  is  right-dense,  then  the  delta- 
derivative  behaves  in  much  the  same  way  as  the  usual  derivative.  It  can  be  regarded 


6 


as  the  slope  of  the  tangent  line  to  the  function  at  t,  although  if  t  is  both  right-dense 
and  left-scattered,  the  limit  is  a  one-sided  limit.  On  the  other  hand,  if  t  is  right- 
scattered,  then  the  f^{t)  is  the  slope  of  the  line  segment  connecting  f{t)  and  f{(r{t)). 
In  this  case,  it  makes  no  difference  what  the  behavior  of  the  function  is  like  to  the 
left  of  t,  beyond  the  requirement  that  /  be  continuous  at  t.  So,  the  delta  derivative 
provides  a  mixture  of  the  discrete  behavior  of  the  forward  difference  operator,  and  the 
continuous  behavior  of  the  usual  derivative.  The  nabla  derivative  can  be  interpreted 
similarly,  but  in  that  case  the  focus  is  on  what  happens  to  the  left  of  a  given  point, 
rather  than  the  right. 

It  is  easy  to  show  that  both  the  delta  and  nabla  derivatives  are  linear.  There  are 
also  product  and  quotient  rules  for  both  derivatives. 

Theorem  10.  Assume  f,g:T-^'Rare  delta-differentiahle  at  t  G  T''.  Then 
(i)  The  product  fg  is  delta- differentiable  at  t,  with 

=  fHt)9{t)  +  nt)g^(t)  =  +  /(()5^(*). 

(ii)  If  g{t)g'^{t)  7^  0,  then  ^  is  delta-differentiable  at  t,  with 

9{t)9''{t) 

Theorem  11.  Assume  are  nabla- differentiable  att  ^  Then 

(i)  The  product  fg  is  nabla-differentiable  at  t,  with 

ifgfit)  =  fmt)  +  /'{«)«’(*)  =  +  /(«)/(<). 

(ii)  If  g{t)gP{t)  ^  0,  then  I  is  nabla-differentiable  at  t,  rvith 

rmt)  -  f{t)g^{t) 

9{t)9''it) 

Of  course,  the  calculus  on  a  time  scale  would  not  be  complete  without  a  concept 
of  integration  to  complement  the  derivative.  It  can  be  shown  that  if  a  function 
/  :  T  ^  M  is  rd-continuous,  then  it  has  a  delta  antiderivative.  Similarly,  if  5^ ;  T  — »■  K 
is  Id-continuous,  then  it  has  a  nabla  antiderivative.  We  then  define  the  delta  and 
nabla  integrals  of  /  and  g  in  terms  of  these  antiderivatives. 

Definition  12.  Assume  /  :  T  — >  M  is  rd-continuous,  and  let  F  be  a  delta  antideriva¬ 
tive  of  /.  That  is,  suppose  F^{t)  =  f{t)  for  all  t  €  T''.  Then  the  indefinite  delta 
integral  of  /  is  given  by 

J  f{t)At  =  F{t)  +  C. 
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The  delta  Cauchy  integral  is  defined  by 

f{t)  At  =  F{s)  —  F{r)  for  all  r,  s  6  T. 

Definition  13.  Assume  5- :  T  ^  E  is  Id-continuous,  and  let  G  be  a  nabla  antideriva¬ 
tive  of  g.  That  is,  suppose  G'^{t)  =  g{t)  for  all  t  €  T^.  Then  the  indefinite  nabla 
integral  of  g  is  given  by 

J'  gif')  Vt  =  G(t)  -b  C. 

The  nabla  Cauchy  integral  is  defined  by 

g{t)  Vt  =  G{s)  —  G{r)  for  all  r,  sgT. 

These  definitions  can  be  generalized  to  apply  to  a  more  broad  category  of  functions 
called  regulated  functions.  See,  for  example,  [6,  Section  1.4].  Since  our  applications 
of  integration  usually  fall  within  the  more  restrictive  definition,  however,  we  will  not 
dwell  on  the  somewhat  technical  details  of  this  more  general  definition. 

As  with  the  derivatives,  both  integrals  are  linear.  Additional  properties  are  given 
in  the  following  theorems. 

Theorem  14.  If  a,  b,  c  eT,  and  f,  g  are  rd-continuous,  then 

(i) 

(ii)  At  =  £f(t)  At  +  m  At. 

M  At  =  {fg){b)  -  (fg){a)  -  At. 

H  Jl  m9^(t)  Ai  =  (fg)(b)  -  (/s)(a)  -  r  (()<(()  At. 

Theorem  15.  If  6,  c  G  T,  and  /,  g  are  Id- continuous ,  then 

(i)  m  = -  s:  m  vt. 

(ii)  i‘  m  vt = /;  m  vt + m  v*. 

(iU)  Slnt)g'’(t)  Vi  =  (fg)(b)  -  (fg)(a)  -  f^r{t)g(t)  Vi. 

W  /(*)«’(*)  Vi  =  (fg){b)  -  {fg){a)  -  Vi. 
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1.2  Generalized  Exponential  Functions 

In  the  study  of  differential  equations,  the  properties  of  the  exponential  function 
are  crucial  to  many  of  the  standard  results.  In  difference  equations,  there  is  a  similar 
reliance  on  functions  of  the  form  a‘.  In  this  section,  we  look  at  generalized  exponential 
functions  on  time  scales,  which  will  play  this  key  role  in  the  time  scales  context. 

Definition  16.  We  say  that  a  function  r  :  T  — R  is  regressive  if 

1  +  p{t)r{t)  ^  0 


for  all  t  €  T*'. 

We  further  define  the  sets  Tt  and  by 

7?.  :=  {r  :  T  — s-  R  I  r  is  rd-continuous  and  regressive} 

and 

TZ'^  :={r  eTZ\  1  +  p{t)r{t)  >  0  for  all  teT} 

We  now  define  the  “circle  plus”  addition,  ©,  on  TZ.  If  p,  q  E  TZ,  then 

p®q:=p  +  q  +  ppq, 

where  yu  is  the  graininess  function.  It  can  be  shown  that  TZ  is  an  Abelian  group 
under  this  circle  plus  addition.  For  this  reason,  we  refer  to  7^  as  the  regressive  group. 
The  family  TZ'^  of  positively  regressive  functions  is  a  subgroup  of  TZ  under  circle  plus 
addition. 

If  p  G  7^,  then  the  additive  inverse  of  p  under  circle  plus  addition  is  denoted  ©p, 
and  is  given  by 


For  p,  q  eTZ,  we  write  pQ  q  to  denote  p  ©  (Qq). 

Ifr  eTZ,  then  the  initial  value  problem 

=  r{t)y,  y{to)  =  1 

has  a  unique  solution.  We  denote  this  solution  by 


6r(’)  fo)) 

and  call  it  the  generalized  (delta)  exponential  function.  It  can  be  shown  that  if 
r  E  TZ^,  then  er{-,to)  >  0  for  all  t  eT. 

Some  of  the  key  properties  of  the  generalized  (delta)  exponential  function  are 
given  here. 


Theorem  17.  Ifp,  q  ETZ,  then 
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(i)  eo{t,s)  =  1  and  ep{t,t)  =  1; 

(a)  ep{a{t),  5)  =  (1  +  fi{t)p{t))ep{t,  s); 

ep(t,s)  ~  ®©p(^> 

(iv)  ep(t,s)  =  =  eep(s,t); 

(v)  ep{t,s)ep{s,r)  =  ep{t,r); 

(vi)  Cp(tf  s)Cq(t,  s)  6p0g(i,  5); 

W  =  epe9(^>s); 

(^)  =  -iffe 

A  generalized  nabla  exponential  function  can  be  developed  in  an  analogous  fashion. 
We  state  the  relevant  definitions  and  properties  here. 

Definition  18.  We  say  that  a  function  r  :  T  — ^  M  is  u-regressive  if 

1  —  v{t)r(t)  7^  0 


for  all  t  G  Tr. 

We  further  define  the  sets  TZi^  and  by 

7?,^  :=  {r  :  T  — »•  M  I  r  is  Id-continuous  and  u  —  regressive} 


and 

TZp  :=  [r  eTlv  \  l-  v{t)r{t)  >  0  for  all  t  eT} 

The  set  71^  is  an  Abelian  group  under  the  operation  0^,  which  is  defined  by 


p®uq-=p  +  q-  J^pq, 


for  p,  q  E  TZv  The  family  TZ'^  is  a  subgroup. 

If  p  e  1Zv,  then  the  additive  inverse  of  p  under  ®u  is  denoted  Qup,  and  is  given  by 


©I/P  — 


-P 

1  —  up 


For  p,  9  G  72.^,  we  write  pQ^q  to  denote  p  0^/  (©i/?)- 
Ur  eTZi,,  then  the  initial  value  problem 


=  r{t)y,  yito)  =  1 
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has  a  unique  solution.  We  denote  this  solution  by 

and  call  it  the  generalized  nabla  exponential  function.  It  can  be  shown  that  if  r  € 
then  er{-,to)  >  0  for  all  t  G  T. 

Some  of  the  key  properties  of  the  generalized  nabla  exponential  function  are  given 
here. 

Theorem  19.  If  p,  q  E  TZv,  then 
(i)  eo{t,  s)  =  1  and  ep{t,t)  =  1; 

(ii)  ep{p{t),s)  =  (1  -  u{t)p{t))ep{t,s); 

ep(t,s) 

(iv)  Cpit^s)  = 

(v)  6p(f,  5)6p(5,  r)  =  6p(f,  r); 

(vi)  6p(i,  5)6g(t,  5)  ^); 

(j;fe)  =  - 


P(t) 
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Chapter  2 

A  Second-Order  Self-Adjoint 
Dynamic  Equation  on  a  Time  Scale 

2.1  Preliminary  Results 

This  chapter  is  concerned  with  the  study  of  the  second-order  self-adjoint  equation 
\p{t)x^]^ +q{t)x  =  0  on  a  time  scale.  Since  the  equation  we  are  interested  in  contains 
both  A  and  V-derivatives,  we  will  want  to  know  how  these  two  different  derivatives 
interact.  Toward  this  end,  we  will  need  an  analog  of  L’Hopital’s  Rule.  A  version  of 
this  crucial  theorem  for  A-derivatives,  is  contained  in  [6],  although  it  is  presented 
here  in  a  slightly  different  form  (Theorem  22).  We  then  develop  L’Hopital’s  rule  for 
V-derivatives  (Theorem  23). 

We  may  want  to  employ  L’Hopital’s  Rule  to  evaluate  a  limit  as  t  — »  ±oo,  so  we 
make  the  following  definitions. 

Definition  20.  Let  e  >  0.  If  T  is  unbounded  above,  we  define  a  left  neighborhood  of 
oo,  which  we  denote  by  Le(oo),  by 

Le(oo)  =  €  T  :  t  >  ^ 

Similarly,  if  T  is  unbounded  below,  we  define  a  right  neighborhood  of  — oo,  denoted 
Rg(-oo)  by 

Re(-oo)  =  G  T  :  i  <  -^ 

We  next  define  right  and  left  neighborhoods  for  points  in  T. 

Definition  21.  Let  e  >  0.  For  any  right-dense  to  G  T,  define  a  right  neighborhood  of 
to,  denoted  Rg{to),  by 


R^ifo')  ■ —  •{tGT;0<t  —  to  c}. 
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Similarly,  for  any  left-dense  to  €  T  ,  define  a  left  neighborhood  of  to,  denoted  Le(to)) 
by 

Lg{to)  :={teT:0<to-t<s}. 

Theorem  22  (L’Hopital’s  Rule  for  A-derivatives).  Assume  f  and  g  are  A- 

differentiable  on  T,  and  let  to  G  TU  {oo}.  If  to  G  T,  assume  to  is  left-dense.  Further¬ 
more,  assume 

lina  f{t)  =  Una  g{t)  =  0, 

t-^tQ 

and  suppose  there  exists  e  >  0  with 

g{t)g^{t)  <  0  for  all  t  e  Le{to)- 


Then  we  have 


lim  inf  <  lim  inf  ^  <  lim  sup  ^ 
t^t-  g^t)  g(t)  t  git) 


<  lim  sup 
t^to 


fHt) 

g^ty 


Theorem  23  (L’Hopital’s  Rule  for  V-derivatives).  Assume  f  and  g  are  V- 
differentiable  on  T  and  let  to  €  T  U  {— oo}.  1/ to  G  T,  assume  to  is  right-dense. 
Furthermore,  assume 

lim  fit)  =  lina  git)  =  0, 

t—^tQ 

and  suppose  there  exists  e  >  0  with 


git)g'^it)  >  0  for  all  t  6  i?e(to). 


Then 


lim  inf  <  lim  inf  <  limsup  ^  <  lim  sup 

t-tj  g^it)  t-.tt  git)  9it)  9^(^) 


Proof  Without  loss  of  generality,  assume  git)  and  g^it)  are  both  strictly  positive  on 


ile(to).  ^  ^ 

Let  5  e  (0,e],  and  let  a  :=  inf^ea,(fo)  ^  —  sup^eRi(to)  To  complete  the 

proof,  it  suffices  to  show 


a  < 


inf  44  < 


sup 


fi^) 


reRs(to)  gir)  r€R,(to)  gi^) 


<b, 


as  we  may  then  let  5  0  to  obtain  the  desired  result. 

We  must  be  careful  here,  as  either  a  or  6  could  possibly  be  infinite.  Note,  however, 
that  since  g^ir)  >  0  on  Rsito),  we  have  a  <  oo.  Similarly,  b  >  —oo.  So  our  only 
concern  is  if  a  =  — oo  or  b  =  oo.  But,  if  a  =  — oo,  we  have  immediately  that 


fjr) 

r^Mto)  gir)’ 


a  <  inf 
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as  desired,  and  if  6  =  oo,  we  have  immediately  that 

sup  /  ^ 

reRsito)  9v) 

as  desired.  Therefore,  we  may  assume  that  both  a  and  b  are  finite.  Then 

^  /^('^)  ^  ’’’  ^  Rsito), 

and  by  a  theorem  of  Guseinov  and  Kaymakgalan  [16], 

/  ag'^{T)  Vr  <  f  /^(t)  Vr  <  f  bg^ (t)  Vr  for  all  s,t  G  Rs{to),  t  <  s. 

Jt  Jt  Jt 

Integrating,  we  see  that 

ag{s)  -  ag{t)  <  f{s)  -  f{t)  <  bg{s)  -  bg{t)  for  all  s,i  €  Rsih),  t  <  s. 
Letting  t  we  get 

o-9{s)  <  /(s)  <  bg{s)  for  all  s  €  Rs{to), 

and  thus 

a  <  inf  — ^  <  sup  <  b. 

seRs(to)  g{s)  seRi(to)  9\^) 

Then,  by  the  discussion  above,  the  proof  is  complete.  □ 

Remark  24.  Although  these  theorems  are  only  stated  in  terms  of  one-sided  limits, 
analogous  results  can  be  established  if  the  limit  is  taken  firom  the  other  direction.  To 
apply  L’Hopital’s  rule  using  A-derivatives  and  a  right-sided  limit,  to  must  be  right- 
dense  (or  — oo  if  T  is  unbounded  below),  and  99^  must  be  strictly  positive  on  a  right 
neighborhood  of  to-  Similarly,  to  apply  L’Hopital’s  rule  using  V-derivatives  and  a 
left-sided  limit,  to  must  be  left-dense  (or  00  if  T  is  unbounded  above),  and  99^  must 
be  strictly  negative  on  some  left  neighborhood  of  to- 

In  order  to  determine  when  the  two  types  of  derivatives  may  be  interchanged,  we 
need  to  consider  some  of  the  points  in  our  time  scale  separately,  so  let 

A  :=  {t  G  T  I  t  is  left-dense  and  right-scattered},  and  Tyi  :=  T  \  A. 
Similarly,  let 

jB  :=  {t  G  T  I  t  is  right-dense  and  left-scattered},  and  :=T\B. 
Lemma  25.  Ifte  then  (r{p{t))  —  t.  IfteTs,  then  p{cT{t))  =  t. 
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Proof.  We  will  only  prove  the  first  statement.  The  proof  of  the  second  statement  is 
similar.  If  t  G  T^,  then  either  t  is  left-scattered,  or  t  is  both  left-dense  and  right-dense. 
If  t  is  left-scattered,  then  p(t)  is  right-scattered  and  it  is  clear  that  cr{p{t))  =  t.  If  t  is 
both  left-dense  and  right-dense,  then  a(t)  =  t  and  p{t)  =  t.  Hence  cr{p{t))  =  a{t)  =  t. 
In  either  case  we  get  the  desired  result.  □ 

Theorem  26.  If  f  :T  ^  is  A-differentiable  on  T"  and  f^  is  rd-continuous  on 
then  f  is  V -differentiable  on  and 


f 


f^m) 

lims_»t-  /^(s) 


tGTA 
te  A. 


If  9  :  T  — >  M  is  V -differentiable  on  and  is  Id-continuous  on  T^,  then  g  is 
A-differentiable  on  and 
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A 


lims_t+  g'^{s) 


t  G 

te  B. 


Proof  We  will  only  prove  the  first  statement.  The  proof  of  the  second  statement  is 
similar.  First,  assume  t  G  T^.  Then  there  are  two  cases:  Either 


(i)  t  is  left-scattered,  or 


(ii)  t  is  both  left-dense  and  right-dense. 


Case  (i):  Suppose  t  is  left-scattered  and  /  is  A-differentiable  on  T''.  Then  p{t)  is 
right-scattered,  and 

^  cripit))  -  Pit) 


Now,  as  /  is  A-differentiable  on  /  is  continuous  on  T.  Then,  since  t  is  left- 
scattered,  /  is  V-differentiable  at  t,  and  we  see  that 


fHpit))  = 


/(t^(p(0))  -  fipjt)) 
o-(p(0)  -  pit) 
fit)  -  fipji)) 

t  -  pit) 

fit)- 


Case  (ii):  Now,  suppose  t  is  both  left-dense  and  right-dense,  and  /  :  T  ^  M  is  contin¬ 
uous  on  T  and  A-differentiable  at  t.  Since  t  is  right-dense  and  /  is  A-differentiable 
at  t,  we  have  that 

lim 
s— t  —  S 

exists.  But  t  is  left-dense  as  well,  so  this  expression  also  defines  fit),  and  we  see 


fit)  -  fis) 
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that 


fit) 


lim 

S—^t 


m  -  m 

t  —  s 


rw 


So,  we  have  established  the  desired  result  in  the  case  where  t  &  T a- 

Now  suppose  te  A.  Then  t  is  left-dense.  Hence  /^(t)  exists  provided 


u„/(*WW 

s^t  t  —  S 


exists. 

As  t  is  right-scattered,  we  need  only  consider  the  limit  as  s  — s-  <  from  the  left. 
Then  we  apply  L’Hopital’s  rule  [6] ,  differentiating  with  respect  to  s  to  get 


lim 


m  -  m 

t  —  s 


lim  f^{s). 

5— 


Since  we  have  assumed  that  is  rd-continuous,  this  limit  exists.  Hence  /  is 
V-differentiable,  and  -  lims_t-  f^it),  as  desired.  □ 

Corollary  27.  If  to  6  T,  and  /  :  T  M  is  rd-continuous  on  T,  then  f{T)AT  is 
V-differentiable  on  T  and 


fipit)) 

limg^t-  /(s) 


teTA 
te  A. 


IftoE'T,  and  p  :  T  — ♦  M  is  Id-continuous  on  T,  then  g{T)'VT  is  A- differentiable 
on  T  and 

[  3(t)Vt 
.Jto 

The  following  corollary  was  previously  established  by  Atici  and  Guseinov  in  their 
work  [4]. 

Corollary  28.  If  f  :  T  M  is  A-differentiable  on  and  if  f^  is  continuous  on 
T**,  then  f  is  V-differentiable  on  and 

fit)  =  fP{t)  forteT,. 

Ifg  :T  is  V-differentiable  on  T**  and  if  g^  is  continuous  on  Tr,  then  g  is 

A-differentiable  on  and 


\  9i<^it))  te^B 
\  lims_t+  g(s)  te  B. 


g^{t)  =  g^^it)  forteT\ 
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Now,  there  are  a  couple  more  integral  formulas  that  will  be  useful,  the  first  two 
of  which  were  established  in  [6]  and  [4] . 

Lemma  29.  The  following  hold: 

(i)  /(«)  =  ><(*)/(«) 

(«)  A)  M  =  ‘'W/W 

fiii; /,"“>/(«)  Vs  =  M«)rW 

/p(i)  /W 

We  conclude  this  section  by  exploring  the  relationship  between  the  generalized 
exponential  functions  associated  with  the  A  and  V  derivatives. 

Lemma  30.  Let  p  :  T  — M.  Then  p  is  regressive  if  and  only  if  — is  u-regressive, 
and  1  +  ij.{t)p{t)  >  0  for  all  t  e  T  if  and  only  if  1  +  >  0  for  all  t  €  T. 

Similarly,  if  q  :  T  R,  then  q  is  u-regressive  if  and  only  if  —q'^  is  regressive,  and 
1  —  u{t)q{t)  >  0  for  allt  eT  if  and  only  ifl  —  p{t)q^{t)  >  0  for  all  t  €  T. 

Proof  We  will  only  prove  the  first  statement.  The  proof  of  the  second  statement  is 
similar. 

First,  assume  p  is  regressive.  We  then  wish  to  show  that  1  +  u{t)(jf{t))  ^  0. 

Case  1:  Fix  t  e  Then  p{t)  €  T,  and  as  p  is  regressive,  we  have  that 

1  +  p{p{t))p{p{t))  ^  0, 

so,  using  the  definition  of  p,{t) , 

1  +  [<T{p{t))  -  p{t)]p{pit))  ^  0. 

But  t  6  Tyi,  so  a{p{t))  =  t,  and  we  get 

1  +  [t  -  p{t)]fit)  ^  0, 
or 

1  +  v{t)f(t)  ^  0 

as  desired. 

Case  2:  Fix  t  E  A.  Then  t  is  left-dense  and  right-scattered,  so  u(t)  =  0.  Hence 
1  -I-  u{t)p^{t)  =  1-1-  0p'’(i)  =  1^0. 

As  1  +  u{t)pA{t)  ^  0  for  any  t  G  T,  we  see  that  —pf  is  i/-regressive. 
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Conversely,  suppose  — p''  is  i/-regressive.  We  then  wish  to  show  that  i=- 

0. 

Case  1:  Fix  t  e  Tb-  Then  a{t)  G  T,  and,  as  —jf  is  i^-regressive,  we  have  that 


1  +  u{a{t))jf{cr{t))  ^  0, 


so,  using  the  definition  of  ^{t), 

1  +  [a{t)  -  p{a(t))]p{p{a{t)))  +  0. 

But  t  G  Tb,  so  p{o(t))  =  t,  and  we  get 

1  +  [<T(t)  -  t]p{t)  ^  0 
or 

1  +  p{t)p{t)  ^  0 

as  desired. 

Case  2;  Fix  t  G  S.  Then  t  is  right-dense  and  left-scattered,  so  p{t)  =  0.  Hence 

1  -t-  p{t)p{t)  =  1-f  0p{t)  =  1^0. 

As  1  +  p,{t)p{t)  ^  0  for  any  t  G  T,  we  see  that  p  is  regressive. 

To  show  1  -t-  p,{t)p{t)  >  0  for  alH  G  T  if  and  only  if  1  +  >  0  for  all  t  G  T, 

simply  replace  0”  by  ”  >  0”  in  the  preceding  proof.  □ 

Theorem  31  (Equivalence  of  delta  and  nabla  exponential  functions).  Ifp  is 

continuous  and  regressive,  then 

ep(t,  to)  (t,  to)  to). 

l+pPl/ 

If  q  is  continuous  and  u-regressive,  then 

~  6  <i‘’  it,  to)  (t,  to). 

Proof.  We  will  only  prove  the  first  statement.  The  proof  of  the  second  statement 
is  similar.  Suppose  that  p  :  T  ^  M  is  continuous  and  regressive,  then  by  Lemma 
30  we  have  that  —jf  is  i/-regressive.  Furthermore,  since  p  is  continuous,  —pP  is  Id- 
continuous.  Hence  —jP  G  Ttv  Then  as  7^^  is  an  Abelian  group  under  0i,,  we  see  that 
Qv{-pP)  =  €  Tlv,  and  therefore  e_^(t,to)  exists. 

To  complete  the  proof,  then,  it  suffices  to  show  that  Cpit,  to)  solves  the  initial  value 
problem 

,v  - 

1  +  p''{t)i'it) 


y,  y{to)  =  1. 
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Let  y{t)  =  ep{t,to).  Then 

y{to)  —  ep(to)  ^o)  =  1- 

Now,  ep{t,to)  —  p{t)ep{t,to),  which  is  continuous.  Hence  by  Corollary  28,  ej (t,to)  = 
to),  and  we  get 

y^{t)  =  eJ(t,to) 

=  P^(Oep(i,io) 

=  pf>(t)  [ep{t,  to)  -  iy{t)e^ {t,  to)]  • 

Rearranging  this  equation  gives 

ej (t,  to)  [1  +  =  P^(*)ep(t,  to), 


2.2  Second-order  Linear  Dynamic  Equations 

Recall  that  we  are  interested  in  the  second-order  self-adjoint  dynamic  equation 

Lx  =  0  where  Lx{t)  =  [p(t)a:^(t)]^  -1-  q{t)x{t).  (2.1) 

Here  we  assume  that  p  ;  T  — >  R  is  continuous,  :  T  ^  R  is  Id-continuous  and  that 

p(t)  >  0  for  all  t  G  T. 

Define  the  set  D  to  be  the  set  of  all  functions  a:  ;  T  — >  R  such  that  ^  R  is 

continuous  and  such  that  [p(t)a;^]^  ;  TJJ  — >  R  is  Id-continuous.  A  function  a:  G  D  is 
said  to  be  a  solution  of  La:  =  0  on  T  provided  Lx{t)  =  0  for  all  t  e  T^. 

Now,  consider  the  second-order  linear  dynamic  equations 

Mix  =  0  where  Mix  =  x^'^  -f  pi{t)x^  +  P2{t)x,  (2.2) 

M2X  =  0  where  M2X  =  x^'^  -f  ai{t)x^  -t-  a2{t)x,  (2.3) 
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and 

M3X  =  0  where  M3X  =  x^'^  +  ri{t)x'^  +  r2{t)x^,  (2.4) 

where  Pi,ai,ri  :  T  — >  R  are  Id-continuous  for  i  G  {1,2}.  Take  Djw'  to  be  the  set 
of  all  functions  a:  :  T  ^  R  such  that  a;  is  A  -differentiable  on  T**,  a;^  :  T''  — >  R  is 
V  -differentiable  on  ,  and  a;^^  :  ^  R  is  Id-continuous.  For  i  =  1,2,3,  we  say 

a:  is  a  solution  of  MiX  =  0  on  T  provided  x  is  in  Dm,  and  MiX  =  0  for  all  t  G  T”. 

Theorem  32.  Ifp2  is  Id-continuous  andpi  G  then  the  dynamic  equation  (2.2) 
can  he  written  in  self-adjoint  form,  with 

p{t)  =  ep^{t,to)  and  q{t)  =  ep^{t,to)p2(t). 

Furthermore,  in  this  case,  if  x  is  a  solution  of  (2.2)  on  T,  then  x  is  also  a  solution 
of  the  self-adjoint  form  of  the  equation. 

Proof.  Suppose  we  have 

3:^^  -1-  Piit)x'^  +  P2{t)x  =  0. 

Assume  p2  is  Id-continuous  and  pi  Ellt-  Then  ep^(t,to)  is  well  defined  and  positive. 
Multiplying  through  by  epi(t,  to),  we  get 

epi(t,to)a;'^'^  +  ep^{t,to)pi{t)x'^  +  epi{t,to)p2{t)x  =  0. 

Then,  since  Cpj  (i,  to)  solves  the  IVP 

y^=Pi{t)y,  J/(«o)  =  l, 


we  have  that 

[epi  (t,  to)]^  =  Pi  (i)epi  {t,  to) . 

So  our  equation  becomes 


epi (i, to)x^'^  -1-  [epi(i, to)\^x'^  +  ep^{t, to)pi{t)x  =  0. 


Furthermore,  x^  is  V-differentiable,  hence  continuous,  so  by  Corollary  28,  a:^  =  x^p 
and  we  get 

{t,  to)x^^  +  [cpi  {t,  to)'^x^P  -f-  gpi  {t,  to)p2{t)x  =  0. 

Then  by  the  product  rule,  we  see  that 

[epi(t,  *0)2:^]^  +  Cpi  (t,  to)p2{t)x  =  0. 


This  equation  is  in  self-adjoint  form  with  p{t)  and  q{t)  as  desired. 

Now  suppose  a;  is  a  solution  of  (2.2),  p2  is  Id-continuous  and  pi  G  R'^.  Based  on 
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the  above  development,  it  is  clear  that  x  satisfies  the  dynamic  equation 

[epi{t,to)x^]'^  +  ep^{t,to)p2it)x  =  0. 

Hence  to  show  a:  is  a  solution  of  this  dynamic  equation,  we  need  only  show  that  x  G  D. 
Note  that  x^  is  V  -differentiable,  and  therefore  continuous.  Also, 

[epi(t,to)x^]^  =  e^^{t,to)x>’ +  ep^{t,to)x^'^ 

=  pi(t)4i  (*)  +  epi  (t,  to)x^'^, 


which  is  Id-continuous,  and  therefore,  x  G  D.  □ 

Corollary  33.  7/o2  is  Id-continuous  and  -oi  G  1Z^,  then  the  dynamic  equation  (2.3) 
can  be  written  in  self-adjoint  form,  with 


and  q{t)  = 


02  (t) 


1  -I-  ai{t)u{t) 


e  <»i 


l+aiu 


{t,to). 


Furthermore,  if  x  is  a  solution  of  (2.3),  then  x  is  also  a  solution  of  the  self-adjoint 
form  of  the  equation. 

Proof  Suppose  we  have 


x^^  -f  oi(t)x^  -t-  02(t)x  =  0. 

Recall  that  if  /  :  T  — »  M  is  V-differentiable,  then  f{t)  =  f^{t)  -\-  u{t)f^{t).  Thus 
u{t)x^^ .  Making  this  substitution,  we  have 

x^'^  ai{t){x^^  -I-  i/{t)x^'^)  -h  a2{t)x  =  0, 


and  hence 

(1  -f  ai{t)u{t))x^'^  -\-  ai{t)x^^  -\-  02(t)x  =  0. 


Now  — Oi(t)  G  'R.'fi  so  the  leading  coefficient  is  positive  and  we  may  divide  through 
by  it.  Furthermore,  as  before,  we  have  that  x^  is  continuous,  so  x^  =  x^^.  Thus,  we 
get 

fli(i)  V  I  ®2(i) 


AV  I 
X  + 


tX''  + 


{1 -\- ai{t)u{t))  (H-ai(t)i/(t)) 


X  =  0. 


This  is  in  the  form  (2.2).  As  oi  and  02  are  Id-continuous,  so  are  and 

Further, 


1 -\- ai{t)u{t)  J 


1 -\- ai{t)u{t)  -  ai{t)u{t)  ^  1  Q 

1  -h  ai{t)u{t)  1  -H  ai(t)z^(t)  ’ 


so  the  coefficient  of  the  x^  term  is  in  TZ'f.  Hence  by  Theorem  32  above,  the  equation 
can  be  written  in  self-adjoint  form,  with  p{t)  and  q{t)  in  the  desired  form,  and  solutions 
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of  equation  (2.3)  are  also  solutions  of  the  self-adjoint  form  of  the  equation.  □ 

Corollary  34.  Ifr2  is  Id-continuous  and  {ri  —i'r2)  €  then  the  dynamic  equation 
(2.4)  can  be  written  in  self-adjoint  form,  with 

p{t)  =  e(ri 

—ur2  ){t,to)  and  q{t)  =  r2{t)e(r^  —vr2  ){t,to). 

Furthermore,  if  x  is  a  solution  of  (2.4),  then  x  is  also  a  solution  of  the  self-adjoint 
form  of  the  equation. 

Proof.  Suppose  we  have 

x^'^  -h  ri(t)x^  -I-  r2{t)x^  =  0. 


Then 

-f-  ri(t)x^  -f  r2{t){x  —  u{t)x^)  =  0, 
or 

x^^  +  (ri(t)  -  u{t)r2{t))x^  +  r2(t)x  =  0. 

This  is  in  the  form  (2.2),  and  the  coefficients  meet  the  requirements  of  Theorem  32. 
Thus  the  result  follows.  D 

2.3  Abel’s  Formula  and  Reduction  of  Order 

We  begin  this  section  by  looking  at  the  Lagrange  Identity  for  the  dynamic  equation 
(2.1).  We  establish  several  corollaries  and  related  results,  including  Abel’s  Formula 
and  its  converse.  We  conclude  the  section  with  a  reduction  of  order  theorem.  Some 
of  the  results  in  this  section  are  due  to  Atici  and  Guseinov.  Specifically,  Theorems  35 
and  43,  and  Corollaries  39  and  42  were  previously  established  in  their  work  [4] .  Our 
conditions  on  p  and  q  are  less  restrictive  than  Atici  and  Guseinov’s,  and  our  domain 
of  interest,  D,  is  defined  more  broadly.  In  spite  of  this,  however,  many  of  the  proofs 
contained  in  [4]  remain  valid.  As  this  is  the  case,  we  have  omitted  the  proofs  of  some 
of  the  following  theorems,  and  refer  the  reader  to  Atici  and  Guseinov’s  work. 

Theorem  35.  If  to  €  T,  and  xq  and  xi  are  given  constants,  then  the  initial  value 
problem 

Lx  =  0,  x{to)  =  Xq,  x^{to)  =  Xi 
has  a  unique  solution,  and  this  solution  exists  on  all  o/T. 

Definition  36.  If  x,y  are  A-differentiable  on  T'*,  then  the  Wronskian  of  x  and  y, 
denoted  W{x,y){t)  is  defined  by 


W{x,y){t)  = 


x{t)  y{t) 
x^{t)  y^{t) 


for  t  e  T'^. 


22 


Definition  37.  If  x,  y  are  A-differentiable  on  T**,  then  the  Lagrange  bracket  of  x  and 
y  is  defined  by 

{x]y}{t)=pit)W{x,y){t)  fortGT". 

Theorem  38  (Lagrange  Identity).  If  x,y  e  then 

x{t)Ly{t)  -  y{t)Lx{t)  =  {x;  y}^{t)  for  t  E  T;^. 

Proof.  Let  x,  y  G  D.  We  have 

{x;y}^  =  Wix^vT 

=  \xpy^  -  ypx^^ 

=  x^jfy^P  +  x\py^]'^  -  y^jPx^P  -  y\px^^ 

=  x^jjPy^  +  x\^^^  -y'^fx'^  -y]jpx^^ 

=  x\py^]^  -y\px^]'^ 

=  x([pj/^]^  +  qy)  -  y{\px^^  +  qx) 

=  xLy  —  yLx, 

where  we  have  made  use  of  the  fact  that  x^  and  y^  are  continuous  and  applied 
Corollary  28.  □ 

Corollary  39  (Abel’s  Formula).  Ifx,y  are  solutions  of  (2.1),  then 

W(x,y){t)  =  ^^  /orier, 


where  C  is  a  constant. 

Proof  If  X,  y  are  solutions  of  (2.1),  they  belong  to  D.  Then,  by  Theorem  38,  we  have 
x{t)Ly{t)  -  y{t)Lx{t)  =  {x;  y}^(t)  for  t  G  T''. 

But  Lx  =  Ly  0,  so 

0  =  {x;j/}^(t)  for^GT^ 

Integrating,  we  see  that 


{x\y)  =p{t)W{x,y){t)  =  C, 
which  gives  the  desired  result. 

Definition  40.  Define  the  inner  product  of  x  and  y  on  [a,  b]  by 


(x,y)  :=  /  x{t)y{t)m. 
J  a 


□ 
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Corollary  41  (Green’s  Formula).  Ifx,ye  D,  then 

{x,Ly)  -  {Lx,y)  =  \p{t)W{x,y)]l. 

Proof.  Integrating  the  expression  in  Theorem  38  gives  the  result  immediately.  □ 
Corollary  42.  Ifx,y  are  solutions  of  (2.1),  then  either 
(i)  W{x^y)  ^  0  for  t  or 
(a)  W{x,y)  =  0  forte 

Case  (i)  occurs  if  and  only  if  x  and  y  are  linearly  independent  on  T,  and  case  (ii) 
occurs  if  and  only  if  x  and  y  are  linearly  dependent  on  T. 

In  the  standard  way,  one  uses  the  uniqueness  theorem  to  prove  the  following  result. 

Theorem  43.  If  xi  and  X2  are  linearly  independent  solutions  of  (2.1)  on  T,  then  a 
general  solution  of  (2.1)  is  given  by 


x{t)  =  CiXi{t)  +  C2X2{t). 


Theorem  44  (Converse  of  Abel’s  Formula).  Assume  u  is  a  solution  of  (2.1) 
with  u{t)  7^  0  for  f  €  T.  If  v  eB>  satisfies 


W{u,v){t) 


C 

Pity 


then  V  is  also  a  solution  of  (2.1). 

Proof  Suppose  that  u  is  a  solution  of  (2.1)  with  u{t)  7^  0  for  any  t,  and  assume  that 
V  €  D  satisfies  W{u,v){t)  =  Then  by  Theorem  38,  we  have 

u{t)Lv{t)  —  v{t)Lu{t)  =  {u;t;}^(t). 


so 


u{t)Lv{t)  =  \pitWi'>^’'^)it)r 

= 

=  0. 


As  u{t)  7^  0  for  any  t,  we  can  divide  through  by  it  to  get 

Lv{t)  =  0  for  t  G 


□ 


Hence  v  is  a  solution  of  (2.1)  on  T. 
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Theorem  45  (Reduction  of  Order).  Let  to  e  T'',  and  assume  u  is  a  solution  of 
(2.1)  unth  u{t)  ^  0  for  any  t.  Then  a  second,  linearly  independent  solution,  v,  of 
(2.1)  is  given  by  ^ 

v{t)  =  u{t)  f  .  . 

K  P{s)u{s)u‘^{s) 

for  t  e  T. 

Proof  By  Theorem  44,  we  need  only  show  that  6  ID  and  that  W{u,  v)(t)  =  ^  for 
some  constant  C.  Consider  first 

W(u,v)(t)  = 


=  u{t)v^{t)  —  v{t)u^{f) 

u{t)  U  + 

U  p{s)u{s)u<^{s)^^ 

-  u{t)u  p(5)„(s)„<T(s)^^  +  p(i) 

«(<)“  p{s)u{s)u'^{s)^^ 


p{f)u{t)u^{t) 


u(t)u‘'{t) 

t)u{t)u"(t) 


Here  we  have  (7=1.  It  remains  to  show  that  v  e  D.  We  have  that 


v^(t)  =  u^{t)  [  .  ,  ,  'A5+  ■  ■  jf} „i.\ 

Jto  P(s)u{s)u'^{s)  p{t)u{t)u^{t) 

r*  1  1 

^  Jto  p{s)u{s)u‘^{s)  p{t)u{t)  ■ 


Since  u  €  D,  u(t)  ^  0  and  p  is  continuous,  we  have  that  is  continuous.  Next, 
consider 


r  1  1 

-\-pP{t)u^^{t)  I  .  ,  — 

^  ^  ^  [Jto  p{s)u{s)u<^(s)  J 


^  n^{t) 


u{t)uP{t) 


Now,  the  first  and  last  terms  are  Id-continuous.  It  is  not  as  clear  that  the  center  term 
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is  Id-continuous.  Specifically,  we  are  concerned  about  whether  or  not  the  expression 


Uo  P{s)u{s)u<^{s)  J 

is  Id-continuous.  Note  that  the  integrand  is  rd-continuous.  Hence  Theorem  27  applies 
and  yields 


\f‘ _ L_ 

Jto  P{s)u{s)u‘^{ 


]  lims_»t_ 


te  A. 


First,  suppose  t  G  T^i.  Then  (T{p{t))  =  t,  so  our  expression  simplifies  to 


r _ 1 

Jto  p{s)u{s 


1  ^  1  1  1 
- Ac  =  -  =  . 

s)u'^{s)  J  pP{t)uP{t)u’^P{t)  pP{t)uP{t)u{t) 


Next,  suppose  that  t  e  A.  Then  t  is  left-dense,  and  therefore, 

lim  a{s)  =  t. 

s—*t— 

Then  as  p  and  u  are  continuous,  we  have 


s^t- p(s)u{s)u^{s)  p{t)v?{ty 

Now  for  t  e  A,  t  is  left-dense,  so,  if  we  like,  we  may  write  this  expression  as 

1  ^  1  ^  1 

s-»t- p(s)u(s)«‘^(s)  p(t)u^{t)  pP{t)uP(t)u{t) 

This  is  the  same  expression  we  got  for  t  G  T^,  so  we  have  that 


/  ,  '  .  . — = — 7^ — ,  '  for  t  G  T. 

[Jto  p{t)u{t)u^{t)  j  pP{t)uP{t)u{t) 

This  function  is  Id-continuous,  and  so  we  have  that  u  G  D.  Hence  by  Theorem  44,  v 
is  also  a  solution  of  (2.1).  Finally,  note  that  as  W{u,v){t)  =  ^  ^  0  for  any  t,  u  and 
V  are  linearly  independent.  □ 


2.4  Oscillation  and  Disconjugacy 

In  this  section,  we  establish  results  concerning  generalized  zeros  of  solutions  of  (2.1), 
and  examine  disconjugacy  and  oscillation  of  solutions. 


Definition  46.  We  say  that  a  solution,  x,  of  (2.1)  has  a  generalized  zero  at  t  if 

x{t)  =  0 
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or,  if  <  is  left-scattered  and 

x{p{t))x{t)  <  0. 

Definition  47.  We  say  that  (2.1)  is  disconjugate  on  an  interval  [a,  6]  if  the  following 
hold. 


(i)  If  re  is  a  nontrivial  solution  of  (2.1)  with  x{a)  =  0,  then  x  has  no  generalized 
zeros  in  (a,  6]. 

(ii)  If  rc  is  a  nontrivial  solution  of  (2.1)  with  :E(a)  ^  0,  then  x  has  at  most  one 
generalized  zero  in  (a,  6]. 

Definition  48.  Let  uj  =  sup  T,  and  if  a;  <  oo,  assume  p{oj)  =  u.  Let  o  6  T.  We  say 
that  (2.1)  is  oscillatory  on  [a,a;)  if  every  nontrivial  real- valued  solution  has  infinitely 
many  generalized  zeros  in  [a.w).  We  say  (2.1)  is  nonoscillatory  on  [a,t<;)  if  it  is  not 
oscillatory  on  [a,a;). 

Lemma  49.  Let  u  —  supT.  If  u)  <  oo,  then  assume  p{u})  =  u).  Let  a  €  T.  Then  if 
(2.1)  is  nonoscillatory  on  \a^uf),  there  is  some  to  G  T,  to  ^  such  that  (2.1)  has  a 
positive  solution  on  [toj^^)- 

Proof.  Assume  (2.1)  is  nonoscillatory  on  [a,  u),  and  then  there  is  a  nontrivial  solution, 
u  of  (2.1)  such  that  u  has  only  finitely  many  generalized  zeros  in  [o,a;).  Let  b  = 
max{t  e  T  ;  u  has  a  generalized  zero  at  t}.  Fix  to  G  T  such  that  to  >  b.  Then  either 
w  >  0  on  [to,a;)  or  — u  >  0  on[to,tc;).  □ 

Theorem  50  (Sturm  Separation  Theorem).  Letu  and  v  be  linearly  independent 
solutions  of  (2.1)  on  T  .  Then  u  and  v  have  no  common  zeros  in  T.  If  u  has  a 
zero  at  ti  G  T,  and  a  generalized  zero  at  t2  >  ti  G  T,  then  v  has  a  generalized  zero  in 
(ti ,  t2] .  If  u  has  generalized  zeros  at  ti  G  T  and  t2  >  ti  G  T,  then  v  has  a  generalized 
zero  in  [ti,t2]. 

Proof.  If  u  and  v  have  a  common  zero  at  to  G  T”,  then 


W{u,v){tQ) 


uito)  v{to) 
n^(to)  v^{to) 


Hence  u  and  v  are  linearly  dependent. 

Now  suppose  u  has  a  zero  at  ti  G  T,  and  a  generalized  zero  at  t2  >  ti  G  T.  Without 
loss  of  generality,  we  may  assume  t2  >  cr(ti)  is  the  first  generalized  zero  to  the  right 
of  ti,  u(t)  >  0  on  (ti,t2),  and  u(t2)  <  0.  Assume  u  is  a  linearly  independent  solution 
of  (2.1)  with  no  generalized  zero  in  (ti,t2].  Without  loss  of  generality,  v(t)  >  0  on 
[ti,t2]. 
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Then  on  [ti,t2], 


u\^  (.s  ^  v{t)u^{t)  -  u{t)v^{t)  ^  C 
v)  ^  v(t)v‘^(t)  p(t)v(t)v‘^(t)’ 


which  is  of  one  sign  on 
that 


[ii,t2)-  Thus  ^  is  monotone  on  [ti,t2].  Fix  G  (<1,^2)-  Note 


ujti) 

v{ti) 


=  0,  and 


^(^3) 

v(t3) 


>0. 


But 


U{t2) 

V{t2) 


<0, 


which  contradicts  the  fact  that  ^  is  monotone  on  [ti,i2]-  Hence  v  must  have  a  gener¬ 
alized  zero  in  (ti,t2]- 


Finally,  suppose  u  has  generalized  zeros  at  ti  G  T  and  <2  >  €  T.  Assume 

t2  >  is  the  first  generalized  zero  to  the  right  of  ii.  If  ti(ti)  =  0,  we  are  in 

the  previous  case,  so  assume  u{t\)  7^  0.  Then,  as  it  has  a  generalized  zero  at  ti,  we 
have  that  t\  is  left-scattered.  Without  loss  of  generality,  we  may  assume  u{t)  >  0 
on  [ti,t2),  u{p{ti))  <  0  and  u{t2)  <  0.  Assume  u  is  a  linearly  independent  solution 
of  (2.1)  with  no  generalized  zero  in  [<1,^2)-  Without  loss  of  generality,  v{t)  >  0  on 
[ti,t2],  and  v{p{ti))  >0.  In  a  similar  fashion  to  the  previous  case,  we  apply  Abel’s 
Formula  to  get  that  ^  is  monotone  on  [p(ti),t2]-  But 


ujpjti))  ^ 

vipih))  ^ 


ujh) 

v{ti) 


>0, 


and 


u{t2) 

vit2) 


<0, 


which  is  a  contradiction.  Hence  v  must  have  a  generalized  zero  in  [<1,^2]- 


□ 


Theorem  51.  If  (2.1)  has  a  positive  solution  on  an  interval  X  C  T,  then  (2.1)  is 
disconjugate  onX.  Conversely,  ifa,b  G  TJJ  and  (2.1)  is  disconjugate  on  [p{a),a{b)]  C 
T,  then  (2.1)  has  a  positive  solution  on  [p{a),a{b)]. 

Proof  Assume  (2.1)  has  a  positive  solution,  w  on  J  C  T.  If  (2.1)  is  not  disconjugate 
on  X,  then  (2.1)  has  a  nontrivial  solution  v  with  at  least  two  generalized  zeros  in  X. 
Then,  without  loss  of  generality,  there  are  ti,t2  in  X  such  that 


v{ti)  <  0,v{t2)  <  0,  and  v{t)  >  0  on  (ti,t2)  with  (^1,^2)  ^  0- 


=  u{t)v^{t)  -  v{t)u^{t) 
u)  '  u{t)u’^{t) 

W{u,v){t) 
u{t)u‘^{i) 


Note  that 
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C 

p{t)u{t)u‘^  {t) 

is  of  one  sign  on  T'^.  Hence  ^  is  monotone  on  J.  But 

©  -  “■ 

This  contradicts  the  fact  that  ^  is  monotone.  Hence(2.1)  is  disconjugate  on  J. 

Conversely,  suppose  that  (2.1)  is  disconjugate  on  the  compact  interval  [p(a),cr(6)]. 
Let  u,v  be  the  solutions  of  (2.1)  satisfying  u{p{a))  =  0,u^{p(a))  =  1  and  v(a-(b))  = 
0,v^(b)  =  —1.  Since  (2.1)  is  disconjugate  on  [p{a),a{b)],  we  have  that  u{t)  >  0  on 
(p(a),cr(6)],  and  v{t)  >  0  on  [p(a),(T(6)).  Then 

x{t)  =  u{t)  +  v{t) 


is  the  desired  positive  solution.  □ 

Theorem  52  (Polya  Factorization).  If  (2.1)  has  a  positive  solution,  u,  on  an 
interval  X  CT,  then  for  any  x  eB,  we  get  the  Polya  Factorization 

Lx  =  a:i(t){Q:2[Q!ix]^}^(i)  for  tel, 


where 


and 


Q!i  :=  —  >  0  on  X, 
u 


02  :=  jmu'^  >0  onX. 


Proof.  Assume  that  u  is  a  positive  solution  of  (2.1)  on  X,  and  let  a:  €  D.  Then  by 
the  Lagrange  Identity  (Theorem  38), 


u{t)Lx{t)  —  x{t)Lu{t) 
u{t)Lx(t) 

Lx(t) 


{u-,x}'^{t) 

{u;x}^(i) 


u{t) 

1 

u{t) 

1 

u{t) 


{u-,x}^it) 


{pWiu,x)}^it) 


(t) 


ai(t){o2[«ix]^}^(t). 


for  t  eX,  where  a\  and  02  are  as  described  in  the  theorem. 


□ 


29 


Theorem  53  (Trench  Factorization).  Let  a  €  T,  and  let  u)  :=  supT.  If  ui  <  oo, 
assume  pioS)  =  u).  If  (2.1)  is  nonoscillatory  on  [a,  w),  then  there  is  to  ^  T  such  that 
for  any  x  G  D,  we  get  the  Trench  Factorization 

fort  G  [to)^);  where  Pi,P2>  ^  on  [to,(jj),  and 

fW  ^ 

/  ^  .  X  At  =  oo. 

Jto  P^it) 

Proof.  Since  (2.1)  is  nonoscillatory  on  [a,a;),  (2.1)  has  a  positive  solution,  u  on  [to,  w) 
for  some  to  G  T.  Then  by  Theorem  52,  Lx  has  a  Polya  factorization  on  [to,w).  Thus 
there  are  functions  a\  and  02  such  that 

Lx{t)  —  Q:i(t){a2[Q!ia:]^}''^(t)  for  t  G  [to,cc>). 


with 


Now,  if 


Oi  =  —  and  02  =  pu‘fjP- 

u 


r  1 


>tQ  Oi2{t) 

then  take  /?i(t)  =  Q:i(t),  and  ^2(i)  =  «2(i),  and  we  are  done.  Therefore,  assume  that 


r — 

Jto  «2(0 


At  <  00. 


In  this  case,  let 


Piit)  =  and  /32(t)  =  0(2 (t)  [  — yrAs  /  — 

It  "2(5)  J,it)  «2(s) 

for  t  G  [to,w).  Note  that  as  0:1,02  >  0,  we  have  ,0i,/32  >  0  as  well.  Also, 


r  1  A*  i--  ^ 

L  W)‘'*  "  ‘JSeri.  a,(t)jr^A5/; 

,.  0:2  (s 


0-2(8)  Jo-(t)  02(8) 

1 

“2  (8_) _  A  . 

^As 


As 


At 


=  lim 


t 

tn  / 


oc2(s)  0(2 (s) 

A 

1 


to  il  012  (S) 


■UJ  1 


As 


At 
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=  lim 

&— >tJ,6€T 

=  oo. 


ft 


CJ  1 


b  a2(s) 


As 


Now  let  a:  G  D.  Then 


ai{t)x{t) 


'^^As 


n  A 


ft  0:2  (s) 
for  t  e  So  we  get 


(t) 


ft  aM^^fc 


■w  1 


t  a2(s)  J<T(t)  02(5) 


As 


p2{tWi{t)A^  =  a2(i)K  ’’’ 

for  t  e  [to,u)).  Taking  the  V  -derivative  of  both  sides  gives 

+  {a2(t)[a.(f)x(«)]^}'  [jT 
-|-[Q;i(i)a:(t)]^ 

for  t  e  We  now  claim  that  the  last  two  terms  in  this  expression  cancel. 

Looking  only  at  these  last  two  terms,  put  the  expression  back  in  terms  of  our  positive 
solution  u.  We  get 


\p{t)u{t)u^{t)Y 


x{ty 

L“W. 


Ap 


■  r _ 1 

Jt  Pis)u{s 


nV 


s)u'^(s) 


As 


+ 


x(t) 


Now  consider  two  cases: 

Case  1:  t  G  T^.  Then  Theorem  26  applies,  and  we  get 


pP{t)u>’{t)u{t) 


'x{ty 

Ap  - 

A't). 

r  1 
Jt  p{s)u{s)u^ 


tV 


(5) 


As 


+ 


x{t) 


nV 


w(t) 


pP{t)uP{t)u{t) 


+ 


x{t) 


tV 


u{t) 


x{t) 


[u(t) 


+ 


u{t) 


=  0. 


31 


Case  2:  t  E  A.  In  this  case  we  have  that  p{t)  =  t,  and  we  get 


\p{t)u{t)u'^{t)Y 


_ _ As 

[Jt  p{s)u{s)u‘^{s)  _ 


+ 


\p{t)u{t)u'^{t)]  u{t)x^{t)  -  x{t)vy{t) 

p{t)u^{t)  u^{t) 

u{t)x^^{t)  -  x{t)u^^{t)  u{t)x'^{t)  -  x(t)v7{t) 

u^{t) 

—u{t)x^{t)  +  x{t)vy{t)  +  u{t)x^{t)  —  x{t)u^{t) 


=  0. 


Here,  we  have  made  use  of  the  fact  that  a;,  u  G  D,  which  gives  us  that  x^p  =  x^  and 
u^p =  . 

In  either  case,  the  last  two  terms  cancel,  and  we  have  that 

It  then  follows  that 

Piit)  {/02(t)[A(i)ic(t)]^}^  =  ai(t)  {a2(t)[aiWaj(i)]^}'^  =  Lx{t), 

for  t  G  [to,u})  and  the  proof  is  complete.  □ 

Theorem  54  (Recessive  and  Dominant  Solutions).  Let  a  G  T,  and  let  u  := 
supT.  If  oj  <  oo,  then  we  assume  p{oj)  =  u.  If  (2.1)  is  nonoscillatory  on  [a,u)), 
then  there  is  a  solution,  u,  called  a  recessive  solution  at  oj,  such  that  u  is  positive  on 
[to,w)  for  some  to  G  T,  and  if  v  is  any  second,  linearly  independent  solution,  called  a 
dominant  solution  at  oj,  the  following  hold. 

(i)  limt^o;-  ^  =  0 

/to  p{t)u(t)u<’(t)^^  ~  °° 

Ib  p(t)v(t)v'^(t)^^  <6oforb<0J,  sufficiently  close,  and 
(iv)  fof  t  <oj,  sufficiently  close. 

The  recessive  solution,  u,  is  unique,  up  to  multiplication  by  a  nonzero  constant. 
Proof.  As  (2.1)  is  nonoscillatory,  by  Theorem  53,  there  is  a  Trench  Factorization: 

Lx{t)=m{Wixf}'’{t), 
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where  ^1,^02  >  0  on  [to,u)),  and 

r  1 

/  -"-T-rAt  =  OO. 

Jto 

Then  if  u{t)  =  tt  is  a  positive  solution  of  (2.1).  Now,  let 


Then, 

Lvo{t)  =  l3iit){l32[f3iVoff{t) 

=  ft(t){i}’ 

=  0. 


So  uo  is  a  solution  of  (2.1).  Note  that 


v  «(*) 

hm  — ^ 

t-w-  Vo{t) 


lim  f., 

to-a; 


-^A<i 
to  02{s)^^ 


\uJ  u(t)u‘^(t)  p(t)u(t)u‘^(t) 

where  C  is  a  constant  by  Theorem  39.  Note  that  C  ^  0,  since  u  and  vq  are  linearly 
independent.  Integrating  both  sides  of  this  last  equation  from  <o  fo  t,  we  get 

_ A. 

u{t)  Jt^  pis)u{s)u‘^{s) 

Taking  the  limit  as  t  o’,  we  get 


lim 

t—mj 


vo{t) 

u{t) 


C 

p{s)u{s)u^{s) 


As, 
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and  we  see  that 


f 


C 


-As  =  oo, 


/to  p(s)u(s)u'^(s) 

as  desired. 

Now  let  V  be  any  solution  of  (2.1)  such  that  u  and  v  are  linearly  independent. 
Then 

v{t)  =  Ciu{t)  +  C2Vo{t),  where  C2  7^  0, 

and 


lim 


u{t) 


lim 


u{t) 


t-^u-  v{t)  t-+w-  Ciu{t)  +  C2Uo(^) 


=  lim 


fo(t) 


-  ^  , 


=  0. 


Now,  let  u  be  a  fixed  solution  of  (2.1)  such  that  u  and  v  are  linearly  independent. 
Choose  ti  e  [to)^^)  such  that  v{t)v°'{t)  >  0  on  [ti,a;).  Then  for  t  6  [fi,a;), 


W{v,u){t) _ Cl 


p{t)v{t)v'^  {ty 


Cl 


where  Ci  7^  0.  Integrating, 

u(t)  u{ti) 
v[t)  v{ti)  pis)v{s)v‘^{s) 

Letting  t  — t  w— ,  we  see  that 

^(fi)  r  c-x 


=/ 


As. 


which  implies  that 


Furthermore,  for  t  E.[ti,uj), 


ti)_  r 

vih)  Jt^  p{s)v{s)v<^{s) 

r  1 

Jti  p{s)v{s)v<^{s) 


As, 


As  <  00. 


p{t)vy  _  p(t)u^(^  ^  pjtmuvm  ^  ,,here  C,  5^  0. 

v(t)  u(t)  u(t)v(t)  u(t)v(t) 

It  remains  to  show  that  C2  >  0.  We  have 

liixx  -77:  =  00, 

t-n^  u{t) 
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/^nA  Wiu,v){t)  C2 

\u/  ^ u{t)u'^{t)  p{t)u{t)u'^{t)’ 

which  implies  that  C2  >  0,  as  desired. 

Finally,  we  need  to  establish  uniqueness,  up  to  multiplication  by  a  nonzero  con¬ 
stant.  Let  ui  be  a  recessive  solution  of  (2.1),  and  suppose  U2  is  another  recessive 
solution.  If  ui  and  U2  were  linearly  independent,  U2  would  be  a  dominant  solution. 
Hence  ui  and  U2  must  be  linearly  dependent,  and  we  see  that  U2  =  ku2  for  some 
nonzero  constant  k.  □ 


2.5  The  Riccati  Equation 


Usually,  linear  dynamic  equations  are  considerably  easier  to  solve  than  nonlinear 
ones.  In  this  section,  we  are  going  to  discuss  the  relationship  between  a  particular 
nonlinear  equation,  called  the  Riccati  equation,  and  our  self-adjoint  equation.  We 
will  see  that  there  is  a  correspondence  between  solutions  of  these  two  equations.  The 
Riccati  equation  is  defined  by 

(zP(t))^ 

Rz  =  0,  where  Rz(t)  ■-  z^{t)  +  q(t)  +  (2-5) 


for  t  €  T”.  Here  we  assume  that  p  :  T  M  is  continuous,  ?  :  T  — >  R  is  Id-continuous 
and  that 

p{t)  >  0  for  all  t  G  T. 

Define  the  set  to  be  the  set  of  all  functions  z  :  T''  — >  R  such  that  :  TJ^  — »•  R  is 
Id-continuous  and  such  that  p^{t)  +  i/(t)z^(t)  >  0  for  any  t  G  TJJ.  A  function  z  G  Dk 
is  said  to  be  a  solution  of  Rz  =  0  on  T'^  provided  Rz{t)  =  0  for  all  t  G  TJ^. 

Then  we  have  the  following  theorem: 


Theorem  55.  Assume  a;  G  D  has  no  generalized  zeros  in  T,  and  z  is  defined  by  the 
Riccati  substitution 

<t)  =  m 


x{t) 


for  t  G  T”.  Then  z  G  D/j,  and 


Lx{t)  =  x{t)Rz{t) 


for  t  G 

Proof  We  first  wish  to  show  that  z  G  D^.  We  have  by  the  quotient  rule 


\t) 


p{t)x^{t) 


x{t) 


V 


x{t)\p{t)x^{t)f  -  pit)x^{t)x^ jt) 
x{t)xp{t)  ' 
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which  is  Id-continuous  on  TJ^,  since  x  €  1D>.  Next,  note  that 

pp{t){xP{t)  +  v{t)x^{t)) 

Xp{t) 

XP{t) 

for  all  t  €  T”,  since  x  has  no  generalized  zeros  in  T.  It  remains  to  show  that 
x{t)Rz{t)  =  Lx{t)  for  t  e  Suppressing  the  arguments,  we  get 


xRz  =  X 


=  X 


=  X 


,p\2 


pP  -I-  UZP 


+  q  + 


pP  +  iy(s^y 


px 


,A\  P\  2 


X 


x{px^)^  —  px^x'^ 


XXP 


+  q  + 


Xp 

pp{xp  +  ux^p)  (xpy 
xpp{x'^y 


.2-] 


,  AxV  X  px^x^ 

=  (px^)^ - \-qx+  , 

^  xP  xP  xP{xP  +  ux^) 


px^x^ 

XP 


+  qx  + 


XpP{x^) 

XPx 


{px^) 


px^x^  ^  pP{x'^)^ 


=  (px^)'^  +  qx  + 

\y  j  -T'i  -T 

{px^y^vx'^  —px^x^  +pP{x'^){x^P) 


=  Lx  d" 

=  Lx  d- 

=  Lx  + 
=  Lx. 


xP 


x^{pPx^P  +  u{px^)'^)  —  px^x'^ 


xP 

,A^V 


x^{px^)  —px^x 
xP 


Again,  we  have  made  use  of  the  fact  that  a:  G  D  and  applied  Corollary  28.  □ 

Theorem  56.  The  self-adjoint  equation  (2.1)  has  a  positive  solution  onT  if  and  only 
if  the  Riccati  equation  (2.5)  has  a  solution  z  on  T**. 

Proof  First,  assume  that  a:  is  a  positive  solution  of  (2.1),  and  let  z  be  defined  by  the 
Riccati  substitution  (2.6).  Then  by  Theorem  55,  z  G  Djj,  and  Lx  =  xRz.  Since  x  is 
a  solution  of  Lx  =  0  and  has  no  generalized  zeros,  it  follows  that  Rz  =  0,  as  desired. 

Conversely,  assume  that  z  is  a  solution  of  the  Riccati  equation,  (2.5)  on  T'^.  Then 
z  G  Dfi,  so  pP{t)  d-  u{t)zP{t)  >  0  for  all  t  G  T^,  and  2:  is  continuous  on  T''.  This  gives 
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us  that  —  and  thus,  by  Lemma  30,  |  €  TV^.  Now,  let  to  G  T,  and  let  x 

be  the  solution  of  the  initial  value  problem 

=  ^x,  x{to)  =  1. 

PW 

Note  that  although  ^  is  only  defined  on  T**,  x  is  defined  on  T.  Furthermore,  as 
x{t)  =  es.{t,to),  X  is  continuous  and  positive  on  T.  Next,  consider 

[p{t)x^{t)]^  =  [z{'t)x{t)f 

=  z'^{t)x^{t)  +  z{t)x'^{t) 

=  z^{t)x^{t)  +  z{t)x^^{t), 


which  is  Id-continuous  on  T^.  Hence  a;  G  D.  Moreover,  we  see  that 

so  by  Theorem  55  Lx  =  xRz  =  0.  Hence  x  is  the  desired  positive  solution  of  (2.1).  □ 

Now,  define  A  to  be  the  set  of  functions 


A  :=  {u  G  Clid{[p{a),a{b)],R)  :  u{p{a))  =  u{a{b))  =  0}. 

Here,  Cpi^  denotes  the  set  of  all  continuous  functions  whose  V-derivatives  are  piece- 
wise  Id-continuous.  Then  we  define  the  quadratic  functional  T  on  A,  by 


ra{b) 

Jp{a) 

Definition  57.  We  say  T  is  positive  definite  on  A  provided  ^(u)  >  0  for  all  «  G  A, 
and  JF(it)  =  0  if  and  only  if  u  =  0. 

Lemma  58  (Completing  the  Square).  Assume  z  is  a  solution  of  the  Riccati 
equation  (2.5)  on  [p(u),  b].  Let  u  E  A.  Then  for  all  t  G  [a,  b],  we  have 


{zu^Y{t)  =  pP{t){u'^{t)f -q{t)u^{t) 

z^{t)u(f) 

y/pf>{t)  +  u\t)zP{t) 


n2 


^pP{t)  +  u(t)zP{t)u'^  (t) 


Proof.  Let  z  be  a  solution  of  the  Riccati  equation  (2.5)  on  [p(o),6],  and  let  u  G  A. 
Then  for  t  E  {a,  b], 


z^{t){u\t))+zP{t){u\t)r 
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=  z{t)'^  {u^{t))  +  z^{t){u^{t)v7{t)  +  u{t)u^(t)) 

+z^{t)u^{t)v7  (t)  +  z^{t)u{t)v7{t) 

2.  ^  (zP(t)fu\t) 

-  W  pp(i)  +  v{t)zP(t) 

+z^{t)u{t)v7  (t)  +  z'^{t)v7  {t){u{t)  —  I/(t)u^(t)) 

=  ,(i]At] 

+2zP{t)u{t)v7  (t)  —  z^{t)u{t){v7  {t)y 

/./  \/  V/  nn2  /  n  2/  \ 

=  ?'(*)(«  W)  -«(*)“  W 

+2z^{t)u{t)u^ (t)  -  {pP{t)  +  z^{t)u{t)){vy {t))^ 


-q{t)u^{t)  - 


f{t){v:^{t)f-q{t)n\t)- 


^pp(t)  +  v{t)zP{t)u^  {t) 


Theorem  59.  Let  x  be  a  solution  of  (2.1)  on  [p(a),  (t(6)],  and  let  c,  d  E  [p(a),o'(6)] 
with  p{a)  <c<  a{c)  <d<  <j{h).  If  c  =  p{a),  assume  x{c)  =  0.  Now,  let 

{0  p(a)  <t  <  c 
x{t)  c<t  <d 

0  d  <t  <  cr(b). 

Then  u  E  A,  and  J^u  =  C  +  D,  where 

{  -p{c)x^{c)x{c)  v{c)  =  0 

C7  -  j  p'-(c)x(yp(c)  ^ 

J  p[d)x^{d)x{d)  v{d)  =  0 

^  -  I  uld)  >  0. 

Proof.  Let  x,  «  be  as  described  in  the  statement  of  the  theorem.  We  first  claim 
that  u  E  A.  It  is  apparent  from  the  definition  that  u  E  Cpj^([p(a),<T(6)],M),  and 
that  u{a{b))  =  0.  The  fact  that  u{p{a))  =  0  is  also  clear  from  the  definition  unless 
p(a)  =  c.  In  this  case,  however,  u{p{a))  =  u{c)  =  x{c)  =  0,  by  our  assumption  on  x. 
So,  «  €  A,  as  desired.  Now  consider 

^<7(6) 

[pP{t){u'^{t)f  -  q{t)u'^{t)]  Vi. 

Jp(a) 
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We  have  u{t)  =  0  on  [p(o),  c)  U  [d,  a{b)],  and  =  0  on  [p{a),c)  U  {d,  (r(6)],  so  we  get 
J^u=  f  [p^{t){u^(t)Y  -  q{t)u^{t)]  Vt. 

Jpic) 

Breaking  up  the  integral,  we  get 


Now,  we  apply  Lemma  29  to  get 

Tu  =  f{c){vy{c))M(^)+ifid){v7{d))Md) 

—q{c)u^{c)i'{c)  —  q{d)xj?{d)v{d) 

rp{d)  pp{d) 

+  /  f(t)(u'=(t)fvt- q{t)u\t)Vt. 

Since  u{d)  =  0,  the  fourth  term  in  this  expression  vanishes.  Furthermore,  u{t)  =  x{t) 
on  [c,d),  and  on  (c,  d),  thus  we  may  substitute  x  for  u  in  the  two 

remaining  integrals.  We  make  this  substitution  and  then  evaluate  the  first  of  the  two 
remaining  integrals  by  parts,  which  yields 

Tu  =  p'’(c)(u^(c))^i/(c) +j/(d)(u^(d))^^'(d) 

—q{c)t^{c)v{c)  —  q{d)v?' {d)v{d) 

/p{d)  pp{d) 

J  q(t)x\t)Vt 

=  /(c)(«’(c))V(c) 

-q(c)v‘{c)lf(c)  +  f(l>{S))x‘‘{p(lI))x(p{lI}) 

rp{d)  — 

—p{c)x^(c)x{c)  —  /  [p(t)x^(t)]  x{t)  Vt 

fpid) 

+  J  vt 

=  p^(c)(u^(c))V(c)  +f{d){v7{d)fu{d) 

-q{c)u^{c)u{c)  +  p{p{d))x^{p{d))x{p{d)) 

rp(d) 

—p{c)x^{c)x{c)  —  j  x{t)Lx{t)  Vt 
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=  f{c){u^{c))Mc)  +f{d){v:^{d))Md) 

-q{c)u^{c)u{c)  +  p{p{d))x^{p{d))x{p{d)) 
—p{c)x^{c)x{c) 

=  C  +  D, 

where 

C  -  iy{c)p'’{c){u^(c)f -u{c)q{c){u{c)f 
-p{c)x^{c)x{c), 


and 

D  =  i^(dK(d)(«'^(d))^  +p(p(d))a:^(p(d))a:(p(d)). 

Note  that  if  u{c)  =  0,  then  C  =  -p{c)x^{c)x{c).  If  v{c)  >  0,  then  c  is  left-scattered, 
so  we  get 


C  =  u{c)jf{c) 


u{c)  -  u{p{c)) 


T  2 


-  ^{c)q{c)u^{c) 


—p{c)x^{c)x{c) 


=  ^  -  v{c)q{c)x^{c)  -  p{c)x^{c)x{c) 

l^\C) 

-t-]/(c)a;^^(c)a:(c)  —  p^(c)x^'’(c)x(c) 

=  ^  “  p'(c)x‘''’{c)x{c) 


—u{c)x{c) 
pf’{c)x^{c) 


I/(c) 


—  p^{c)x^  {c)x{c) 


I/(c) 

-u{c)x{c)  {q(c)x{c)  +  [pa;^]^(c)) 

p^{c)x^{c)  —  pP{c)x^{c)  +  jf{c)x{c)xf*{c) 
u{c) 

_  p'’{c)x(c)x'’{c) 
u{c) 


so  C  is  as  described  in  the  statement  of  the  theorem.  Now  note  that  if  v{d)  =  0,  then 
D  =  p{p{d))x^{p{d))x{p{d))  =  p{d)x^{d)x{d).  If  i/(d)  >  0,  then  d  is  left-scattered, 
so  we  get 


D 


u{d)pp{d){vy{d)f + ifid)x'^{d)x»{d) 
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=  i^{d)pP{d) 


u{d) 


v{d) 


pP{d){xP{dyf  pP(d)x^{d)x(d)  pP{d){x^(d))^ 

l^d)  lid)  Kd) 

p^{d)x^{d)x{d) 

v{d) 


Thus  jD  is  as  desired,  and  the  proof  is  complete. 


□ 


Theorem  60  (Jacobi’s  Condition).  The  self-adjoint  equation  (2.1)  is  disconjugate 
on  [p{a),o{h)]  if  and  only  if  T  is  positive  definite  on  A. 


Proof.  First,  suppose  (2.1)  is  disconjugate  on  [p(a),cr(6)].  Then  there  is  a  positive 
solution,  X,  of  (2.1)  on  [p{a),(r{b)].  Let  z{t)  :=  Then  by  Theorem  56,  z  is  a 

solution  of  =  0  on  [/9(a),  6].  Thus  by  Lemma  58,  for  any  u  E  A, 


{zit)u\t))^ 


j/{t){u'^{t))^  -  q{t)u^{t) 
Zp{t)u{t) 

Vp^{t)  +  ^i^)z'’i^) 


2 


^/pP{t)  +  u{t)zP{t)u^  (t) 


for  t  E  [a,  6].  In  fact,  it  can  be  shown  that  this  equation  holds  at  t  =  a{b)  as  well.  As 
the  equation  holds  on  [a,  cr(6)],  we  may  integrate  from  p{a)  to  cr{b),  and  noting  that 
u{p{a))  —  u{a{b))  =  0,  we  get 


zP{t)u{t) 


n  2 


^/pp^)+u(J)zP(l) 


y/pP{t)  +  u{t)zP{t)u^  (t) 


Vt, 


so  Tu  >  0  for  all  u  E  A.  Furthermore,  it  is  clear  that  if  «  =  0,  then  Tu  =  0.  Now 
suppose  !Fu  =  0.  Then 


=  y/pP{t)  +  u(t)zP{t)u"^{t), 

■s/pP{t)  +  v{t)zP{t) 


so  u  solves  the  initial  value  problem 


u 


V  _ 


pP  +  VZP 


u,  u{a{b))  =  0 


on  [a,  a{b)]  Since  ^  solution  of  this  IVP  is  unique,  and  gives  u{t)  =  0 

on  [a, <7(6)].  As  u{p{a))  =  0  as  well,  we  get  u{t)  =  0  on  [p{a),a{b)].  Hence,  is 
positive  definite  on  A. 

We  will  prove  the  converse  of  this  statement  by  contrapositive.  Suppose  (2.1) 
is  not  disconjugate  on  [p{a),a{b)].  Then  there  is  a  nontrivial  solution  x  of  (2.1) 
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such  that  either  x(p{a))  =  0  and  x  has  a  generalized  zero  in  (p(o),cr(6)],  or  x  has 
two  generalized  zeros  in  {p{a),a{b)].  In  either  case,  let  c  <  d  be  the  two  smallest 
generalized  zeros  of  x  in  [p{a),a{b)].  Then,  let 

{0  p{a)  <t<c 
x{t)  c  <t  <  d 
0  d  <t<  cr{b). 

Applying  Theorem  59,  we  then  have  ^  =  C  +  £>  <  0.  As  w  is  not  identically  0,  this 
tells  us  that  !F  is  not  positive  definite.  By  contrapositive,  the  proof  is  complete.  □ 

Theorem  61  (Sturm  Comparison  Theorem).  Let 

L\x  -  \pi{t)x^]'^  +  qi{t)x, 

L^x  =  [p2(t)a;^]^  +  q2{t)x- 

Assume  qi{t)  >  q2{t)  and  0  <  pi{t)  <  P2{t)  for  t  €  [p{a),a{b)].  If  Lix{t)  =  0  is 
disconjugate  on  [p{a),a(b)],  then  L2x{t)  =  0  is  disconjugate  on  [p(a),o-(6)]. 

Proof.  Let 

p<T{b) 

/  [rf(t)(«’W)"-«l(«y(f)]V(, 

Jp(a) 

pa{b) 

^2{u)  :=  /  -  q2{t)u^{t)]  Vt. 

Jp(a) 

Assume  that  Lix{t)  =  0  is  disconjugate  on  [p{a),a{b)].  Then  by  Theorem  60,  the 
quadratic  functional  is  positive  definite  on  A.  Then,  for  u  €  A,  we  have 

T2U  =  /  [p^{t){u'^{t)f -q2{t)u^it)]Vt 

Jp{a) 
path) 

Jp(a) 

=  J^iU. 

Hence  T2  is  positive  definite  on  A.  Then,  again  by  Theorem  60,  L2X  =  0  is  disconju¬ 
gate  on  [p(a),cr(6)].  □ 

Theorem  62.  Let  Lix,  L2X  be  as  in  the  Sturm  Comparison  Theorem.  Then  if  LiX  = 
0  is  disconjugate  on  [p(a),cr(fe)]  fori  =  1,2,  and  if 

p{t)  =  Aipi(i)  -1-  X2P2{i),  and 

q{t)  =  Xiqi{t)  +  X2q2{t), 

where  Xi  >  0,  A2  >  0,  then  Lx  =  0  is  disconjugate  on  [p(o),  £r(6)]. 
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Proof.  Suppose  LiX  =  0  is  disconjugate  on  [p(a),  tr(&)]  for  i  =  1,2.  Then  the  quadratic 
functionals  P'1  and  P2  are  positive  definite  on  A.  Then  for  u  €  A,  we  have 

fa(b) 

■'<■(») 

f<T(b) 

Jpifl) 

=  /  it)f  -  qi{t)u  (t)]  Vt 

Jp{a) 

pa{h) 

+  / 

Jp{a) 

=  PiU  4-  p2'^- 

Therefore,  P  is  positive  definite  on  A,  and  hence  La:  =  0  is  disconjugate  on  [p(a),  17(6)]. 

□ 

We  summarize  some  of  the  major  results  in  the  following  theorem. 

Theorem  63  (Reid  Roundabout  Theorem).  The  following  are  equivalent: 

(i)  Lx  =  0  is  disconjugate  on  [p(a),a(i>)]; 

(a)  Lx  =  0  has  a  positive  solution  on  [p{a),o’{b)]; 

(Hi)  The  quadratic  functional  P  is  positive  definite  on  A; 

(iv)  the  Riccati  differential  inequality  Rz  <  0  has  a  solution  on  [p(a),6]. 

Proof.  By  Theorem  51,  (i)  and  (ii)  are  equivalent.  By  Theorem  60,  (i)  and  (iii)  are 
equivalent.  By  Theorem  56,  (ii)  implies  (iv).  It  remains  to  show  that  (iv)  implies  (i). 
So,  assume  Rz  <0  has  a  solution,  z  on  [p(a),  6],  and  let 

w{t)  ;=  Rz{t)  for  t  E.  [a,  b]. 

If  p{a)  <  a,  let  w(p(a))  =  0,  and  if  a{b)  >  b,  let  w{(7{b))  =  0.  Then  z  is  a  solution  of 
the  Riccati  equation 

Z^(t)  +  (n(t)  -  W(t))  +  ,  =  0 

-rygy)  \  jj-r  ^  u{t)zP{t) 

on  [p(a),6].  This  implies  that  the  self-adjoint  dynamic  equation 

{p{t)x^)^  +  {q{t)  -  w{t))x  =  0 
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is  disconjugate  on  [p{a),a{b)].  But 

q{t)  -  w{t)  >  q{t) 

on  [p{a),a{b)],  so  by  Theorem  61,  Lx  =  0  is  disconjugate  on  [p{a),cr{b)],  and  the 
proof  is  complete.  D 

Theorem  64  (Wintner’s  Theorem).  Assume  supT  =  oo,  a  G  T,  u{t)  >  K  >  0, 
and  0  <  p{t)  <  M  for  some  constants  K  and  M.  If 


poo 

/  q{t)  Vt  =  oo, 
J  a 


then  (2.1)  is  oscillatory. 

Proof.  By  way  of  contradiction  assume  the  hypotheses  of  the  theorem  hold,  but  (2.1) 
is  nonoscillatory.  Then  by  Lemma  49  there  is  some  to  €  T  such  that  (2.1)  has  a 
positive  solution  on  [to,  oo)-  Define  z  :  T  — »  M  by 

p{t)x^{t) 

■ 

Then  by  Theorem  56,  z  is  a  solution  of  the  Riccati  equation  Rz  =  0,  and  z  G  D/j,  so 
pP{t)  +  u{t)zP{t)  >  0,  and  therefore, 

^(t)  >  -€i!l  >  .E 

!/(«)-  K' 

If  we  integrate  both  sides  of  the  Riccati  equation,  we  see  that 

{z»{t)f 


poo  poo 

(O-z(io)  =  -f  9WV<-y^ 

poo 

<  -  q{t)  Vt. 

Jtn 


Vt 


□ 


This  implies  that  limf_,oo  z{t)  —  —  oo,  which  is  a  contradiction. 

Theorem  64  is  a  special  case  of  the  Leighton- Wintner  Theorem,  which  we  now 
state. 

Theorem  65  (Leighton-Wintner  Theorem).  Assume  supT  =  oo,  o  G  T,  and 
p{t)  >  0  for  all  t  G  T.  If 

poo  1  poo 

then  (2.1)  is  oscillatory  on  [a,oo). 
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(2.1)  is  nonoscillatory.  Then  there  is 
[to,  oo).  Then  by  Theorem  54,  there  is 


At  <  oo 

generality,  we  may  take  v{t)  >  0  on 

[&,ooj. 

Define  z{t)  by  the  Riccati  substitution.  Then  by  similar  reasoning  to  the  proof  of 
Theorem  64,  we  have 

z{t)  -  z{b) 

so  limt_^oo  z{t)  =  —  oo.  Then  we  can  find  ti  G  T  such  that  ti  >  b,  and  z{t)  <  0  on 
[ti,oo).  Then  as  2:  =  p{t)v^{t)/v{t),  v  is  positive  and  decreasing  on  [ti,oo),  so 


Proof.  By  way  of  contradiction,  assume  that 
some  to  eT  such  that  (2.1)  is  disconjugate  on 
a  dominant  solution,  v{t)  such  that 

1 

Jb  p{t)v{t)v<’{t) 

for  6  G  T,  sufficiently  large.  Without  loss  of 


1 

p{t)v{t)v^  {t) 


At 


1  f°°  _1_ 

“  v^ih)  Jtx  P{t) 
=  00, 


which  is  a  contradiction. 


□ 
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Chapter  3 

Linear  Dynamic  Equations  in 
Factored  Form 

3.1  Linear  Equations  with  Delta-Derivatives 

3.1.1  The  General  Case 

In  this  section,  we  look  at  solution  techniques  for  linear  equations  containing  delta- 
derivatives  which  can  be  written  in  factored  form.  Fix  n  G  N,  and  for  1  <  i  <  n  let 
Oj  :  T  — >  IR  be  rd-continuous.  Furthermore,  for  1  <  i  <  n,  assume  ai(t)  ^  0  for  any 
t  €  T.  Define  the  operators  1  <  i  <  n  by 

For  1  <  i  <  n,  we  take  the  domain  of  Da^  to  be  the  family  of  all  delta-differentiable 
functions.  Clearly  each  Da^  is  a  linear  operator. 

Now  consider  the  dynamic  equation 

Ly  =  0  where  Ly  =  -  Ai) (Da^  “  ^2)  •••  (jDa„  -  A„)y.  (3.1) 

Here  Aj,  1  <  i  <  n  are  (possibly  complex)  constants.  Note  that  the  order  in  which  the 
factors  are  written  down  is  important,  as  the  factors  do  not  necessarily  commute  with 
one  another.  We  define  the  domain  of  the  operator  L,  which  we  will  denote  by  D,  to 
be  the  set  of  all  functions,  y  :  T  — »  E,  such  that  {Da^  —  M)  •  •  •  {Dan—\n)y  '■  T'*”  ^  ^  E 
is  defined,  and  is  delta-differentiable.  In  the  case  n  =  1,  we  understand  this  to  mean 
that  our  domain  is  just  the  set  of  all  delta-differentiable  functions.  We  say  that  y  €  D 
is  a  solution  of  (3.1)  provided  Ly{t)  =  0  for  all  t  €  T*^". 

Definition  66.  We  say  the  dynamic  equation  (3.1)  is  regressive  provided  ^  e  7^  for 
1  <  i  <  n. 
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Theorem  67.  The  dynamic  equation  (3.1)  is  equivalent  to  the  system 


where 


and 


=  A{t)x, 


Xi 

X2 

Xfi 


1 

On  (^) 

*^71— 1 
1  (^) 


0 

0 


0 

1 

^'n— 1(^) 


^2 

a2{t) 

0 


0 

0 

0 

1 

02  (t) 
Ai 


(3.2) 


Proof.  To  see  that  this  system  is  actually  equivalent  to  om  original  equation,  suppose 
y  is  a  solution  of  (3.1),  and  let 


xi  =  y 

X2  =  (.^On  ^n)y 

Xn  ~  i^a.2  "^2)  •  •  •  {T)an  ^n)y- 


In  other  words, 

xi  =  y,  and  Xi+i  =  -  Xn+i-i)xu  l<i<n-l. 

Then  for  1  <  i  <  n  —  1,  we  have 

2^1+1  —  {Dan+i-i  ~  ^n+l-i)Xi 

—  On+1— -^n+l— 

Solving  for  xf',  we  see  that 

A  -^n+l— i  ,  1 

Xi  = - Xi  H - a:i+i. 
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Furthermore,  by  our  dynamic  equation,  we  have 

Ly  =  (^Dai  ^l)^n  ~  0- 

Solving  for  x^,  we  get 

A  -^1 

<  =  —Xn, 
fll 

and  we  see  that  a;  is  a  solution  of  (3.2). 

Conversely,  suppose  a:  is  a  solution  of  (3.2).  We  claim  that  y  ;=  xi  is  a  solution 
of  (3.1).  Note  that 


and  thus 

Now  we  also  have  that 


A 

<  =  —Xn, 
Oi 


{Dai  ^l)Xn  —  0- 


A  _  ^2 

Xn—l  —  Xn-1  "h  Xn- 
0,2  02 


Solving  this  for  a:„,  we  get 

Xn  ~  ^^{fyXn—x  ^2Xn—l  {Da^  ^2')Xn—l, 

and  therefore 

{Dai  ^l){Da2  ^2)Xn—l  —  0- 

Continuing  in  this  fashion,  we  get 

{Dai  ^l){Da2  -^2)  •  •  •  {Dan  ^n)Xl  =  0, 


which  is  the  desired  result.  □ 

Theorem  68.  If  (3.1)  is  regressive,  to  &  T ,  yi  E  M.,  1  <  i  <  n,  then  the  initial  value 
problem 

Ly  =  0, 

y{to)  =  yu  {Dan-K)y{to)  =  y2,  {Da2-)^2)  ■--{Dan-K)y{to)  =  yn  (3.3) 

has  a  unique  solution. 

Proof  Suppose  (3.1)  is  regressive,  and  consider  the  equivalent  system  (3.2).  By  [6, 
Theorem  5.8],  it  suffices  to  show  that  the  matrix-valued  function  Aft)  is  rd-continuous 
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and  regressive.  Recall  that  A{t)  was  defined  by 


^  ^n— 1 _ 1 

Ctn—l(t')  ®n— 1(^) 


A{t)  := 


X2  i 

«2(t)  a2(^) 

Xi 


Then,  as  each  function,  Oj,  1  <  i  <  n,  was  assumed  to  be  rd-continuous,  we  have 
that  A{t)  is  rd-continuous  as  well.  To  show  A{t)  is  regressive,  we  must  show  that 
det(/  +  fi{t)A{t))  0  for  any  t  €  T.  We  get 


det(/  +  fi{t)A{t)) 


1  +  M^)dr( 


=  det 


1  + 


^0, 

since  (3.1)  is  regressive.  Thus  the  IVP  (3.1),  (3.3)  has  a  unique  solution.  □ 

Now  that  we  have  established  that  solutions  of  IVPs  exist  and  are  unique,  we  turn 
our  attention  to  actually  finding  these  solutions,  and  examining  their  properties. 

Definition  69.  For  1  <  i  <  n,  let  yj  €  D.  We  define  the  Wronskian  associated  with 

Ly  =  0,  W  =  W(yi,...,yn):T«”"^  ^Rby 


yi  2/2  •  •  •  Vn 

DanVi  Da„y2  •  •  •  Da„yn 

W{yi,...,yn)  =det  D  dn—l  (D  (In  yi)  D  1  {D  (Xn  3/2)  ■■  D„_,{D  dn  2/n) 

■Da2(' •  *  (-^anyi))  Da2{' *  •  i^anV^))  •••  'Oa2(- *  •  (‘^anyn))_ 
Adding  (nonzero)  multiples  of  one  row  to  another  row  does  not  change  the  deter- 
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minant  of  a  matrix,  so  we  may  write  the  Wronskian  in  the  equivalent  form 

yi  • • •  Vn 

(Don  ~  ^n)yi  ■  •  •  i^an  ~  ^n)yn 

Vr  =  det 

{Da2  ~  -^2)  •  •  •  {Dan  ~  ^n)yi  ■  ■  •  (-^02  ~  -^2)  •  •  •  {Don  ~  ^n)yn 

We  then  get  the  following  theorem. 

Theorem  70.  If  (3.1)  is  regressive,  and yi,. . .  ,yn  are  solutions  of  (3.1),  then  either 
(i)  W{yi,...,yn)  =  0,  or 
(a)  W{yi, . .  .,yn){i)  ^  0  for  any  t  e  T. 

The  second  case  occurs  if  and  only  if  the  functions  yi,  1  <i  <n  are  linearly  indepen¬ 
dent  on  T. 

Proof  Suppose  (3.1)  is  regressive  and  yi, . . .  ,y„  are  solutions  of  (3.1).  Define  the  n 
vector  valued  functions,  xi,...,Xn  by 

yi 

{Dan  ~  ^n)yi 

Xi  =  . 

{Da2  ~  ^2)  •  •  •  {Da„  —  ^n)yi 

Then  we  see  that  for  1  <  i  <  n,  Xj  is  a  solution  of  the  system  (3.2).  Next,  define  the 
column  matrix,  X  by 

X  =  [Xi  X2  . . .  x„]  . 

Then  X  solves  the  matrix  dynamic  equation  X^  =  A{t)X.  Furthermore,  we  see  that 

W{yi,...,yn)  =  detX. 

Since  (3.1)  is  regressive,  A{t)  is  regressive,  and  we  may  apply  a  theorem  developed 
by  Cormani  in  [8]  to  conclude  that 

=  q{t)W,  (3.4) 

where  q  eTZ.  In  her  paper,  Cormani  describes  the  function  q  in  detail.  Here,  however, 
we  are  only  interested  in  the  fact  that  qE'R,  and  since  the  expression  describing  q  is 
quite  complicated,  we  choose  to  omit  it. 

Now,  since  q  €  TZ,  we  can  solve  equation  (3.4),  and  we  see  that 


W{t)  =  e,{t,to)W{to). 
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The  exponential  function  eq{t,to)  is  never  zero,  and  we  see  that  if  W{tQ)  =  0,  then 
W  =  0,  and  if  W{to)  ^  0,  then  W{t)^0  for  any  teT. 

It  remains  to  show  that  W(t)  ^  0  for  any  t  6  T  if  and  only  if  yi, . . . ,  yn  are  linearly 
independent. 

Suppose 

ClVl  +  C2y2  +  .  •  .  +  CnVn  =  0, 
where  Ci,  1  <  i  <  n  are  constants.  Then 

Cl{Da„  —  An)yi  +  C2{Dan  ~  An)y2  +  •  .  •  +  Cn{Da„  —  ^n)yn  —  0 

Cl (£>02  “  A2)  .  .  .  (£>a„  “  A„)yi  +  .  .  .  +  Cn(Da2  —  A2)  .  .  .  (Da^  —  An)yn  =  0, 
and  we  see  that 

Cl 

C2 

c  := 

Cn 

satisfies 

X{t)c  =  0, 

where  X  was  defined  in  the  first  part  of  this  proof.  This  holds  for  any  t  6  T**”  ^ ,  so 
we  may  evaluate  X  at  an  arbitrary  point,  ti  G  T  ,  and  we  get 

X(ti)c  =  0. 

Now  recall  that  H^(yi, . . .  ,y„)  =  det  X.  Hence  if  VF(ti)  7^  0,  then  X{ti)  is  invertible, 
and  we  see  that  c  =  0,  which  implies  that  yi, . . . ,  y^  are  linearly  independent.  If,  on 
the  other  hand,  W  =  Q,  then  X{ti)  is  not  invertible,  and  the  system  X{ti)c  =  0  has  a 
nontrivial  solution.  In  other  words,  there  are  constants,  ci, . . .  ,c„,  not  all  zero,  such 
that 

y  :=  ciyi  +  C2y2  + . . .  +  4yn 

satisfies  y(ti)  =  0.  Therefore,  by  the  linearity  of  the  operators  £>oj,  I  <n, 

y(ti)  =  {Da„  -  An)y(ti)  =  • . .  =  (£>02  -  A2)  •  •  •  (£>a„  -  A„)y(ti)  =  0. 

Furthermore,  since  yi, . .  .yn  are  solutions  of  (3.1),  y  is  also  a  solution  of  (3.1).  By 
Theorem  68,  however,  solutions  of  initial  value  problems  are  unique,  and  we  see  that 
this  implies 

y  =  ciyi  +  C2y2  +  . . .  +  CnVu  =  0. 

As  the  constants  q  were  not  all  zero,  this  proves  that  yi, . . . ,  y„  are  linearly  dependent 
on  T.  n 
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The  solutions  of  (3.1)  turn  out  to  be  somewhat  complicated.  To  simplify  things, 
we  introduce  the  following  notation.  For  i<j<n,  define  the  doubly- 

indexed  family  of  functions  {vij}  by 

yi,iit)  =  e^{t,to), 


and  for  i  +  1  <  j  <  n, 

y,At)  =  <=^(Mo)£  As. 

Speaking  informally,  the  first  subscript  tells  you  which  exponential  function  to 
start  with,  then  you  proceed  using  the  recursive  formula  described  above.  The  second 
subscript  tells  you  when  to  stop.  Note  also  that  the  expression  niay  be 

rewritten  as  follows: 

I _ = 

Ai  =  0. 


1 

.  ai{t) 


Occasionally,  we  will  make  use  of  this  simplification. 
Example  71.  Suppose  n  =  3.  Then 

y3,3  = 


03 


and 


2/1,3  = 

®3 


Theorem  72.  Suppose  (3.1)  is  regressive.  Then  {yi,n}f=i  linearly-independent 
solutions  of  (3.1)  are  ,  and  a  general  solution  of  (3.1)  is  given  by 


^  ''j  Oliyi,ni 


i=l 


where  Oj  are  constants,  1  <i  <n. 

Proof.  As  (3.1)  is  regressive,  is  well  defined  for  1  <  i  <  n,  and  thus  the  functions 

“t 

{2/i,n}r=i  exist.  It  is  relatively  straightforward  to  demonstrate  that  these  functions 
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are  in  D.  We  will  prove  the  remainder  of  the  statement  by  induction  on  n.  (Note, 
the  proof  is  by  induction  on  n,  the  order  of  the  equation,  not  by  induction  on  i.) 
Base  case,  n  =  1:  Then  our  dynamic  equation  has  the  form 


{Day,  -  Ai)y  =  0. 


Since  we  only  have  one  (nontrivial)  function,  linear  independence  is  clear.  It  remains 
to  show  that  yi^i  satisfies  the  dynamic  equation.  Substitution  gives 

{Day  -  \\)y\,\{t)  =  ai(t)yu(t)-Aiyi,i(t) 

=  ai  (t)e^  {t,  to)  -  Aie^  {t,  to) 

oi 


—  oi(t)  /,\^hx(^i^o)  Aie^(i,  to) 

flj  (t )  ay  “1 

—  Aie^(t,  to)  Ai  e a,  (t,  to) 

oi  ay 


=  0, 


and  we  have  successfully  established  the  base  case. 

Induction  step:  Now,  assume  that  linearly  independent  solutions 


of 


{Day  -  Xi){Dayy  -  As)  .  .  .  (I>a„_i  “  An-l)y  =  0, 


and  consider  {yi,n}?=i-  We  will  first  show  that  these  functions  are  solutions  of 


Ly  =  {Day  -  Xi){Da,  -  As) . . .  (Pa„  -  Xn)y  =  0. 


Note  that 

Ly  =  {Day  —  Ai)  .  .  .  {Don-y  ~  A^-l)  [{Da„  —  Aji)y] . 

Now,  our  functions  take  on  different  forms  depending  on  whether  i  =  n,  ox  1  <  i  < 
n  —  1.  Therefore,  we  must  consider  these  cases  separately. 

Case  —  Then,  ignoring  the  argument,  t, 


Lyn,n  — 


{Day-Xl). 

{Day-Xi). 

{Day-Xi). 

{Day  -  Al)  . 

{Day-X,). 


•  {Dun-y  An-l)  [(T^On  Xn)yn,n] 

•  {Don-y  ~  A„_i)  [Unyn.n  ~  Any„_„] 

•  {.Da„-y  Aji— i)  An  A^e A 


•  iPon-y  A„_l) 

.  (C.,-.  -  A._i)(0) 


_  \ 

Clrt  ^  An  An.  6  An 
dr) 
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Case  2,  l<i<n  —  1:  Then  yi,n  has  the  form 


2/i,n(i)  -  (a„(s)  +//(s)A„)  ^e^(^>'to)yi,n-i{s)  As 


and  hence 


Vtnit)  = 


Aji 


(t 


,to)  f  (-7^ 

Jto  \fln(s)  + 


^  \  /'c'l  J.  /  eQ^(s,to)yi,n-i(s)  As 


)u(s)A„ 


+e la.  (^>  ^o) 


1 


ar,{t)+n{t)\n) 

//(s)A 


(^)  ^o)yt,n— 1  (^) 


n(t)  +/x(t)An~\  _ 1 

/x(t)An, 


On  (t) 

An 


«,.w + ( I 


Cljx  (t) 


^n(^) 

Therefore,  we  see  that 


O^.C*)**'”**’ 

('^nyi,n(^)  "1“  2/i,n~l(^))  * 


1 


{Da^-\n)vyn{t)  =  “ '^nyi,nW 

=  dnif)  TTT  ('^nyijnC^)  “I"  yi,n-l(^))  ^nVi^nify 

an{t) 

“  '^nyi,n(^)  "h  yi,n— 1(^)  '^nyijTiC^) 

=  yi,n-l(i)) 


and  finally,  we  have 

Lyi^n  =  (^ai  “■  ^l)  •  •  •  {Dan-1  “  [(-^an  ““  ^n)yi,n] 

=  {Dai  '^l)  *  *  *  {Dan-1  ~  ^n-l)yi,n-l 

=  0, 

and  we  see  that  the  dynamic  equation  (3.1)  is  satisfied.  It  remains  to  show  that  the 
functions  {yi^n}i=.i  are  linearly  independent.  To  see  this,  consider  the  Wronskian, 
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W^(yn,n,---,2/i,n)-  We  have 
W  =  det 


yn,n 

i^an  n  )yn  ,n 


3/n-l,n 

{^an  “^n)3/n-l,n 


2/1, n 

(■C^an  ~^n)yiyn 


^{Da2~^2)''-{Dan^^n)yn,n  (-^02 “'^2)-*‘(^an  ^n)2/n~l,n  {^0.2  ^2)'*-(-C^an  An)2/l,n 

Now,  by  previous  calculations,  we  know  that  {Da^—\n)yn,n  =  0,  and  for  1  <  i  <  n— 1, 
(£)o„  -  K)yi,n  =  yi,n-i-  Using  these  facts  to  simplify  our  matrix,  we  get 


M^  =  det 


Vn.n 

0 


2/Ti~l,n 
2/n-l,n— 1 


2/1, n 
2/1, n-1 


0  (Z)a2^^2)‘'‘(-^a„_i  ~^n— l)2/n~l,n  — 1  •••  (-^02  ^2)***(-^a^_l  ^n— l)2/l,n— 1  J 

We  now  expand  this  determinant  about  the  first  column,  to  get 

2/n  — l,n  — 1  “•  2/1, n—1 

(£)a2“'^2)***(-^^ayi_>l  — An-l)2/Ti-l,n— 1  •••  (^02  “■A2)...(X)a„_l  ~‘An-l)2/l,n-l 

Looking  closer,  we  see  that  what  we  really  have  is 

W{yn,n,  •  •  •  ,  3/1, n)  =  3/n,nW (yn-l,n-l,  •  •  • ,  3/1, n-l)- 


W  =  y„,n  det 


Now  W(y„_i,n-i,  •  •  •  ,3/i,n-i)  is  never  zero,  since,  by  our  induction  hypothesis,  we 
have  assumed  these  functions  are  linearly  independent.  And,  as  exponential  func¬ 
tions  are  never  zero,  we  have  that  yn,n  is  also  never  zero.  Thus  we  can  conclude 
that  . . .  ,yi,n)  is  never  zero,  and  hence  the  functions  {yi,n}i=i  are  linearly 

independent.  By  induction,  the  proof  is  complete.  □ 


3.1.2  Special  Cases 

In  this  section,  we  will  look  at  how  this  theory  simplifies  in  the  case  where  all  of  the 
functions,  Oj,  1  <  i  <  n  are  the  same.  We  will  denote  this  common  function  simply 
by  a.  Formally,  let  a  :  T  ^  M  be  rd-continuous,  and  assume  that  a{t)  ^  0  for  any 
t  €  T.  Define  the  operator.  Da,  by 

[Day]{t)  :=  a{t)y^{t). 

Clearly,  Dq  is  a  linear  operator.  Now  consider  the  dynamic  equation 

Liw  =  0,  where  Liw  =  CnD^w  -I-  Cn-iD2~^w  -j - 1-  ciDaW  +  CqW,  (3.5) 

and  Ci  is  constant,  1  <  i  <  n,  c„  ^  0.  As  equation  (3.5)  can  be  multiplied  by  a  nonzero 
constant  without  changing  the  solutions,  we  will  usually  take  c„  =  1.  Here  we  take 
our  domain  Dl,  to  be  the  set  of  functions  lu  :  T  M  such  that  D2~^w  '■  T'*”  ^  —*  R 
is  defined  and  is  delta-differentiable,  and  we  say  that  a  function  w  €  Di,  is  a  solution 
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of  (3.5)  provided  Liw{t)  =  0  for  all  t  G  T''".  This  is  consistent  with  the  domain  and 
definition  of  solution  which  we  used  in  the  general  case. 

Definition  73.  The  characteristic  polynomial  associated  with  L\w,  which  we  will 
denote  by  p,  is  given  by 

p(A)  =  c„A”  +  Cfi-iX^  ^  +  •  •  •  +  CiA  +  cq. 


The  equation 


p{X)  =  0 


is  called  the  characteristic  equation. 


Definition  74.  Denote  the  n  roots  (including  multiplicity)  of  the  characteristic  equa¬ 
tion,  p(A)  =  0,  by  At,  1  <  i  <  n.  We  say  the  dynamic  equation  (3.5)  is  regressive 
provided  ^  G  7^  for  1  <  i  <  n. 

This  definition  agrees  with  the  earlier  definition  of  regressivity  of  (3.1).  Note, 
then,  that  equation  (3.5),  with  c„  =  1,  can  be  written  in  factored  form: 


Liw  = 


i=\ 


w  =  0. 


(3.6) 


It  is  easy  to  show  that  these  factors  commute  with  one  another. 

Although  the  results  presented  below  are  valid  for  any  rd-continuous  function,  a, 
without  zeros,  there  are  two  specific  choices  of  a  which  merit  further  comment.  First, 
if  we  choose  a{t)  =  1,  then  (3.5)  is  just  a  constant  coefficient  dynamic  equation. 
Second,  if  a(t)  =  t,  then  (3.5)  is  an  Euler-Cauchy  dynamic  equation.  These  equations 
have  been  studied  extensively  by  Bohner  and  Akin-Bohner  in  [7,  2],  and,  although 
our  notation  is  slightly  different,  the  results  in  this  section  should  be  considered  to 
be  a  generalization  of  their  work. 

One  of  the  drawbacks  of  our  definition  of  regressivity  is  that  the  roots  of  the 
characteristic  equation  must  be  known  before  it  can  be  determined  whether  or  not 
the  equation  is  regressive.  In  the  following  lemma,  we  present  an  equivalent  condition 
for  regressivity  which  does  not  require  the  roots  of  the  characteristic  equation  to  be 
determined  explicitly. 


Lemma  75.  Equation  (3.5)  is  regressive  if  and  only  if 


^  Ci(a(t))»(-)u(t))”"V0  for  all  t  G  T. 

i=0 

The  proof  essentially  mirrors  the  proof  presented  by  Bohner  and  Akin-Bohner  for 
Euler  equations  [7,  Theorem  2.15],  and  we  do  not  include  it  here. 
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Theorem  76.  Let  m  be  the  number  of  distinct  roots  of  the  characteristic  equation 
associated  with  LiW  =  0.  Denote  these  roots  by  ri,  1  <  i  <  m.  For  1  <  i  <  m,  let 
m{ri)  denote  the  multiplicity  of  the  root,  r^.  Then,  for  1  <i  <m,  defined 

by 


where 


Z,,,(t)  =  l,  and  for 

are  m{ri)  linearly  independent  solutions  of  L\w  =  0  associated  with  r^. 

Moreover,  if  aij  are  arbitrary  constants,  then  a  general  solution  of  LiW  =  0  is 
given  by 

m  m{ri) 

1=1  j=l 

Before  embarking  on  the  proof  of  this  theorem,  we  will  introduce  some  notation 
that  will  be  needed,  and  we  will  establish  a  lemma,  which  will  also  be  useful  in  the 
proof. 

Recall  that  our  dynamic  equation  can  be  written  in  factored  form.  (See  equation 
(3.6).)  Although  the  factors  commute,  for  the  time  being,  we  will  fix  the  order  of  the 
factors,  and  write  the  equation  as 


{Da  -  Xl){Da  -  As)  .  .  .  {Da  -  Xn)w  =  0. 

Note  that  although  the  constants  Aj  need  not  be  distinct,  we  have  that  for  1  <  i  <  n, 
Ai  =  rj  for  some  l<j<m.  Now  for  1  <  i  <  m,  1  <  A:  <  n,  define  the  set  (pk{ri)  by 

Tk{ri)  =  {j  ■  Xj  =  ri,l<j  <  k}, 

and  let 

'mk{ri)  =  \^k{ri)\, 

where  |  •  |  denotes  the  order  of  the  set.  In  other  words,  mk{ri)  is  the  number  of  times 
ri  appears  in  the  first  k  factors,  and  mn{ri)  =  m(ri).  Furthermore,  for  1  <  A;  <  n, 
define 

m  mk{ri) 

w't  =  U  U  {“’«>■ 

j=l 

(If  mk{ri)  =  0  for  some  i,  we  understand  ^e  empty  for  that  value  of  i.) 

Note  that  contains  exactly  k  elements,  and  that  Wk  C  Wfc+i  for  1  <  A:  <  n  —  1. 
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Lemma  77.  For  \  <k  <n  —  l,  if  Wij  €  Wk,  then 

belongs  to  the  span  ofWk+i,  denoted  by  S'fVFfc+i]. 

Proof.  Fix  k  with  1  <  A:  <  n  —  1.  We  must  show  that  Mjj.fc,  as  defined  in  the  statement 
of  the  theorem,  is  in  S'[Wfc+i]  for  all  Wij  €  Wk-  That  is,  we  must  show  Uij^k  €  S'fWfe+i] 
whenever  1  <  i  <  m,  and  1  <  j  <  mk{ri).  Not  that  if  mk{ri)  =  0  for  some  i,  then 
the  statement  is  vacuous.  Therefore,  fix  i  with  1  <  i  <  m,  and  assume  mk{ri)  >  1. 
We  now  proceed  by  induction  on  j. 

Base  case,  j  =  1:  If  j  =  1,  then 

Wij{t)  =  Wi^i{t)  =  ea(i,to)^i.i(i)  =  ea(t,to)- 

So,  we  get 

Ui,j,k{t)  =  e^{t,to)  £  0(5)  +  ^(5) Afc+i 

Now  this  expression  simplifies  differently,  depending  on  whether  or  not  A^+i  =  n,  so 
we  consider  these  cases  separately. 

Case  1,  Ajfc+i  =  rj:  In  this  case,  it  is  clear  that  mk+iin)  >  2,  and  we  get 


u. 


jM  =  ea((.«o)  -  »«W  e  M'wi- 


Case  2,  A^+i  ^  rf.  Choose  q  with  1  <  g  <  m  such  that  A^+i  =  rg.  In  this  case, 
we  see  that  mk+\{ri)  >  1,  and  mk+i{rq)  >  1.  Furthermore, 


eri(s,to) 

a 

era  (5,  to) 


As 


=  era(t,io)  / 
®  Jtc 


1  eii(5,to) 


As 


a 

(t,  to)  [e^e^(s,to)]  As 


=  - era  (t,  to)  [eri©ra  (5,  to) 

rj  —  Tq®  '  I  a  J  s=to 

—  I  ^era(t,  to)  ^eri0ra  (t,  to)  1 

=  — —  eri(t,to)  -  era(t,to)l 

r<  —  L  “  “  J 


Ti-  iq 
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n-Vg 


Wi,lit)  -  Wg^l{t)] 


As  either  Case  1  or  Case  2  must  occur,  we  see  that  this  establishes  the  base  case. 

Induction  step:  Now  suppose  Uij^k  €  and  assume  Wjj+i  G  Wk-  (If 

Wij+i  ^  Wk,  there  is  nothing  to  show.)  Then,  we  must  show  Uij+i.fc  G  5'[Wfc+i]-  In 
this  case,  we  have 

Wij+i  (t)  =  ea  {t,  to)zij+i  (t) . 

Thus, 

Uij+i,k{t)  =  e^(t,to)£  o(s)  +  )u(s) Afe+i  io)zi,j+i  (s)  As 

As  before,  this  expression  simplifies  differently  depending  on  whether  or  not  A^+i  =  n, 
so  we  consider  these  cases  separately. 

Case  1,  Afc+i  =  n:  In  this  case, 

/■*  1 

Uij+i,k{t)  =  ea{t,to)  a(s)  + 

Now,  as  Wij+i  G  Wk,  we  have  mk{ri)  >  j  +  I,  and,  since  A^+i  =  Vi,  we  see  that 
rrik+iiri)  >j +  2.  Hence  Wij+2  G  Wk+i,  which  gives  Uij+i.fc  G  S'fHfc+i]. 

Case  2,  Afc+i  ^  ri:  In  this  case. 


tj+l,fc(f)  —  (^)  ^o)  / 

“  Ju 


‘  'TT~, — Tvi — e  is,to)en(s,to)zu+iis)  As 


Ti  -  \k 


'k+l  (^1  ^o)  / 

“  Jto 


Now,  integrating  by  parts,  we  get 


X  ^^k+l  (^?  ^o) 

’j  —  Afc-j-l  a  \  a  ^  o  / 

-  f  ■ 

Jto  a  J 


Now,  note  that 
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so 


4-m(*)  = 


1 


a{t)  +  ^i{t)ri 

Additionally,  we  get 

■2i,j+l(to)  —  0- 

Substituting  the  expression  for  into  our  expression  for  and  making  some 

other  simplifications,  we  get 


Ti  “  L  a  a 

”  /o  ^  ^  a(s)  +  ii{s)ri 

e!i{t,to)zij+i{t) 


ri  -  Afc+i  [  “ 


/  1 

eAfc+i(t,to)  /  -7-re^^(s,io)ea(s,to)^i,j(s)  As 

~  ^  a(s)  ®  a  “  J 


n  -  Afc+i 


to 

Wi,j+l{t) 


f*  1 

^  .  V. — e  ^(s,to)eri(s,to)^i,i(s)  As 


n  -  Afc+i 


[Wij+i{t) 


Ti  —  Afc+i 


!^i,j+l{t)  '^i,j,k{t)] 


By  our  induction  hypothesis,  Uij^k  €  5[M4+i],  so  we  see  that  «tj+i,ifc  € 

Thus,  by  induction,  we  have  shown  that  for  1  <  j  <  mk{ri),  Uij^k  €  S'[lTfc+i].  Since  i 
and  k  were  arbitrary,  this  completes  the  proof  of  the  lemma.  □ 

Now,  armed  with  this  lemma,  we  are  ready  to  prove  Theorem  76. 

Proof  of  Theorem  76.  By  Theorem  72,  we  know  that  the  n  linearly  independent  so¬ 
lutions  of  Liw  =  0  are  {j/i,n}"=i)  as  defined  in  the  previous  section.  Let  Y  denote  this 


set.  That  is,  put  Y  =  {yj,n}r=i-  Then,  what  we  are  really  trying  to  show  is  that 

S-IF]  =  5[W„], 
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Where,  as  before,  S'[V]  denotes  the  span  of  the  set  V.  Note  that  5[y]  is  finite 
dimensional.  Specifically,  since  the  elements  of  Y  are  linearly  independent  on  T, 
S'[y]  has  dimension  n.  We  also  know  that  W„  contains  n  elements,  and  therefore 
dim(6'[Wn])  <  n.  Thus  we  need  only  show  that  Y  C  S'[W„],  which  will  yield  S'fy]  = 
5(^4],  and,  furthermore,  will  establish  the  linear  independence  of  the  elements  of 
on  T. 

Actually,  we  are  going  to  show  more  than  Y  C  We  are  going  to  show  that 

for  1  <  i  <  n,  i  <  i  <  n,  yij  €  ^IWn].  Fix  i,  with  1  <  i  <  n.  We  now  proceed  by 
induction  on  j.  We  are  going  to  show  that  for  i  <  j  <n,  yi^j  €  Then,  since 

Wj  C  Wj+i  G...CWn, 


this  will  establish  the  desired  result 

Base  Case,  j  =  i:  We  must  show  that  yi,i  G  5'[Wi].  By  definition,  yi^i{t)  = 
e^{t,to).  But  Aj  =  rk  for  some  1  <  A;  <  m,  and,  clearly,  mi{rk)  >  1.  Hence 

a 

a  ^ 


SO  yi,i  G 

Induction  step:  Assume  yij  G  S[Wj].  We  must  show  that  yij+i  G  S'fWj+i].  We 
have 

pt 

yi.i+iji)  —  (^i^o)  /  /  \  ,  /  \  \  (^i^o)yi,i(^) 

Jto  ®  “ 

Now,  as  yij  G  -S'fW}],  we  can  write  it  as  a  linear  combination  of  elements  of  Wj: 

m  mj(rk) 
fc=l  p=l 


where  ak  p  are  constants.  Thus 


ft  I 

yij+iit)  =  e^{t,to)J^^  ^  «fc,p«^fe,p(s)  As. 


These  sums  are  finite,  so  we  may  move  them  outside  the  integral  without  any  trouble. 
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yielding 


ViJ+l  (^) 


TO  rnj{rk)  -f 

g  g  ak,p 


m  "rrijirk) 

EE  0!k,p'^k,p,j{t)  1 

k=l  p=l 


where  Uk,pj  is  as  defined  in  Lemma  77.  By  Lemma  77,  Uk,pj  6  5[Vl^+i],  and  therefore, 
yij+i  is  a  linear  combination  of  elements  of  Hence  Pij+i  €  S'fHj+i].  By 

induction,  then,  we  have  shown  that  for  i  ^  J  ^  n,  pij  €  ‘S'fVKj],  which  implies  that 
for  i  <  j  <  n,  pij  G  <S[W„].  Our  choice  of  i  was  arbitrary,  so  we  have  established 
that  for  1  <  i  <  n,  i  <  i  <  n,  Pij  G  S'fWn].  Hence  Y  C  S'fWn],  and,  by  our  earlier 
discussion,  the  proof  is  complete.  □ 


3.2  Linear  Equations  with  Nabla-Derivatives 

In  this  section,  we  examine  linear  equations  containing  nabla-derivatives  which  can 
be  written  in  factored  form.  As  one  would  expect,  the  results  are  analogous  to  the 
results  developed  above  for  delta-derivatives,  and  the  proofs  are  essentially  the  same. 
Therefore,  we  will  simply  state  the  theorems  here,  and  not  repeat  the  proofs. 

Fix  n  G  N,  and  for  1  <  i  <  n  let  6j  :  T  — >  R  be  Id-continuous.  Furthermore,  for 
1  <  i  <  n,  assume  bi{t)  ^  0  for  any  i  G  T.  Define  the  operators  Ai,.,  1  <  i  <  n  by 

WbiVW)  ■= 

For  1  <  i  <  n,  we  take  the  domain  of  A,,  to  be  the  family  of  all  nabla^differentiable 
functions.  Clearly  each  Ni,^  is  a  linear  operator. 

Now  consider  the  dynamic  equation 

My  =  0  where  My  =  -  \i){Nb^  -  A2) . . .  (A,„  -  A„)y,  (3.7) 

and  Ai,  1  <  i  <  n  are  (possibly  complex)  constants.  We  define  the  domain  of  the 
operator  M,  which  we  will  denote  by  Dm,  to  be  the  set  of  all  functions,  y  :  T  — »  R, 
such  that  (A,2  —  A2) . . .  (A,„  —  A„)y  :  T^n-i  — >  R  is  defined,  and  is  nabla-differentiable. 
In  the  case  n  =  1,  we  understand  this  to  mean  that  our  domain  is  just  the  set  of  all 
nabla-differentiable  functions.  We  say  that  y  G  Dm  is  a  solution  of  (3.7)  provided 
My{t)  =  0  for  all  t  G 

Definition  78.  We  say  the  dynamic  equation  (3.1)  is  regressive  provided  ^  G  TZu 
for  1  <  i  <  n. 
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Theorem  79.  The  dynamic  equation  (3.7)  is  equivalent  to  the  system 


where 


and 


=  B(t)x, 


Xi 

X2 


1 

bnit) 

■^n-1 

K-lii) 


0 

0 


0 

1 

bn-l(t) 


0 

0 


A2 

b2{t) 

0 


0 

1 

b2(t) 

bi(t). 


(3.8) 


Theorem  80.  If  (3.7)  is  regressive,  io  €  T,  y*  e  M,  1  <  *  <  n,  then  the  initial  value 
problem 

My  =  0, 

y(to)  =  yi,  (Nb„-K)y(to)  =  y2,  ,•••,  (Nb2->^2)---(Nb„-K)y(to)  =  yn  (3-Q) 

has  a  unique  solution. 


Definition  81.  For  1  <  i  <  n,  let  j/j  6  Dm-  We  define  the  Wronskian  associated 
with  My  =  0,  Wn  =  Wiv(yi, . . . ,  y„)  by 


Wjv(yi,  • .  • ,  yn)  =  det 


yi 

Nbnyi 

Nbn-l(^bnyi) 


y2 

Nbny2 

Nbn-l(Nbny2) 


yn 

■^6„yn 

^6n-i(-^6„yn) 


Nb2(..-(Nb„yi))  Nb^(- .  ■  (Nb„y2))  •••  Nb2(- ■  ■  (Nbnyn)) 


Adding  (nonzero)  multiples  of  one  row  to  another  row  does  not  change  the  deter¬ 
minant  of  a  matrix,  so  we  may  write  the  Wronskian  in  the  equivalent  form 

yi  •  •  •  Vn 

(^bn  ~  ^n)yi  •  •  •  (^b„  ~  ^n)yn 

(Nb,-X2)...(Nb,,-Xn)yi  ...  (Nb,-X2)...(Nb„-Xn)yn_ 


Wn  =  det 
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Theorem  82.  Ifyi,...,yn  o,re  solutions  of  (3.7),  then  either 

(i)  WNiyi,...,yn)  =  0,  or 

(ii) WN{yi,---,yn)¥^OforanyteT. 

The  second  case  occurs  if  and  only  if  the  functions  yi,  1  <  i  <  n  are  linearly  indepen¬ 
dent. 

Now,  introduce  the  following  notation.  For  l<i<n,  i<j<n,  define  the 
doubly-indexed  family  of  functions  {yij}  by 


and  for  i  -|- 1  <  <  n, 


Theorem  83.  Suppose  (3.7)  is  regressive.  Then  {j/i,n}"=i  linearly-independent 
solutions  of  (3.7),  and  a  general  solution  of  (3.7)  is  given  by 


n 

i=l 


where  are  constants ^  1  <  i  <n. 

As  before,  we  now  look  at  the  special  case  where  all  of  the  functions,  bi,  1  <  i  <  n 
are  the  same.  We  will  denote  this  common  function  simply  by  b.  Let  6  :  T  ^  M  be 
Id-continuous,  and  assume  that  b{t)  ^  0  for  any  t  G  T.  Define  the  operator,  Ai,,  by 

[Nby]{t)  :=  b{t)y'^{t). 

Clearly,  Nb  is  a  linear  operator.  Now  consider  the  dynamic  equation 

Miw  =  0,  where  Miw  =  CnN^w Cn-iN^~^w -\ - h  ciNbW cqw,  (3.10) 

and  Ci  is  constant,  1  <  i  <  n.  Again  we  will  usually  take  c„  =  1.  Here  we  take  our 
domain  Dmi  to  be  the  set  of  functions  u;  ;  T  ^  R  such  that  Nb~^w  :  T^n-i  — >  R  is 
defined  and  is  nabla-differentiable,  and  we  say  that  a  function  w  G  DjUi  is  a  solution 
of  (3.10)  provided  Miw{t)  =  0  for  all  t  6  T^n. 

Definition  84.  The  characteristic  polynomial  associated  with  Miw,  which  we  will 
denote  by  q,  is  given  by 

^(A)  =  CnA”  +  Cn-lX^  ^  +  •  •  •  d"  CiA  -|-  Co- 
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The  equation 


9(A)  =  0 


is  called  the  characteristic  equation. 


Definition  85.  Denote  the  n  roots  (including  multiplicity)  of  the  characteristic  equa¬ 
tion,  9(A)  =  0,  by  Ai,  1  <  i  <  n.  We  say  the  dynamic  equation  (3.10)  is  regressive 
provided  ^  6  Tlv  for  1  <  f  <  n. 

Note,  then,  that  equation  (3.10),  with  Cn  =  1,  can  be  written  in  factored  form: 


MiW  = 


n(^‘  -  •'<) 


ii;  =  0. 


Li=i  J 

It  is  easy  to  show  that  these  factors  commute  with  one  another. 
Lemma  86.  Equation  (3.10)  is  regressive  if  and  only  if 


(3.11) 


for  all  tel. 

t=0 

Theorem  87.  Let  m  be  the  number  of  distinct  roots  of  the  characteristic  equation 
associated  with  MiW  =  0.  Denote  these  roots  by  rj.  For  1  <  i  <  m,  let  m(ri)  denote 
the  multiplicity  of  the  root,  ri.  Then,  for  1  <i  <m,  defined  by 


Wij{t)  =  er^(t,to)Zij{t), 

where 


Zi,i{t)  =  1,  and  Zij+i 


Jto 


b{s)  —  i'{s)r. 


ZiJ 


(s)  Vs  for  1  <  j  <  m{ri)  —  1. 


are  m{ri)  linearly  independent  solutions,  of  Miw  =  0  associated  with  ri-  Moreover,  if 
aij  are  arbitrary  constants,  then  a  general  solution  of  Miw  =  0  is  given  by 


m  Tn(ri' 


i=l  j-l 


3.3  Equations  with  Both  Delta  and  Nabla  Deriva¬ 
tives 

Here,  we  turn  our  attention  to  linear  equations  containing  both  delta  and  nabla 
derivatives.  Things  become  quite  a  bit  more  complicated  when  these  two  types  of 
derivatives  interact,  so  here  we  will  only  look  at  second-order  dynamic  equations.  Let 


65 


oi  :  T  ^  R  be  continuous,  let  02  :  T  ^  R  be  rd-continuous,  and  assume  that  oi  and 
02  do  not  vanish.  Now  consider  the  dynamic  equation 


{Na,-Xi){Da,-X2)y  =  0.  (3.12) 

Definition  88.  We  say  (3.12)  is  regressive  provided  ^  is  i/-regressive  and  ^  is 
regressive. 

So,  suppose  (3.12)  is  regressive.  Now,  if  we  look  at  the  left-most  factor,  we  notice 
that  this  factor  annihilates  the  function  e^(t,  to)-  As  oi  is  continuous,  we  may  apply 

Theorem  31  to  see  that 

“  6  Ai 

ai  a^-Xifi 

and  this  second  function  is  annihilated  by  —  Ai).  We  then  theorize  that  if 

we  replace  the  factor  {Naj,  -  Ai)  by  the  factor  {Da^-x^^  -  Ai),  perhaps  we  will  get  a 
related  equation  which  contains  only  delta  derivatives.  This  new  equation  can  then 
be  solved  using  techniques  from  previous  sections.  We  formalize  this  result  in  the 
following  theorem. 

Theorem  89.  Suppose  (3.12)  is  regressive,  and  let  bi{t)  :=  Oi(t)  —  Ai/u(t).  If  y  is  a 
solution  of  the  dynamic  equation 

{Dt,-Xi)iDa,-X2)y  =  0,  (3.13) 

then  y  is  also  a  solution  of  (3.12). 

Proof.  Assume  (3.12)  is  regressive,  and  let  bi{t)  :=  ai{t)  —  Ai^(t).  Furthermore,  let 
y  be  a  solution  of  (3.13).  Then  by  Theorem  72,  y  has  the  form 

y(t)  =  Cie^{t,t(j)  +  C2e^ (t, to)  [  ©  (  /  \ O0^(^(to)6^(s, to)  As. 

“2  “2  Jto  /  “2  <>1 

Substituting  this  expression  into  (3.12),  we  see  that 

(iVa,-Ai)(D„,-A2)y(t)  =  (A„,-Ai)e^(t,to) 

=  {Nai  —  Xi)e^{t,to) 

“1 

=  0, 


and  the  proof  is  complete.  □ 

Note  that  we  are  not  asserting  that  these  equations  are  equivalent,  nor  do  we  make 
any  uniqueness  claims  regarding  the  solutions  of  initial  value  problems  associated  with 
(3.12).  These  issues  will  require  further  exploration  at  a  later  date. 
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Of  course,  there  is  no  particular  reason  the  factor  with  the  nabla  derivative  was 
written  first.  We  could  just  as  easily  consider  a  dynamic  equation  of  the  form 


(A,-Ai)(iV6,-A2)2/  =  0.  (3.14) 

In  this  case,  we  assume  bi  is  continuous,  62  is  Id-continuous  and  neither  bi  nor  62 
vanish,  and  we  say  that  (3.14)  is  regressive  if  ^  is  regressive  and  ^  is  i/-regressive. 
Again,  we  will  replace  the  left-most  factor  by  a  related  factor  containing  a  nabla 
derivative,  and  obtain  the  following  result. 

Theorem  90.  Suppose  (3.14)  is  regressive,  and  let  ai(t)  :=  l^{t)  +  If  y  is  a 

solution  of  the  dynamic  equation 

(iV„,-Ai)(iV6,-A2)y  =  0,  (3.15) 

then  y  is  also  a  solution  of  (3.14). 

The  proof  is  analogous  to  that  of  the  previous  theorem. 

Clearly,  the  results  contained  in  this  section  are  only  a  beginning.  We  plan  to 
continue  studying  the  mixed  derivative  case,  and  hope  to  include  broader  results  in 
future  publications. 

3.4  More  Equations  with  Mixed  Derivatives 

In  final  section  of  this  chapter,  we  turn  our  attention  to  an  equation  which  cannot 
itself  be  written  in  factored  form,  but  which  is  equivalent  to  one  which  can  be  written 
in  factored  form.  The  second-order  dynamic  equation 

Lx  =  0,  where  Lx  =  x^^  -I-  ax^  -I-  ^x,  and  a,  €  K,  (3.16) 

has  been  studied  extensively.  (See,  for  example,  [6,  Section  3.2,  3.3].)  Equation  (3.16) 
is  said  to  be  regressive  if  1  -  Oiix{t)  -|-/?/x^(f)  7^  0  for  t  e  T,  or  in  other  words,  if  -  a 
is  regressive. 

In  this  section,  we  look  at  second-order  linear  dynamic  equations  of  the  form 

Mx  =  0,  where  Mx  =  x^^  -f  ax^  +  (dx^  and  a,  /3  €  M.  (3.17) 

We  take  the  domain  of  this  operator  to  be  the  set  of  all  functions  x  :  T  — »  R  such  that 
X  is  A-differentiable  on  T'^,  x^  is  V-differentiable  on  T]^,  and  is  Id-continuous 
on  TJJ.  We  say  a  function  x  ;  T  M  is  a  solution  of  (3.17)  if  x  is  in  the  domain 
of  M,  and  Mx{t)  =  0  for  all  t  G  Equation  (3.17)  is  said  to  be  regressive  if 
1  —  au{t)  +  /du^it)  7^  0  for  t  G  T,  or  in  other  words,  if  a  —  /?!/  is  z/-regressive.  It  was 
demonstrated  in  Chapter  2  that  equation  (3.17)  can  be  written  in  self-adjoint  form, 
and  hence  solutions  of  initial  value  problems  for  this  equation  are  unique. 
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Lemma  91.  The  dynamic  equation  (3.16)  is  regressive  if  and  only  if  the  dynamic 
equation  (3.17)  is  regressive. 

Proof.  First  suppose  that  (3.16)  is  regressive.  We  then  wish  to  show  that  (3.17)  is 
regressive.  That  is,  we  want  to  show  that  a-  Pu  is  i/-regressive,  or  that 

1  —  u{t){a  —  ^^{t))  ^  0  for  t^T 

Case  1:  t  is  left-dense.  Then  i/(t)  =  0,  and  we  have 

1  —  iy{t){a  —  Pv{t))  =  1^0. 

Case  2:  t  is  left-scattered.  Then  (T{p{t))  =  t.  As  (3.16)  is  regressive,  we  have  that 
P/j,-  a  is  regressive.  Thus  by  Lemma  30,  -{Pp  -  aY  is  ^-regressive,  and  hence 

1  -  u{t){a  -  Pp{p{t)))  i-  0. 


Then  we  have 


1  -  u{t){a  -  Pu{t))  =  1  -  u{t){a  -  P{t  -  p{t))) 

=  1  -  u{t){a  -  P{a{p{t))  -  p{t))) 

=  l-iy{t){a- Pp{p{t))) 

^  0. 

Hence  a  —  Pu  is  i/-regressive,  and  thus  (3.17)  is  regressive. 

Conversely,  suppose  that  (3.17)  is  regressive.  This  time,  we  seek  to  show  that 
Pp.  —  a  is  regressive,  or  that 

1  -1-  p{t){Pp(t)  —  a)  ^0  for  t  G  T. 

Case  1:  <  is  right-dense.  Then  p{t)  =  0,  and  thus 

1  +  p{t){Pp{t)  -  a)  =  1  0. 

Case  2:  t  is  right-scattered.  Then  p{cr{t))  =  t.  Since  (3.17)  is  regressive,  we  have 
that  a  —  /3i/  is  t'-regressive,  and  thus  by  Lemma  30,  -{a-  P^Y  is  regressive.  This 
give  us  that 

1  + p{t){Pv{a{t))-a)i^i). 


1  +  p{t){Pp{t)  -  a) 


=  1  +  p{t){p{a{t)  - 1)  -  a) 

=  1  +  -  «) 
=  1  +  p{t){pu{<T{t))  -  a) 

+  0. 


Then 
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Hence  Pfi  —  a  is  regressive,  and  thus  (3.16)  is  regressive.  □ 

Theorem  92.  If  (3.16)  is  regressive  and  x  is  a  solution  of  (3.16),  then  (3.17)  is 
regressive  and  x  is  a  solution  of  (3.17). 

Proof.  Assume  (3.16)  is  regressive  and  a:  is  a  solution  of  (3.16).  The  fact  that  (3.17) 
is  regressive  follows  immediately  from  Lemma  91.  It  remains  to  show  that  a;  is  a 
solution  of  (3.17). 

Now,  as  a;  is  a  solution  of  Lx  =  0,  we  know  that  x^  is  A-differentiable,  and  hence 
continuous.  Furthermore,  since  a;  is  a  solution  of  Lx  =  0,  we  have  that 

a:^'^  =  —ax^  —  fix. 

In  other  words,  x^^  is  a  linear  combination  of  continuous  functions  and  is  therefore 
continuous.  Thus,  applying  Corollary  28  we  see  that 

Mx  =  x^^  +  ax'^  +  flx^ 

=  x^^^  +  ax^^  +  fix'’ 

=  {x^^  +  ax^  +  ^xY 
=  {LxY 
=  0. 


□ 

Corollary  93.  The  solution  of  (3.17)  is  determined  by  the  roots  of  the  characteristic 
equation  as  follows: 

(i)  If  —  4/3  >  0,  then  the  characteristic  equation  has  distinct  real  roots,  \\  and 
A2,  and  a  general  solution  of  (3.17)  is  given  by 

x{t)  =  cieAi(t,to)  +  O2ex^{t,to). 

(ii)  If  —  Afl  <  0,  then  the  characteristic  equation  has  complex  conjugate  roots 
p  ±  iq,  and  a  general  solution  of  (3.17)  is  given  by 

aj(t)  =  Ciep{t,  to)  ^0)  "b  ^2^p(t,  to)  ®1^9/n-/tp(t)  to). 

(Hi)  If  —  4f3  =  0,  then  the  characteristic  equation  has  a  repeated  real  root,  X,  and 
a  general  solution  of  (3.17)  is  given  by 

nt  ^ 

x{t)  =  ciex{t,to)  +  C2ex{t,to)  ^  As. 

Corollary  94.  The  regressive  equations  (3.16)  and  (3.17)  are  equivalent. 
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Proof.  As  demonstrated  in  Chapter  2,  equation  (3.17)  can  be  written  in  self-adjoint 
form,  and  therefore,  solutions  of  IVPs  for  (3.17)  are  unique.  Therefore,  although 
Theorem  92  only  establishes  one  direction  of  this  equivalence,  the  converse  follows 
from  uniqueness  of  solutions.  □ 


Chapter  4 

A  Self-Adjoint  Matrix  Equation  on 
a  Time  Scale 


In  this  chapter,  we  are  concerned  with  the  self-adjoint  matrix  equation  [P(t)A^]'^  + 
Q{t)X  =  0.  This  is  a  generalization  of  the  scalar  equation  studied  in  Chapter  2,  and 
many  of  the  results  we  obtain  are  analogous  to  results  obtained  in  the  scalar  case. 
The  added  complexity  of  the  matrix  case,  however,  does  provide  a  couple  of  surprises 
along  the  way. 

4.1  Preliminary  Results 

As  in  the  scalar  case,  contained  in  Chapter  2,  we  begin  by  establishing  some  results 
regarding  the  interaction  of  the  A  and  V  derivatives.  We  first  note  that  differentiation 
and  integration  of  matrices  is  defined  in  terms  of  differentiation  and  integration  of 
the  matrix  entries,  and  therefore  Theorem  26  and  its  corollaries  carry  through  in  the 
matrix  case  without  further  effort. 

The  definition  of  regressivity  (or  i/-regressivity),  however,  changes  slightly  in  the 
matrix  case,  so  extending  Lemma  30  and  Theorem  31  to  the  matrix  case  requires 
some  justification. 

Lemma  95.  Let  P  be  an  nx  n  matrix-valued  function  on  T.  Then  P  is  regressive 
if  and  only  if  —P^  is  u -regressive.  Similarly,  ifQisannxn  matrix-valued  function 
on  T,  then  Q  is  u-regressive  if  and  only  if  is  regressive. 

Proof.  We  will  only  prove  the  first  statement.  The  proof  of  the  second  statement  is 
similar. 

First,  assume  the  nxn  matrix-valued  function,  P,  is  regressive.  We  then  wish  to 
show  that  I  +  u{t){P^{t))  is  invertible. 

Case  1:  Fix  t  G  T^.  Then  p{t)  G  T,  and  as  P  is  regressive,  we  have  that 


I  -H  p{p{t))P{p(t)) 
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is  invertible.  Using  the  definition  of  /x(t), 

I  +  [(T{p{t))  -  p{t)]P{p{t)) 
is  invertible.  But  t  €  T>i,  so  (T{p{t))  —  t,  and  we  get 

I  +  [t-p{t)]PPit)  =  l  +  u{t)P^{t) 


is  invertible,  as  desired. 

Case  2:  Fix  t  e  A.  Then  t  is  left-dense  and  right-scattered,  so  u{t)  =  0.  Hence 

I  +  v{t)PP{t)  =  1  +  0PP{t)  =  /, 

which  is  invertible.  As  /  +  v{t)PP{t)  is  invertible  for  any  t  €  T,  we  see  that  —P^  is 
iz-regressive. 

Conversely,  suppose  — is  i/-regressive.  We  then  wish  to  show  that  I  +  p,{t)P{t) 
is  invertible. 

Case  1:  Fbc  t  6  T^.  Then  <j{t)  €  T,  and,  as  -F^  is  i/-regressive,  we  have  that 

I  +  u{a{t))PP{a{t)) 
is  invertible.  Using  the  definition  of 

I  +  [(T(i)  -  p{(T{t))]P{p{a{t))) 
is  invertible.  But  t  €  Tb,  so  p{(T{t))  =  t,  and  we  get 

I  +  [o-(t)  -  t]P{t)  =  I  +  p{t)P{t) 


is  invertible,  as  desired. 

Case  2:  Fix  t  £  B.  Then  t  is  right-dense  and  left-scattered,  so  p,{t)  =  0.  Hence 

I  +  fi{t)P{t)  =  I  +  0P{t)  =  I, 

which  is  invertible.  As  1  -f-  p{t)P{t)  is  invertible  for  any  t  €  T,  we  see  that  F  is 
regressive.  d 

Theorem  96  (Equivalence  of  delta  and  nabla  exponential  functions).  If  the 

nx  n  matrix-valued  function  P  is  continuous  and  regressive,  then 


ep(t,  to)  =  ee„(_pp)(t,to)- 

If  the  nx  n  matrix-valued  function  Q  is  continuous  and  u -regressive,  then 


&Q{t,to)  —  e0(_Q<7)(t,to)- 

Proof.  We  will  only  prove  the  first  statement.  The  proof  of  the  second  statement  is 
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similar.  Suppose  that  the  n  x  n  matrix-valued  function  P  is  continuous  and  regres¬ 
sive,  then  by  Lemma  30  we  have  that  — is  i/-regressive.  Furthermore,  since  P  is 
continuous,  — P^  is  Id-continuous.  Hence  — P^  €  IZv.  Then  asTZ^ism  Abelian  group 
under  0,,,  we  see  that  Qu{—P^)  €  TZu,  and  therefore  e©,^(_p,»)(t,to)  exists. 

To  complete  the  proof,  then,  it  suffices  to  show  that  ep{t,to)  solves  the  initial 
value  problem 

Y^  =  eA-p'’)y,  y{to)  =  i. 


Note  first  that 


ep(^0)*o)  =  I- 


Furthermore,  ep{t,to)  =  P{t)ep{t,to),  which  is  continuous.  Hence  by  Corollary  28, 
ep(i,  to)  =  to),  and  we  get 


ep{t,to)  =  ep'^{t,to) 

=  Pf’{t)e^p{t,to) 

=  P^it)  [ep{t,  to)  -  i'{t)ej{t,  to)]  • 


Rearranging  this  equation  gives 

[I  -1-  i/(t)P^(t)]  ep(t,  to)  =  P^(t)ep(t,  to), 


so 

4(t,to)  =  [I  +  u{t)PP{t)]-^P^{t)ep{t,to) 
=  ©i/(~P^)ep(t,  to), 


and  the  proof  is  complete. 


□ 


4.2  Abel’s  Formula  and  Reduction  of  Order 

Let  P  and  Q  be  Hermitian  nxn  matrix-valued  functions  on  a  time  scale  T.  Further¬ 
more,  assume  that  P  is  continuous,  Q  is  Id-continuous,  and  P(t)  is  invertible  for  all 
t  e  T.  Here  we  are  concerned  with  the  second-order  matrix  equation 

LX  =  0,  where  LX  =  [Pit)X^]^  +  Q{t)X.  (4.1) 

We  take  the  domain  of  L,  denoted  by  D,  to  be  the  set  of  all  n  x  n  matrix-valued 
functions  X  defined  on  T  such  that  X  is  delta-differentiable  on  T'',  X^  is  continuous, 
PX^  is  nabla-differentiable  on  T",  and  [PA^]^  is  Id-continuous.  We  say  X  is  a 
solution  of  the  matrix  equation  (4.1)  on  T  provided  X  €  D  and  LX{t)  =  0  for  all 
teTl- 
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Next,  put 
Z 


■  X 

PX^ 


onT'' 


and 


S  = 


'u{PT'^Q 

-Q 


0 


on  T. 


Lemma  97.  The  matrix-valued  function  X  solves  (4.1)  if  and  only  if  Z  is  nabla- 
differentiable  and  solves 

Z^  =  S{t)Z. 

Proof.  Suppose  X  solves  (4.1).  We  first  show  that  Z  is  nabla-difFerentiable.  Since 
X  is  a  solution  of  (4.1),  we  have  X  €  D,  and  therefore  PX^  is  nabla-differentiable. 
Furthermore,  X^  is  continuous,  so  by  Corollary  28,  we  have  X  is  nabla-differentiable. 
Hence  Z  is  nabla-differentiable. 

It  remains  to  show  Z'^  =  S{t)Z.  We  again  apply  Corollary  28  which  gives  X^  = 
X^p.  Then 

PPX'^  =  PPX^P  =  {PX^)P  =  PX^  -  u{PX^)^. 

Additionally,  LX  =  0,  and  thus  we  have  (PX^)^  =  -QX.  Substituting  into  the 
above  expression,  we  see  that 

PPX'^  =  PX^  +  uQX. 

Multiplying  on  the  left  by  (F^)“^  then  gives 

XV  =  (pp)-i(PX^)  -b  u{PP)-^QX. 


So 


Z^  = 


■  x^  ■ 

(PX^)V_ 

'{pp)-\px^)  -b  uippy^Qx 

-QX 


u{PP)-'^Q 

-Q 


(pp)-r 

0 


X 

PX^ 


sz. 


Conversely,  suppose  Z  is  nabla-differentiable  and  solves 

Z^  =  S{t)Z. 

Let  Zi  and  Z2  denote  the  components  of  Z.  We  claim  X  =  is  a  solution  of 
(4.1).  We  first  show  that  Zi  6  D.  As  Z  is  nabla-differentiable,  both  Zi  and  Z2  are 
nabla-differentiable,  hence  continuous.  Furthermore,  since  Z^  =  S(t)Z,  we  have 

z7  = 


u{PP)-^QZi  +  {PP)-^Z2 
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=  '  u{PP)-\-Z^)  +  {PP)-^Z2 
=  {PP)-\Z2-UZ^) 

=  {P»)-\Z^2). 


which  is  Id-continuous.  Applying  Theorem  26,  we  see  that  Zi  is  delta-differentiable, 


and 


Z^{a{t))  t  e  Tb 
lims^t+  ZY (5)  teB. 


To  simplify  this,  first  consider  t  eTb-  Recall  that  B  denoted  the  set  of  points  in  T 
which  were  both  right-dense  and  left-scattered.  Hence  for  t  €  T^,  we  have  p{a(t))  =  t, 
and  therefore,  at  these  points. 


Zt{t)  =  Z^(<t(())  =  ((P>')-'Zi){c(t))  =  (F-'ZJW- 


Now,  ifteB,  then  t  is  right-dense  and  left-scattered,  and  we  get 

Zf  (f)  =  lim  Zf  (s)  =  \im  {(P')-'ZS){s)  =  (P-'Z2)(t). 

S— s— 

Since  we  got  the  same  expression  in  both  cases,  we  conclude  that 

ZY  =  P~'^Z2. 

Since  both  P  and  Z2  are  continuous,  Z^  is  continuous.  Furthermore,  PZ^  =  Z2 
is  nabla-differentiable,  and  =  -QZi,  which  is  Id-continuous.  Hence 

Zi  G  1^. 

Finally,  consider  LZi.  We  get 

LZi  =  [pzYT^  +  qz^ 

=  zY  +  Qz, 

=  —QZi  -f-  QZi 

=  0, 

which  completes  the  proof.  □ 

We  have  already  pointed  out  that  S  is  Id-continuous.  In  addition,  it  is  also  v- 
regressive.  To  see  this,  note  that 


7  -  zz2(p^)-ig 

-p{PP)-^' 

I  iy{PP)-^ 

uQ 

I 

-uQ  /-i/2g(pp)-i 

Hence  [I  —  i/S]  is  invertible,  and  therefore  S  is  i^-regressive. 

In  [3]  and  [7,  Section  3.9]  it  was  shown  that  under  these  conditions,  the  initial 
value  problem 


Z^  =  S{t)Z,  Z{to)  =  Zo 
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has  a  unique  solution.  By  the  equivalence  developed  above,  then,  we  have  shown  that 
if  Xo,  Xq  are  any  n  x  n  matrices,  then  the  initial  value  problem 

LX  =  0,  Xih)  =  Xo,  X^(to)  =  X^ 


has  a  unique  solution. 

Definition  98.  The  unique  solution  of  the  IVP 

=  X{to)  =  0,  X^{to)  =  P-\to) 

is  called  the  principal  solution  of  (4.1)  at  to-  The  unique  solution  of  the  IVP 

LX  =  0,  X(to)  =  -/,  X^{to)  =  0 
is  called  the  associated  solution  of  (4.1)  at  to- 

Definition  99.  For  X,  V  €  D,  we  define  the  Wronskian  matrix  of  X  and  Y  by 
W{X,Y){t)  :=  X*{t)P{t)Y^{t)  -  [P(t)X^(t)]*y(t), 


for  t  e  T*". 

We  now  establish  the  Lagrange  Identity. 

Theorem  100  (Lagrange  Identity).  If  X,Y  G  D,  then 

X*{t)LY{t)  -  [LX{t)]*Yit)  =  [W(X,F)]^(t) 

forteT'^. 

Proof.  Let  X,  F  G  D.  Then 

[W(X,F)f  =  {X*PY^  -[PX^YY)"^ 

=  X*[PY^]'^  +  {X*)^{PY^Y  -  ([PX"^]*)^  Y  -  ([PX'^]*)^y^ 

=  X*[PF^]^  +  {X^yP^’Y^P  -  ([PX^]^)*  F  -  ((X^)*P)^  F'^ 

=  X*[PF^]^  +  {X^ypPY^  -  ([PX^]^)*  F  -  (X^)*P^F^ 

=  X’iPF"^]^- ([PX^n^F 
=  X*[PF^]^  +  X*QF  -  ([PX^]^)*  F  -  X*QF 
=  X*  ([PF^^]^  +  QY)  -  ([PX^]^  +  QX)*  F 
=  X*LF  -  [LX]*F 

on  ¥«.  □ 
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Corollary  101  (Abel’s  Formula).  If  X,  Y  are  solutions  of  (4.1)  on  T,  then 

WiX,Y)it)  =  C, 

for  t  €  T'^;  where  C  is  a  constant  matrix. 

Proof.  Suppose  that  X  and  Y  are  solutions  of  (4.1)  on  T.  Then  by  the  Lagrange 
Identity, 

W'^{X,Y)(t)  =  0 

for  all  teT^.  Hence  W(X,  Y)  must  be  a  constant  matrix  on  T”.  □ 

Theorem  102  (Converse  of  Abel’s  Formula).  Suppose  U  is  a  solution  of  (4.1) 
such  that  U (t)  is  invertible  for  a/Z  t  €  T.  IfV  €  D  satisfies 

W{U,V){t)^C 

on  T'^,  where  C  is  a  constant  matrix,  then  V  is  also  a  solution  of  (4.1). 

Proof.  Suppose  Z7  is  a  solution  of  (4.1)  such  that  U{t)  is  invertible  for  alH  €  T,  and 
assume  that  F  G  B  satisfies  W{U,  V){t)  =  C  for  all  t  G  T'*.  Then  by  Theorem  100, 
we  have 

u*{t)LV{t)  -  [LU{t)Yv{t)  =  [wiu,  y)]^(t)  =  =  0. 

But  17  is  a  solution  of  (4.1),  so  LU  =  0.  This  gives 

U*{t)LVit)  =  0. 

Multiplication  on  the  left  by  (C/“^)*  then  gives  the  desired  result.  □ 

Abel’s  formula  shows  that  if  X  is  a  solution  of  (4.1),  then  Vr(X,  X)(t)  =  C  for 
t  G  T'^.  This  leads  to  the  following  definition. 

Definition  103.  (i)  If  A  is  a  solution  of  (4.1)  such  that 

W{X,X){t)  =  0  for  tGT'", 

then  we  say  A  is  a  prepared  solution  (or  conjoined  solution)  of  (4.1). 

(ii)  If  X  and  Y  are  two  conjoined  solutions  such  that 

W{X,Y){t)  =  I  for  tGT'', 

then  we  say  that  X  and  Y  are  normalized  conjoined  bases  of  (4.1). 

Theorem  104.  Suppose  X  is  a  solution  of  (4.1)  on  T.  Then  the  following  are 
equivalent: 

(i)  X  is  a  prepared  solution; 
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(ii)  X*{t)P{t)X^{t)  is  Hermitian  for  all  t  €  T''; 

(in)  X*{to)P{to)X^{to)  is  Hermitian  for  some  to  E  T'^. 

Proof.  First,  we  will  show  that  (i)  (ii).  Assume  A  is  a  prepared  solution  of 

(4.1)  on  T.  Then 

t 

WiX,X){t)  =  X*(t)P{t)X^{t)  -  [P{t)X^{t)]*X{t)  =  0, 
for  t  G  T'^.  This  gives 


x*{t)P(t)x^{t)  =  {x^y{t)p{t)x{t), 

for  t  eV,  and  therefore  X*{t)P{t)X^(t)  is  Hermitian  for  alH  G  T'^. 

Next,  note  that  (ii)  (Hi)  is  trivial.  Therefore,  to  complete  the  proof,  we 

need  only  establish  (Hi)  (i).  Suppose  A  is  a  solution  of  (4.1)  on  T,  and  assume 
there  is  some  to  G  T"'  such  that  X*{to)Pito)X^{to)  is  Hermitian.  Then  we  have 

WiX, X)ito)  =  X*{to)P{to)X^{to)  -  [Pito)X^ito)]*X{to)  =  0. 

Since  Abel’s  formula  tells  us  that  W {X,  X){t)  is  constant,  we  then  have  W {X,  X){t)  = 
0  for  all  t  G  so  A  is  a  prepared  solution.  □ 

Remark  105.  In  the  previous  theorem,  we  assumed  that  X  was  a  prepared  solution  of 
(4.1).  This  assumption  was  necessary  in  order  to  apply  Abel’s  formula  and  establish 
the  implication  (Hi)  (i).  To  establish  (i)  =>  (ii),  however,  we  did  not 

actually  need  to  know  that  X  was  a  solution  of  (4.1),  we  only  needed  A  G  D  satisfies 
W{X,X){t)  =  0  on  T.  Therefore,  in  the  following  corollary,  we  include  only  this 
weaker  assumption,  and  do  not  require  LX  =  0  to  hold. 

Corollary  106.  // A  G  D  satisfies  W{X,X){t)  =  0,  then  the  following  matrices  are 
Hermitian  on  T''. 

(i)  X*PX^, 

(ii)  {Xf’yppx^, 

(Hi)  {X‘^ypx,  and 

(iv)  X*PPXf‘. 

If,  in  addition,  X{t)  is  invertible  for  all  t  E  T**,  then  the  following  matrices  are 
Hermitian  on  T*^. 


(i)  PA-"A-^ 

(ii)  px{xn-\ 
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(in)  PX^X-\ 

(iv)  PPX{XP)-^, 

(v)  PPXPX-\  and 

(vi)  PPX^{XP)-\ 

Proof.  Suppose  X  G  D  satisfies  W{X,X)  =  0.  Then  by  the  same  proof  given  in 
Theorem  104,  X*{t)P{t)X^{t)  is  Hermitian  for  all  t  e  T**.  Evaluating  this  expression 
at  p{t),  we  see  that 

x-Mt))P(p(t))x^(p(t)) 

is  Hermitian.  But  X  €  D,  and  therefore,  X^p  =  X^  which  gives  the  second  result. 
To  get  the  third  result,  note  that 

{xypx  =  (X + px^y  px  =  x*px + pix^ypx, 

which  is  Hermitian  by  earlier  parts  of  this  proof. 

Similarly,  to  get  the  fourth  result,  note  that 

X*PPXP  =  {XP  +  uX^y  PPXP  =  (X^)*P^X^  +  u{X^)*PPXP, 

which  is  Hermitian  by  earlier  parts  of  this  proof. 

Now,  assume  that,  in  addition,  X(i)  is  invertible  for  all  t  €  T''.  Then  by  earlier 
parts  of  this  proof,  we  have  that 

(X<T)*px  =  X*PX‘^.  (4.2) 

Multiplication  by  (X"^)*  on  the  left,  and  X~^  on  the  right  gives 

(X-^)*(X‘")*P  =  PX^X-^ 

which  shows  that  PX'^X"^  is  Hermitian.  Similarly,  multiplying  equation  (4.2)  on  the 
left  by  ((X*^)"^)*,  and  on  the  right  by  (X‘^)“\  we  get 

PX{X‘^)-^  =  ((X‘")-^)*X*P, 

which  shows  that  PX(X‘^)“^  is  Hermitian. 

Furthermore,  by  Theorem  104,  we  have 

x*px^  =  {x^ypx. 

Multiplying  on  the  left  by  (X“^)*,  and  on  the  right  by  X“\  we  get 

px^x-^  =  (X-^)*(X^)*P, 

which  shows  that  PX^X“^  is  Hermitian. 


The  remaining  three  results  are  derived  similarly  from  the  expressions 
{XPypPX^  =  {X^ypPXP  and  X*PPXf’  ^{XPypf’X. 
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□ 

Lemma  107.  The  principal  solution,  X,  of  (4.1)  at  to  and  the  associated  solution 
Y ,  of  (4.1)  at  to  are  normalized  conjoined  bases  of  (4.1). 

Proof.  Note  that 

X*{to)P{to)X^{to)  =  Y*{to)P{to)Y^{to)  =  0, 

which  is  Hermitian.  Thus  by  Theorem  104,  X  and  Y  are  both  prepared  solutions  of 

(4.1) .  It  remains  to  show  that  W{X,Y){t)  =  /.  We  get 

W{X,  Y)ito)  =  X*ito)Pito)Y^{to)  -  [P(to)X^(^o)]*5^(to) 

=  0-[P{to)P-\to)n-I) 

=  1. 

Then  by  Abel’s  formula  we  have  W{X,Y){t)  =  7,  as  desired.  □ 

Theorem  108  (Reduction  of  Order).  Assume  that  U  is  a  prepared  solution  of 

(4.1)  such  that  U{t)  is  invertible  for  all  t  eT,  and  let 

V{t)  :=  U{t)  f  {U*PU’^)-^  (s)  As. 

Jto 

Then  V  is  a  second  solution  of  (4.1).  Moreover,  U  and  V  are  normalized  conjoined 
bases  of  (4.1). 

Proof.  Assume  f/  is  a  prepared  solution  of  (4.1)  such  that  U{t)  is  invertible  for  all 
t  eT.  We  will  first  show  that  V  as  defined  above  is  a  solution  of  (4.1).  By  Theorem 
102,  we  need  only  show  that  F  6  D,  and  W{U,  17)  =  (7  on  T'^,  where  (7  is  a  constant 
matrix. 

In  the  following  calculations,  we  will  suppress  the  argument  t  to  make  things  easier 
to  read,  except  in  the  integral  where  we  include  the  arguments  for  the  sake  of  clarity. 
Consider 


w{u,v)  =  u*pv^  -[pu^yv 

=  U*P  U‘^{U*PU'^)-'^  +  U^  [  {U*PU‘^y'^  (s)  As 

Jto 

-{U^ypu  [  {U*puy~^  (s)  As 

Jto 
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=  7  +  U*PU^  f  {U*PU^)~^  (s)  As 

Jto 

-{u^ypu  f  {u*pu‘^y\s)  as 

Jto 


Now  U  satisfies  the  conditions  of  Corollary  106,  so  {U^)*PU  is  Hermitian,  and  the 
last  two  terms  therefore  cancel.  This  gives 

W{U,V)^I, 

and  we  have  that  the  Wronskian  of  U  and  V  is  a  constant  matrix,  as  desired. 

To  show  that  V  €  D,  consider 

=  U^{U*PUy-'^  +  U^  [  {U*puy~^  (s)  As 

Jto 

=  w^iuy-^p-ywy^  +  u^  f  (t7*P17T^  (s)  As 

Jto 

=  p-\U*)-'^  +  f  {U*PU'^)~^  (s)  As, 

Jto 


which  is  continuous.  Then 

[PV^y  = 


(7/*)-^  +  PU^  j*  {U*PUy-^  (s)  Asj 

=  {{U*y)-^  +  [PU^f  f  {U*PU'^)~^  (s)  As 

Jto 

+[PU^Y  j\u*PU’^)~ys)  As^  . 

The  first  two  terms  axe  clearly  Id-continuous,  but  we  need  to  simplify  the  third  term 
further.  By  Corollary  27,  we  have  that 


PU‘^y^  (s)  As 


V 


r  {u*puy-yp{t))  tejA 

\  lim,^(-  {U*PU‘^y^  (s)  t  €  A. 


Recall  that  A  is  the  set  of  points  which  are  left-dense  and  right-scattered.  Thus  if 
t  e  Tyi,  (T{p{t))  =  t.  This  gives 


V 


PU‘^y^  (s)  As 


u-\a{pmp-Hpmur\p{t)) 

u-\t){pryt){{ury\t) 
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for  t  e  T^.  On  the  other  hand,  if  t  €  ^,  then  t  is  left-dense  and  right-scattered,  and 
we  get 


(s)  As 


V 


lim  (s) 

s— 

u-\t)(prHmur)-'{t), 


for  IgA.  Looking  at  this  again,  we  see  that  in  either  case,  we  get 


r  {U*PU‘')-^  (s)  As 
yto 


and  thus. 


{U*PU^y^  (s)  As 


V 


ppu^pu-\p'’)-\{u*Y)-'^ , 


which  is  Id-continuous.  So,  V  €  D,  and  therefore  by  Theorem  102,  V  is  a  solution  of 
(4.1). 

To  show  that  U  and  V  are  normalized  conjoined  bases,  we  must  show  that  V  is 
a  prepared  solution.  Earlier  calculations  have  already  established  that  W {U,  V)  = 
I.  Recall  that  by  Theorem  104,  F  is  a  prepared  solution  if  and  only  if  V*PV^  is 
Hermitian.  Again  suppressing  the  arguments,  (except  in  the  integral)  we  get 


r  ft  1* 

V*PV^  =  /  {U*PU‘^)~ys)  As  U*P 

Jto  . 


[;^(^<r)-lp-l(^*)-l 

+U^  f  {U*PU‘')~'^  (s)  As 

Jto 

ft  1* 

/  (C7*PC/T^  (s)  As 

Jto 

ft  1*  ft 

/  {U*PU‘^)~'^  (s)  As  U*PU^  /  {WPU^^y^  (s)  As. 

Jto  .  io 


Now,  since  17  is  a  prepared  solution  of  (4.1),  we  have  that  U*PU^  and  U*PU'^  are 
Hermitian.  Therefore  V*PV^  is  Hermitian,  and  thus  F  is  a  prepared  solution  of 
(4.1),  which  completes  the  proof.  □ 
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4.3  Riccati  Techniques 

Next,  we  consider  the  matrix  Riccati  equation 

RZ  =  0,  where  RZ  =  Z'^  +  Q{t)  +  {Zf’)*  [PP{t)  +  (4.3) 

Here,  as  with  the  self-adjoint  matrix  equation,  we  assume  P  and  Q  are  Hermitian 
n  X  n  matrices,  P{t)  is  invertible  for  all  t  e  T,  P  is  continuous,  and  Q  is  Id-continuous. 
We  will  take  the  domain  of  R,  denoted  D/j,  to  be  the  set  of  all  functions  .Z  :  T"  ^  M 
such  that  Z  is  nabla-differentiable,  — >  M  is  Id-continuous  and  such  that 

P^{t)  -h  v{t)ZP{t)  is  invertible  for  any  t  G  T”.  A  function  Z  G  D/j  is  said  to  be  a 
solution  of  (4.3)  on  T**  provided  RZ{t)  =  0  for  all  t  G 

The  following  theorem  shows  us  how  the  self-adjoint  matrix  equation,  (4.1)  and 
the  Riccati  matrix  equation  are  related. 

Theorem  109.  Assume  X  G  D  satisfies  W{X,X)  =  0,  and  that  X{t)  is  invertible 
for  all  t  G  T.  Let  Z  be  defined  by  the  Riccati  substitution 

Z{t)  =  P{t)X^{t)X-Ht),  (4.4) 

for  t  Then  Z  G  Dij  is  Hermitian,  and 

LX{t)  =  [RZ{t)]X{t) 


forteT'^. 

Proof.  Suppose  X  is  as  described  in  the  theorem,  and  define  Z  by  the  Riccati  sub¬ 
stitution,  (4.4).  Then,  by  Corollary  106,  Z  is  Hermitian.  To  see  that  Z  eBr,  note 
that  X  G  D,  which  implies  PX^  and  X~^  are  nabla-differentiable,  so  Z  is  nabla- 
differentiable.  Furthermore, 

Z^{t)  =  [PX^X-^f{t) 

=  [PX^]^(t)(X0-^(t)  -f-  [PX^](t)[X-']^(t) 

=  [PX^]^(t)(X0-'(t)  -  P{t)X^{t)X-'^it)X^{t)iXP)-\t) 

=  [PX^]'^(t){XP)-^it)-P{t)X^it)X-\t)X^f>{t){Xf’)-'^it), 


which  is  Id-continuous.  Next,  note  that 

P<’{t)  +  I/WZ'W  =  P'{t)  +  u(t)lPX‘^X-']i‘{t) 

=  pi’{t)[I  +  u{t)X‘^i’{t){X'’)-'{t)] 

=  P<’(t)|X'’(«)  +  v{t)x'’mx<')-'(t) 

=  P'{t)X{t)(X’‘)-'(t), 
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which  is  invertible,  and  hence  Z  e  Br.  Finally,  to  complete  the  proof,  we  need  to 
show 

LX{t)  =  [RZ{t)]X{t) 

for  t  €  TJJ.  Suppressing  the  arguments,  we  get 

RZ  =  Z^ +  Q  +  {ZP)*[PP  +  I^ZP]-\ZP) 

=  +  Q  +  ZP[Pf’  +  pP{i/X^P){X'’)-'^]-\ZP) 

=  Z^  +  Q  +  ZP[PP[I  +  {X-  XP){XP)-%\ZP) 

=  Z^ +  Q  +  ZP[I  +  X{XP)-^ -I]-\Pf’)-'^{ZP) 

=  Z^ +  Q  +  Z'’[X{X>’)-'^]-\Pf’)-^Z'’ 

=  Z'^ +  Q  +  ZPXPX-\PP)-^Z>’ 

=  [PX^X-'^]'^  +  Q  +  {PX^X-yX^X-'^  {PP)-^  [PX^X-y 
=  [PX^fX-'^  +  [PX^Y{X-y  +  Q 

+pPX^p{Xf’)-^XPX-\pP)-'^P^X^p{XP)-^ 

=  [PX^fX-^  -  [PX^]P[X-^X^{XP)-'^]  +  Q  +  ppx^PX-'^X^P{XP)-'^ 

=  [PX^fX-'^+QXX-^ 

=  [LX]X-\ 

Multiplying  on  the  right  by  X  then  gives  the  desired  result.  □ 

Theorem  110.  The  self-adjoint  matrix  equation,  (4.1),  has  a  prepared  solution,  X, 
such  that  X(t)  is  invertible  for  all  t  €.  T  if  and  only  if  the  Riccati  matrix  equation, 
(4.3),  has  a  Hermitian  solution  Z.  In  this  case,  X  and  Z  are  related  by  the  Riccati 
substitution  (4.4). 

Proof.  First,  suppose  that  X  is  a  prepared  solution  of  (4.1)  such  that  X (t)  is  invertible 
for  all  t  G  T,  and  define  Z{t)  by  the  Riccati  substitution,  (4.4).  Then  by  Theorem 
109,  Z  6  Dii  is  Hermitian,  and 

RZ{t)  =  [LX]{t)X-'^{t)  =  0. 

Conversely,  suppose  the  Riccati  matrix  equation,  (4.3)  has  a  Hermitian  solution 
Z.  Then  —{Pp)~^Zp  is  Id-continous.  Furthermore,  Z  €  D/j,  so 

I  +  u{PP)-^ZP  =  {PP)-^[PP  +  i^ZP] 

is  invertible,  and  hence  —{PP)~^ZP  €  Rv  Then,  by  Lemma  95,  P~^Z  G  R.  Now,  let 
X  be  the  solution  of  the  initial  value  problem 

X^  =  P~\t)Z{t)X,  X{to)  =  I, 

or,  in  other  words,  let  X  =  ep-iz{t,to).  Note  that  although  P~^Z  is  only  defined  on 
T'',  X  is  defined  on  all  of  T.  We  immediately  see  that  X{t)  is  invertible  for  alH  G  T, 
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that  X  is  delta-differentiable,  and  that  is  continuous.  Next,  consider 

=  m)x{t)f 

=  Z^{t)XP{t)-\rZ{t)X^{t) 

=  Z^{t)XP{t)  +  Z{t)X^P{t), 

which  is  Id-continuous  on  T^.  Hence  X  g  D.  Suppressing  the  arguments, 

W{X,X)  =  X*PX^  -[PX^]*X 

=  X*P[P-'^ZX]  -  [PP-^ZX]*X 

=  x*zx-x*zx 

=  0. 

Furthermore,  Z{t)  =  P{t)X^{t)X~^{t).  Hence  the  conditions  of  Theorem  109  are 
met,  and  we  have 

LX{t)  =  [RZ]{t)X-\t)  =  0. 

Therefore  X  is  a  prepared  solution  of  (4.1).  □ 

Next,  we  wish  to  develop  Jacobi’s  condition  for  the  self-adjoint  matrix  equation 
(4.1).  Toward  that  end,  we  want  to  establish  Picone’s  Identity,  which  will  be  crucial. 
The  calculations  required  to  establish  Picone’s  Identity  are  quite  lengthy,  however,  so 
we  break  things  up  by  first  establishing  a  computational  lemma. 

Lemma  111  (Completing  the  Square).  Let  X  be  a  prepared  solution  of  (4.1)  such 
that  X(t)  is  invertible  for  all  t  €  T,  and  define  Z  via  the  Riccati  substitution  (4.4). 
Fix  t  €  T'',  let  u  be  an  n-dimensional  vector  function  which  is  nabla- differentiable  at 
t,  and  let  D  =  X'^X~^{P^)~^ .  Then  at  t,  we  have 

{u*Zuf  =  {u'^ypPu'^  -  u*Qu  -  [Pf’u^  -  Zf’uP]*D[PPu'^  -  Z^u^]. 

Proof.  Let  X,  Z,  and  u  be  as  described  in  the  lemma.  Note  that  Theorem  110  applies 
and  thus  Z  is  Hermitian,  and  is  a  solution  of  (4.3).  Furthermore,  on  T**,  we  have 

pp  +  „ZP  =  PP  +  u[PX^X-Y 
=  pp  +  uPPX^P{Xf>)-'^ 

=  pp[l  +  uX^{XP)~^] 

=  pp  [I  +  {X  -  x^yx^)-^] 

=  pp[i  +  x{Xp)-'^-i] 

=  pPX{XP)-\ 
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Additionally,  note  that  the  Riccati  substitution  tells  us  that 

Z{t)  =  P{t)X^{t)X-\t), 

and  therefore, 

ZP{t)XP(t)  =  pp{t)x^p{t){xp)-\t)xp{t) 

=  pp{t)x'^{t). 

Now  fix  f  €  T  such  that  u  is  nablardifferentiable  at  t.  Then  at  t,  we  have 

{u*Zuf  +  [Pfyy  -  ZPuP]*XfX-\PP)-^[PPvy  -  Z^uf^] 

=  {u*){zuf  +  {v7y{zuy + {ppvy)*x>>x-\pp)-\p^vy) 

-  {ppvyyxpx-\p>’)-^zpup  -  (zpupyxf'x-^ppy^p^u^) 

+  {zpupyxpx-\p>’)-'^{zpu>’) 

=  u*[z^u  +  zPvX]  +  {v7yzf’uf’  +  {v7yp<^xpx-^u'^ 

-  {u'^yppxpx-^ppy^zf’up  -  (upyzpxpx-^'^ 

+  {upyzpxpx-\pp)-^zpup 

=  u*  [-Q  -  zp[pp  +  uzpy^zp]  u + u*zpv7  +  {v7yzpup 
+  {vXyppxpx-'^v7  -  {v7yppxpx-'^ipp)-'^zpup 

-  {upyzpxpx-'^vy  +  {upyzpxpx-\pp)-'^zpup 

=  -u*Qu  -  u*ZPXPX-^{PP)-^ZPu  +  u*ZPy7 
+  iv7yZPuP  +  {vyypp[x  -  iyX^]X-'^vy 

-  {vXypp[x  -  ux'^]x-\pp)-'^zpup  -  {upyzpxpx-^vy 
+  {upyzpxpx-\pp)-'^zpup 

=  -u*Qu  -  u*ZPXPX-'^{PP)-'^ZPu  +  u*ZPv7  +  (vXyZPuP 
+  {u^yppv7  -  u{v7yppx^x-'^vX  -  {u^yzpup 
+  u{v7yppx^x-\pp)-^zpup  -  {upyzpxpx-'^v7 
+  {upyzpxpx-'^{pp)-'^zpup 

=  (v^yppy^  -  u*Qu  -  u*ZPXPX-\PP)-^ZPu  +  u*ZPv7 

-  u{v7yppx'^x-'^v7  +  u{u^yppx'^x-'^{pp)-'^zpup 

-  {upyzpxpx-'^vy  +  {upyzpxpx-\pp)-'^zpup 
=  (vyyppvy  -  u*qu  -  u*zpxpx-\pp)-^zpu + u*zpvy 

-  u{v7yzpxpx-^vy  +  u{vyyzpxpx-\pp)-'^zpup 

-  {upyzpxpx^'^v7  +  {upyzpxpx-\pp)-'^zpup 
=  {vXypPu^  -  u*Qu  -  u*ZPXPX-'^{PP)-^ZPu  +  u*ZPv7 
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-  [uP  +  uv:^YZPXPX-^u^  +  [uP  +  vvy]*ZPXPX-\pP)-^ZPuP 
=  {vFypPvX  -  u*Qu  -  u*ZPXPX-^{PP)-^ZPu  +  u*ZPvX 

-  u*ZPXPX-'^vX  +  u*ZPXPX-'^{pp)-^ZP[u  -  uvX] 

=  {vXypPvX  -  u*Qu  -  u*ZPXPX-'^{PP)-^ZPu  +  vTZPu^ 

-  u*ZPXPX-^vy  +  u*ZPXPX-\pp)-^ZPu  -  vu*ZPXPX-\PP)-^ZPvX 
=  {u^yppy^  -  u*Qu 

+  u*  [ZP  -  ZPXPX-^  -  uZPXPX-\PP)-^ZP]  u^. 

We  now  claim  that  the  last  term  cancels,  which  will  complete  the  proof.  Looking 
only  at  the  center  expression  of  the  last  term,  we  have 

ZP-ZPXPX-'^  -  uzpxpx-'^{pp)-^zp 

=  ppx^pixpy^  -  ppx^p{xp)-^xpx-'^ 

-  v{pPX^P{XP)-^)XPX-^  {PP)-^PPX^P{XP)-^ 

=  pPX'^iXP)-^  -  PPX^X-^  -  jyPPX^X-'^X'^iXP)-^ 

=  PPX^  -  u{X-'^X^iXP)-^] 

=  PPX^  -  X-^[X  -  XP]{XP)-'^] 

=  PPX'^  -  X-^  -[I-  X-^XP]{XP)-^] 

=  PPX^  [(^0”^  - 
=  pPX^iO) 

=  0. 

So,  the  last  term  does,  in  fact  cancel,  which  yields  the  result.  □ 

The  above  lemma  will  be  extremely  useful  in  establishing  the  following  result. 

Theorem  112  (Picone’s  Identity).  Let  a  G  C”  and  assume  X  and  Y  are  normal¬ 
ized  conjoined  bases  of  (4.1)  such  that  X(t)  is  invertible  for  all  t  €  T.  Put 

Z  =  PX^X-^  and  D  =  XP X-'^  {PPy^  on  T^ 

Furthermore,  fix  t  e  T'',  and  assume  u  is  an  n-dimensional  vector-valued  function 
such  that  u  is  nabla-differentiable  at  t.  Then  at  t,  we  have 

{u*Zu  +  a*X-^u  +  n*(X-')*a  -  o*X-^Fa)^  -  {u'^ypPu^  +  u*Qu 
=  -[ppyX  -  ZPyP  -  {{XP)-^ya]*D{PPu'^  -  ZPyP  -  ((X0“^)*a]. 


Proof  We  will  look  at  each  of  the  terms  individually,  then  put  it  all  together  to  get 
the  desired  result. 


First,  consider 


=  -X-^X^{X^)-^YP  +  X-'^Y^ 

=  -X-\P>’)-^PPX'^{XP)-'^YP  +  X-^Y^. 

Now,  by  Corollaxy  106,  P^X^(X^)“^  is  Hermitian,  so  we  see  that 

ppX^iXP)-'^  =  ((X^)-^)*(X'^)*P^ 

Substitution  then  gives 

(j^-iy)V  ^  _x-i(pp)-i((x^)-i)*(X^)*ppyp  +  X-\pP)-\{XP)-^y{XP)*PPY 
=  x-'^{pp)-^{{x>’)-^)*  [-{x^yp^Yp  +  [XpyppY^] 

=  x-^ppyHix^y^y  [x^fpr"^]  -  [px^yY 

=  x-'^{pp)-\{Xp)-'^yw'’{x,Y) 

=  x-\pp)-\{x>’)-'^y. 

Next,  look  at 

{a*X-'^uY  =  a*{X-yuP  +  a*X-\^ 

=  a*X-'^X^{XP)-'^u‘’  +  a*X-'^v7 
=  -a*{X>')-'^XPX-\Pf‘)-'^P^X^{XP)-'^uf’ 

+ a*{xp)-'^xpx-ypp)-'^ppvy 
=  a*{XP)-'^  [DPPy^  -  DPPX^{Xf)-'^uf’] 

=  a^Xpy^  [DPPyy  -  DZPXP{XP)-'^uf>] 

=  a*{XP)-'^D  [pPy7  -  ZPuP]  . 

Now,  note  that  since  P^X(X^)"^  is  Hermitian,  D  =  [PpX{Xp)~^]~'‘-  is  Hermitian 

well.  This  gives 

{u*{X-^yaY  =  ((a*X-^u)^)* 

=  [a*{Xf’)-'^D  [PPu^  -  ZPuP]]* 

=  [pPy^  -  ZPuy  D{{XP)-^y a 

Now,  putting  this  all  together,  and  applying  Lemma  111,  we  see  that 

{u*Zu  +  a*X-^u  +  u*(X-^)*a  -  a^X-^Fa)^  -  +  vTQu 

=  -[PPy^  -  ZPuP]*D[PPvy  -  ZPuP] 

+  {a*X-'^uy  +  {u*{X-yaY  -  a*{X-'^YYa 
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=  - -  ZPuf]*D[PPv7  -  Z^uP] 

+  a*{XP)-^D  [Ppy!^  -  ZPuP]  +  [P>’u^  -  Z^u^]*  D{{XP)-y a 
-a*X-\pP)-\{XP)-^)*a 
=  -[PPyy  -  ZPuP]*D[Pf>v7  -  ZPu»  -  ((X^)-^)*a] 

+  a*{Xf>)-^D  [P<^v7  -  Z>>uP]  -  a\XP)-^XPX-\pP)-^{{XP)-ya 
=  -[P^u^  -  ZPuPYDlPf^u^  -  ZPuP  -  ((X'’)-^)*a] 

+  a*{XP)-^D  [pPv!^  -  ZPuP]  -  a*{XP)-^D{{XP)-^)*ot 
=  -[ppy^  -  ZPuP]*D{PPu^  -  ZPuP  -  ((XO"^)*a] 

+  oc*{XP)-^D\pPu^  -  ZPuP  -  ((XO‘^)*a] 

=  -[PPy^  -  ZPuP  -  {{XP)-'^)*a]*D[PPu^  -  ZPuP{{XP)-^)*a], 

which  completes  the  proof.  □ 

Remark  113.  In  Lemma  111  and  Theorem  112,  we  assumed  that  X  was  a  prepared 
solution  of  (4.1)  and  that  X{t)  was  invertible  for  all  t  €  T.  In  the  work  that  follows, 
we  are  going  to  consider  a  prepared  solution,  X,  of  (4.1)  on  an  interval  [a,  6],  but  we 
are  going  to  allow  the  possibility  that  X(a)  is  singular.  Of  course,  if  X{a)  is  singular, 
then  the  identities  are  no  longer  valid  &t  t  =  a.  A  careful  review  of  the  proofs  shows 
that  there  may  also  be  problems  at  t  =  cr(a).  For  t  €  (o-(a),  b],  however,  the  identities 
hold.  Therefore,  in  the  proofs  of  the  next  results,  we  will  take  care  to  apply  Lemma 
111  and  Theorem  112  only  on  (<7(0),  b]. 

Definition  114.  Let  P  be  the  quadratic  functional  defined  by 


—  u*Qu  (t)  Vt. 


We  say  P  is  positive  definite  provided  P\p\  >  0  for  all  u  6  Cpjj([o,  6],R")  such  that 
u{a)  =  u{b)  =  0,  and  P[v\  =  0  if  and  only  if  u  =  0. 

Here  Cpid([a,  6])  denotes  the  set  of  continuous  n-dimensional  vector- valued  func¬ 
tions  whose  nabla-derivatives  are  piecewise  Id-continuous. 

Definition  115.  A  prepared  solution,  X,  of  (4.1)  is  said  to  have  no  focal  points  on 
(a,  6]  provided 

X  is  invertible  on  (a,  6],  and  XPX~^{PP)~^  >  0  on  (o,  5]. 


We  continue  to  work  towards  Jacobi’s  condition.  The  next  result  establishes  a 
sufficient  condition  for  the  quadratic  functional  P  to  be  positive  definite. 

Theorem  116  (Sufficient  Condition  for  Positive  Definiteness  of  P).  If  there 
exist  normalized  conjoined  bases  X  and  Y  of  (4.1)  such  that  X  has  no  focal  points 
in  (a,  6],  then  P  is  positive  definite. 


89 


Proof.  Let  X  and  Y”  be  as  described  in  the  theorem,  and  let  D  =  on 

(a,  6].  Fix  u  e  Cpij([a,6])  with  u{a)  =  u{b)  =  0. 

We  will  consider  two  cases,  a  right-scattered  and  a  right-dense.  So,  first,  assume 
a  is  right-scattered.  Then,  appl3dng  Lemma  111,  we  get 


=  f  [iu'^ypf>v7  -  u*Qv]  (t)  m 

Ja 

pcr{a)  fh 

=  /  [{u^yppu^  -  u*Qu]  {t)  vt+  [{vyypf’u^  -  u*Qu]  (t)  m 

Ja  J  cr(a) 


[{u^yP^vr  —  u'Qu]{G{a))v{a{a)) 


ph 

-t-  /  [{u*Zuf  -t-  [PPu^  -  ZPuf>]*D[PPu'^  -  ZPuP]]  (t)  Wt 
Ja{a) 

[{vvPypP{vvX)v~^  —  vu*Qu]{cF{a))  -f  {u*Zu){h)  —  {u* Zu){a{a)) 

ph 

/  [PPu^  -  ZPuP]*D[PPu'^  -  ZPuP]{t)  Vt 
Ja(a) 


Recall  that  X  has  no  focal  points  in  (a,  6],  and  therefore  D  >  0  on  (cr(a),6].  Addi¬ 
tionally,  u{b)  =  0,  so  one  of  our  terms  cancels,  and  we  get 


T[u]  >  {{uvyypp{yu^)v-'^  -  uu*Qu  -  {u*Zu)}  (a(a)) 

=  {[u-  uP]*PP[u  -  uP]u-'^  -  uu*QXX-'^u  -  u*Zu)  {a{a)) 

=  ^u*PPuv-^  +  vu*  [PX^]  ^  X-'^u  -  u*PX^X-'^u^  (aia)) 

=  {u*PPuu-^  +  u*[PX^  -  PPX^P]X-'^u  -  u*PX^X-^u}  {a{a)) 

=  {u*PPuu-'^  +  u*PX^X-'^u  -  u*PPX^X-^u  -  u*PX^X-'^u}  {a{a)) 
=  {u*PPuu-^  -  u*PPX'^X-^u}  {a{a)) 

=  {u*pp[uu-^  -  X'^X-^u]}  {a{a)) 

=  {u*PP[I  -  i^X^X-^]uiy-^ }  (cr(a)) 

=  {u*PP[I  -[X-  XP]X-^]uu-'^}  (a(a)) 

=  {u*PP[I  -I  +  XPX-'^]uu-^}  (<T(a)) 

=  {u*PPXPX-^uu-^}  (a(a)) 

=  {u*PPXPX-^(PP)-^PPui^-^}  (a(a)) 

=  {(PPuyn(PPu)}(a(a))iy-\a(a)) 

>  0. 


Clearly,  .F[0]  =  0,  therefore  to  complete  the  first  case,  we  must  show  that  J^[v\  = 


0  w  =  0.  So,  assume  :F[u]  =  0.  By  the  work  we  did  above,  we  have 
:f[u]  =  -  uu*Qu  -  {u*Zu)}  {cr{a)) 


+ 


fb 

/  [PPv7  -  ZPuf‘]*D[P'’vy  -  ZPuf>]  (t)  Vi, 
J(T(a) 


and 


{{i'v7yP'^{uv7)u  ^  —  vu*Qu  —  {u*Zu)}  (cr(a))  >  0. 
Therefore,  since  P\v\  =  0,  we  must  have 

{{vu^ypp{vu^)v-^  -  uu*Qu  -  (u*Zu)}  {a{a))  =  0, 


and 


Thus 


pb 

/  [PPyy  -  ZPu^’YDlPPyy  -  Z^u^] (i)  Vi  =  0. 

J  <T{a) 


pPyP  _  ZPUP  =  0 

on  (<T(a),6].  Hence  on  (cr(a),6],  we  have 

0  =  P^vP  -  ZPuP 
=  ppv7  - 

=  pp[vP  -  x^ix^y^uP]. 

Since  PP  is  invertible,  we  then  get 

0  =  •u'^  -  x'^ixpyup 
=  vy  -  x'^{xp)yu  -  uvF] 

=  u^  -  X"^{XP)-\  +  vX'^{XP)-\'^. 


Rearranging  this,  we  see  that 

X"^{XP)-\  =  u^  +  uX'^{XP)-^u^ 

=  \I  +  uX'^{XP)yu^ 

=  [I  +  [X  -XP]{XPy]u^ 

=  [i+ x{xpy  -  i]vy 
=  x{xpyv7. 


And  finally,  we  have 


vP  =  xpx-'^x^{xpyu 
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on  (cr(a),6].  Now  XPX~^X'^ {X^)  ^  is  Id-continuous,  and 

/  -  vlXPX-^X^iXf)-^]  =  7  -  XPX-'^IX  -  XP]{XP)-^ 

=  7  -  XPX-^X{XP)-^  +  XPX-'^XP{XP)-^ 

=  1-1  +  xpx-'^ 

=  xpx-\ 

which  is  invertible  for  t  G  (cr(o),6].  Thus  XpX~^X'^{Xp)~^  G  TZu,  and  therefore  the 
initial  value  problem 

v7  =  XPX-'^X'^{XP)-\  u{b)  =  0 


has  a  unique  solution  (the  trivial  solution).  Thus  u{t)  =  0  for  t  G  (o' (a),  &]•  Moreover, 
u(a)  =  0.  So,  u{t)  =  0  on  [o,  6]  except  possibly  at  cr{a).  Now  if  cr{a)  is  right-dense, 
then  by  continuity,  u{a{a))  =  0.  The  only  remaining  possibility  is  that  cr{a)  is  right- 
scattered,  and 

u(t)  = 


0  a{a) 
c  t  =  a{a). 


In  this  case. 


u^{t) 


"O  t  ^  {<T(o),cr^(a)} 

^  TO  ^ 

,  i/(ct2(q))  t  =  a^{a). 


But  0-^(0)  G  {a{a),b],  so  we  have 


=  vXia\a))  =  [XPX-^X^{XP)-^]{a^{a))u{a^{a))  =  O! 

and  thus  c  =  u{a{a))  =  0.  This  gives  u{t)  =  0  on  [a,  6],  and  therefore  >  0,  which 
completes  the  proof  in  the  case  where  a  is  right-scattered. 

We  now  turn  our  attention  to  the  second  case.  Assume  a  is  right-dense.  Select 
a  decreasing  sequence,  {0^}“=!)  points  in  T  such  that  limTO-»ooflm  =  cij  and  for 
m  G  N,  let  am  =  -  {iX^)*Pu)  (am)-  Note  that  as  X^,  P,  and  u  are  continuous,  we 
have 

lim  =  -{{X^)*Pu){a)  =  0. 

m-^oo 

Now  applying  Theorem  112  with  a  =  ajn,  we  have 


pb 

/  [{vyyppv7  -  u*Qu]  (t)  vt 

J  dm 

=  f  {u*Zu  +  alX-'^u  +  u*{X-y*am-oclX-'^Yam)"^  Vt 

**  dm. 
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+  /  [PPy!^  -  ZPuf  -  {{XP)-^yamYD[PPv:^  -  Vi 

Jam 

fb 

>  /  {u*Zu  +  a*^X-\  +  u*{X-'^yam-a*^X-'^Yamf  Vt 

Jam 

=  {a*^X-^Yamm  -  (u*Zu  +  a*^X-^u  +  -  a*^X-^Yam)(am) 

=  (a*^X-^YaJ(b) 

-  (u*Zu  -  u*PX^X-^u  -  u*(X-^y(X^yPu  -  u*PX^X-^Y(X^yPu)(aJ 
=  (a*^X-^Yam)(l>) 

-  (u*Zu  -  u*ZXX-^u  -  u*(X-^yX*Zu  -  u*ZXX-^Y(X^yPu)(arn) 

=  (a*^X-^Yam)(b) 

-  (u*Zu  -  u*Zu  -  u*Zu  -  u*ZY{X^yPu){am) 

=  «X-'ya„)(6)  +  {u*Zu  +  u*ZY{X^yPu){am) 


Now  X  and  Y  are  normalized  conjoined  bases  of  (4.1),  and  therefore, 

/  =  X*[PY^]  -  [PX^]*Y 
=  X*[PY^]  -  [ZX]*Y 
=  X*[PY^]  -  X*ZY. 

Rearranging,  and  multiplying  on  the  left  by  (X*)“\  we  get 

ZY  =  PY^  -  {X*)-\ 

We  now  substitute  this  into  the  expression  from  above  to  get 
r  [{u'^ypPv7  -  u*Qu]  (t)  Vi 

Jam 

>  {alX-^Yam){b)  +  {u*Zu  +  u*[PY^  -  (X*)-'](X^)*Pu)(a J 
=  {alX-'^Yam){b)  +  {u*Zu  +  u*PY\X^yPu  -  u*iX*)-^[PX^]*u){am) 
=  {alX-'^Yam){b)  +  {u*Zu  +  u*PY^{X^yPu  -  u*{X*)-^[ZX]*u){am) 
=  [alX-^YamM  +  {u*Zu  +  u*PY^{X^yPu  -  u*Zu){am) 

=  {a*^X-^Yamm  +  {u*PY^{X^yPu){am) 

=  {alX-'^Ya,n){b)  -  {u*PY^am)M 


P[v\  =  lim 

m— »oo 


Jam 


ppy^  -  M*Qu](i)  Vi 


Then 
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>  lim  [{a*^X-'^Ya^){b)  -  {u*PY^ara)iam)] 

m—*oo  *• 

=  0. 


Finally,  suppose  T[u]  =  0.  We  wish  to  establish  that  this  forces  w  :=  — 

Z^u'’  =  0  on  (a,  b].  If  we  do  so,  then  by  the  same  reasoning  used  in  the  case  where  a 
was  right-scattered,  we  will  have  u{t)  =  0,  and  we’ll  be  done. 

So,  let 

Wm  =  P^y^-ZPuf>-{iXf)-^ya^. 


Then  we  have 


lim  Wm  =  w 


uniformly  on  [ofe,  6]  for  each  A;  €  N.  Therefore, 

lim  /  w'^Dwm  Vt  =  /  li 
Jau  Jak 


'a/s,  •'0(5 

Next,  note  that  if  m  >  A:,  then 

oh 


lim  w’^Dwm  Vt 

OO 


Jak 


w*Dw  Vt. 


j  W^^DWm  vt  >  f  wl,DWm  Vt. 

J  dm 


m^^rn 


Letting  m  oo,  we  then  have 


,/ 

J  dm 


lim  /  w'^Dwm  Vt>  I  w*Dw  VI 


/ 

Jak 


But  P[v\  =  0,  and  therefore 

lim  [  wl^Dwm  Vi  =  0, 

J  am 


which  gives 


ph  pb 

=  lim  /  w^Dwm  Vi  >  /  ‘ 
Jam  Jak 


w*Dw  Vi  >  0. 


Finally,  letting  A;  — »■  oo  gives 


/ 


w*Dw  Vi  =  0. 


Then  we  must  have  w  =  0  on  (a,  6],  which  completes  the  proof.  □ 

Definition  117.  We  say  that  the  self-adjoint  matrix  equation  (4.1)  is  disconjugate 
on  [a,  6]  if  the  principal  solution  of  (4.1)  at  a,  X,  has  no  focal  points  in  (a,  b]. 

And  now,  finally,  we  are  ready  to  prove  Jacobi’s  condition. 
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Theorem  118  (Jacobi’s  Condition).  The  self-adjoint  matrix  equation  (4.1)  is 
disconjugate  if  and  only  if  the  quadratic  functional  T  is  positive  definite. 

Proof.  First  suppose  (4.1)  is  disconjugate.  Then  the  principal  solution  X,  and  the 
associated  solution  Y  satisfy  the  conditions  of  Theorem  116,  and  therefore  P"  >  0. 
Now,  suppose  (4.1)  is  not  disconjugate.  Then  there  is  to  £  T  such  that  either 

(i)  to  £  aiid  X  is  invertible  on  (a, to),  but  X{to)  is  singular,  or 

(ii)  to  G  {a,b],  and  X  is  invertible  on  (a,  6],  but 

D(to)  =  (to)  is  not  positive  definite. 


Note  that  exactly  one  of  these  cases  occurs.  Now  define  the  vector  d  €  M”  \  {0}  as 
follows.  If  (i)  occurs,  choose  d  such  that 

X{to)d  =  0. 

If  (ii)  occurs,  then  there  is  a  nonzero  vector  v  such  that 

v*D{to)v  <0. 


Put 

d  =  X-\to){PT\to)v. 

Then  we  have,  at  to 

d*(XpypPXd  =  v*{{P‘’)-^)*{X-yXP*P^XX-\P^)-\ 
=  v*{{PP)-^)*{X-y{Xf^)*v 
=  v*{Xf>X-^{PP)-yv 
=  v*D*v 
=  v*Dv 
<0. 

Next,  put 

fx(t)d  ift<to 

U\i/j  ——  A 

10  otherwise. 

Then  by  the  way  u  is  defined,  we  have 

vX  it)  =  0 

for  t>  to,  and  therefore  for  t  >  to, 

{u'^yP^’u^  —  u*Qu  =  0. 
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Now,  as  X  is  the  principal  solution  of  (4.1),  we  have  X{a)  =  0,  and  therefore 
u{a)  =  u{h)  =  0.  Note  further  that  if  either  to  is  left-scattered  or  to  is  left-dense  and 
(i)  occurs,  then  u  €  Cpi^.  We  will  deal  with  these  possibilities  first.  If  however,  to  is 
left-dense  and  (ii)  occurs,  u  may  not  be  continuous  at  to,  and  we  leave  this  scenario 
for  the  end  of  the  proof. 

So,  assume  first  that  to  is  left-dense  and  (i)  occurs.  Then 


fto 

m  =  / 

Ja 


{{v^yppyy  -  u*Qu)  (t)  Vt 

(d*{X^)*pPX^d  -  (rX*QXd^  (t)  vt 
^d*{X^)*P>’X'^d  +  crX*[PX^fd^  (t)  vt 
(d*[{x*)^{px^y  +  x*{px^)^]d^  (t)  vt 
(d*[X*PX^]'^d'^  (t)  vt 


=  {d*X*PX^d){to)  -  {d*X*PX^d){a) 


=  0, 


0  and  X{a)  =  0.  Next  suppose  to  is  left-scattered.  Then  we  have 
-  u*Qu]  (t)  Vt 

rp{to)  f^o 

=  /  -  u*Qu]  (t)  Vt  -h  /  [{v7)*pf’v7  -  u*Qu]  (t)  Vt 

Ja  J  p{to) 

The  first  integral  simplifies  in  similar  fashion  to  the  case  above  and  we  get 

P[u]  =  {d*X*PX^d){p{to))  -  {d*X*PXH){a)  +  -  u*Q«](to)i/(to) 

=  {d*{XP)*PPX^d){tQ)  +  {{vu^ypp{uvy)u-\tQ) 

=  {d^xpyppx^d)^^)  +  [(«  -  upypp{u  -  up)u-^]{to) 

=  {d*{xf>yp'’x^d){to)  +  [{u>>yppupu-^]{to) 

=  {d*{XpyP'’X^d  +  d*{Xf’yP^X^du-^){to) 

=  {d*{XPyP'’luX^  -  X‘’]du-^){to) 

=  {d*{Xf‘ypPXdu-'^){to)  <  0. 


since  X{to)d  = 

fto 

P[u]  =  / 
Ja 


So,  in  either  of  these  cases,  we  have  P  is  not  positive  definite.  We  now  consider 
the  remaining  case,  (ii)  occurs  and  to  is  left-dense.  Recall  that  we  showed  above  that 


when  (ii)  occurs,  there  is  a  nonzero  vector  d  such  that 

{d*{XP)*P<’Xd){to)  <  0. 
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Then,  since  to  is  left-dense,  we  have 

{d*X*PXd){to)  <  0. 

Now  put 

c  =  X(to)d. 

Since  X{to)  is  invertible,  c  7^  0,  and  furthermore,  c*P{to)c  <  0.  In  fact,  since  P{to) 
is  invertible,  we  have  c*P{to)c  <  0.  We  seek  to  show  that  P  is  not  positive  definite, 
so,  by  way  of  contradiction,  assume  it  is.  First  assume  to  is  right-scattered,  and  let 
^  (®)^o)  be  an  increasing  sequence  of  points  with  limm-»oo  tm  =  to-  Now,  for 
m  €  N,  put 

otherwise. 

Then  since  to  is  right-scattered,  we  have  to  <  b,  so  Umia)  =  Um{b)  =  0.  Also, 
Um  e  Cpj(j,  and  u{to)  =  \/to  -  tmC  7^  0,  so  u  is  nontrivial.  Then  we  see  that 


r(T(to) 

0  <  =  /  {(^m)  (t)  Vt 


fcr{lo) 

+  l  {{«lrP''ul-u^Qu„](t)m 

=  /*°  c*  f  *  )  ?<■(*)  ( -=2=)  c  Vt 

Jtm  V-s/to  —  tm)  \\/to  —  tmj 

=  r  C*  ( - - - )  PP{t)c  Vt  -  r  C*  ( 

Jtm  \t0—tmj  Jtm  \  ^0 

+  K(<^(t))]*-P(to)  [vZitrito))]  I'icrito)) 


Q{t)c  Vt 


We  now  want  to  look  at  what  happens  to  each  of  these  terms  as  m  00.  Looking 
at  the  third  term,  we  see  that 


^m(^(to))  '^m(to)  _  ^  tn^  _  "v/to  tjn  C 

//(to)  /i(to)  K^o) 
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So,  the  third  term,  then,  is 


tp  tm 

Hito) 


P{to)c, 


which  goes  to  0  as  m  — »•  oo.  Next,  we  consider  the  second  term,  the  integral  containing 
Q.  It  is  easy  to  see  that  this  term  also  goes  to  0  as  m  oo.  Finally,  we  turn  our 
attention  to  the  first  term,  the  integral  containing  P'’.  Using  L’Hopital’s  rule  we  get 


P^(t)c  Vt  =  lim 

m~>oo 


£  c>P'>{ty  m 

to  tjTi 


=  c*P{to)c. 


P^{tm)c 

-1 


Thus,  we  have 

0  <  fim  P[ujn]  =  c*P{to)c  <  0, 

m—^oo 

which  is  a  contradiction.  The  remaining  case,  to  right-dense,  is  handled  similarly.  Let 
tie  a  decreasing  sequence  with  limm^oo  Sm  =  to,  then  choose  an  increasing 
sequence,  {tTn}m=i  hnim-»ooim  =  to  and  such  that  for  each  m  €  N,  to-^m  < 

Sm  —  to-  Applying  the  same  technique  as  above  using 


«m(t)  =  S 


t-t 


m=c 


0 


if  t  6  \tm,  to] 

if  t  e  (to,  Sm] 
otherwise 


results  in  a  similar  contradiction. 

Definition  119.  We  call  a  solution,  X,  of  (4.1)  a  basis  whenever 


□ 


rank 


■  X(a) 
P{a)X^{a) 


=  n. 


We  now  get  the  following  version  of  Sturm’s  separation  theorem. 

Theorem  120  (Sturm  Separation  Theorem).  If  there  is  a  conjoined  basis  of 
(4.1)  with  no  focal  points  in  {a,b],  then  equation  (4.1)  is  disconjugate  on  [o,  6). 

Proof  Let  X  be  a  conjoined  basis  of  (4.1)  with  no  focal  points  in  (o,  6).  Since  X  is 
a  basis,  we  have  that  the  matrix 

A  =  X*{a)X{a)  +  {X^y{a)P^{a)X^{a) 


is  invertible.  Now,  let  Y  be  the  solution  of  the  IVP 

LY  =  0,  Y{a)  =  -P{a)X^{a)A-\  Y^{a)  =  p-\a)X{a)A-\ 


98 


Then  by  Abel’s  formula,  we  have 

W{Y,Y){t)  =  W{Y,Y){a) 

=  [Y*PY^  -  {Y^yPY]  {a) 

=  -{A-y{X^y{a)P^{a)p-^a)X{a)A-^ 

+  {A-^yX*{a)p-\a)P\a)X^{a)A-^ 

=  {A-^y  [(X^)*PX  -  X*PX^]  {a)A-^ 

=  0. 

Furthermore, 

W{X,Y){t)  =  W{X,Y){a) 

=  [X*PY^  -  {X^yPY]  (a) 

=  X*{a)P{a)p-'^{a)X(a)A-^  +  {X^y(a)P^{a)X^{a)A-^ 

=  [X*X  +  {X^yP^X^]  {a)A-^ 

=  AA-^ 

=  I. 

Therefore,  we  see  that  X  and  Y  are  normalized  conjoined  bases  of  (4.1).  Then  by 
Theorem  116,  P  is  positive  definite.  Applying  Jacobi’s  condition  then  gives  us  that 
(4.1)  is  disconjugate,  as  desired.  O 

We  conclude  by  establishing  an  analogue  of  Sturm’s  Comparison  Theorem.  Con¬ 
sider  the  equation 

^PX^Y  +  QX  =  0,  (4.5) 

where  P  and  Q  satisfy  the  same  assumptions  as  P  and  Q,  respectively. 

Theorem  121.  Suppose  that  for  all  t  E.T,  we  have 

P{t)  <  P{t)  and  Q{t)  >  Q{t). 

Then  if  (4.5)  is  disconjugate  then  (4.1)  is  also  disconjugate. 

Proof  Suppose  (4.5)  is  disconjugate,  let  u  E  Cpj^  with  u{a)  =  u{b)  =  0,  and  assume 
u  is  nontrivial.  Then  by  Jacobi’s  condition, 

P[u]  :=  t  [(m^)*P%^  -  it)  Vt  >  0. 


Therefore 


r\u] 


f\u]  >  0. 


u*Qu]  {t)  Vt 
(t)  Vt 


Therefore  is  positive  definite,  and  by  Jacobi’s  condition,  (4.1)  is  disconjugate. 
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